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THE INCLINATION TO THE GALACTIC EQUATOR OF THE GENERAL 
MAGNETIC FIELD OF THE GALAXY IN THE SOLAR VICINITY 


G, A, Shain” 


The distribution of the position angles of the planes of preferential vibrations 
of stellar light is discussed (Figs. 1 and 2). The uniform distribution of @' at 1 = 
= 35-65° and Z = 345-355" corresponds to the directions of the spiral arms and leads 
to a more precise determination of the positions of the latter, Some deviations of 
6} from 90° are a result of local fields (Figs. 4, 5 and 6). 


However, a large part of the deviations is systematic and evidently character- 
izes the orientation of the general field of the Galaxy in the solar vicinity. In order 
to determine the latter typical regions with a large number of stars were chosen and 
the mean values of 0' calculated (Table 1 and Fig. 3), These values are well repres- 
ented by a sine curve; therefore the plane of the field is inclined to the galactic 
plane. The angle between the planes is ~ 18°; the planes intersect at 1 ~ 15° and: 
1 ~ 195°, The coordinates of the pole, corresponding to the plane of the magnetic 
field, are 2 © 105°, b ¥ 72 or Q * 12h, 5) ~+ 46°. These coordinates do not coin- 
cide with those of the pole of the local system (Gould's belt), nor with those of any 
other system. The derived plane is not the plane of the magnetic field of the Galaxy 
as a whole, It evidently characterizes the fluctuations in the direction of the field on 
a scale of the order of 1000 parsecs (which is much larger than the scale of the local 
fields), similar to the fluctuations of the spiral arms, 


Immediately following the discovery of polarized stellar light it was found that in low latitudes the plane 
of vibration shows considerable tendency to be parallel to the plane of the galactic equator but that in Scutum 
(1 = 345-355") and in Cygnus (2 = 35-65") this rule is not followed [1], It then was seen that these two directions 
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Fig. 1. Distribution of 6" for stars in the zone/bl < + 6°. 


correspond to longitudes where the line of sight is directed approximately along spiral arms (Sagitarrius and 
Cygnus), and the distribution of angles @" must be more or less a random one. However, in addition it can 
sometimes be seen in some of Hiltner's figures [2] that @' deviates systematically from 90°, An investigation 
of the slope of the plane of vibration in considerable portions of the sky and of large scale nonuniforimities in 
this slope is of interest for study of the magnetic field structure of the Galaxy, the relationship between diffuse 
matter and the magnetic field and, finally, the origin of the interstellar magnetic field. 


100" 


Fig. 2. Distribution of 6" for stars outside of the zone b= + 8°. 


The first question which can arise is whether (within the limits of observational error) the pole which is 
determined from the magnetic field of the Galaxy coincides with the pole that is determined from the stars, 
The answer will obviously be of importance in connection with the foregoing problems. For this purpose we shall 
consider the distribution of position angles @" which characterize the orientation of the plane of electrical 
vibration of stellar light with respect to the galactic equator. For study of the longitudinal distribution of polari- 
zation it is convenient, as usual, to employ diagrams in which the vectors (P, @') are plotted. To represent 
the general picture of distribution of @' for stars with latitudes 4 8° over 3/4 of the galactic zone we were 
able to compare only the values of @’ with longitude. In Figure 1 each point represents a star with the coordi- 
nates galactic longitude and @'. For this purpose stars were taken from both of Hiltner's catalogues [3], from 
Hall and Mikesell’s catalogue [4] if not included in [3] and stars of the southern sky from Smith's list in Refer- 
ence[5]if not included in References[3]and [4], Stars for which P = 0.3% were not included, In addition, Fig. 2 
represents stars of latitudes b = 8° with the previously known random distribution of the position angles @’. 


The nonuniform longitudinal distribution of the points in Figure 1 results simply from the real nonuniform 
longitudinal distribution of hot stars and also possibly from selectivity in the observational data. Of greater 
interest is the distribution with respect to the ordinate @' which evidently is entirely independent of the distri- 
bution of stars in the sky or of selectivity inthe data. Figure 1 clearly reveals two peculiarities. 


1) There is more or less uniform distribution of the angles @' in the interval 0-180° for the longitude 
intervals 327-347° and 40-50°. In these cases the line of sight appears to be approximately along the spiral 
arms of the Galaxy. The star density here is considerable and it can be stated that polarization observations 
definitely reveal spiral arms in Cygnus and Sagittarius. The very marked reduction in the number of stars on 
both sides of these two directions is another independent argument for the reality of the indicated spiral arms. 


2) Inthe majority of regions besides those 
which correspond to spiral arms in Cygnus and 
Sagittarius the points are concentrated around 0 = 
ce OE ee = 90°; this is the known tendency of the plane of 
. vibration to be parallel to the galactic equator. 
However, for the larger part there are systematic 
shifts to one side or the other. By taking only charac- 


sl teristic regions with a considerable number of stars 
“| (for which an average value of @" is justified) we 
- obtain the table and Figure 3. 
Figure 1 and the average values of 9" provide 
Fig. 3. Average value of @' for star groups as a a serious basis for suspecting that the deviations 
function of J @'-— 90° can be represented roughly by a sine curve 
on which are superimposed the individual deviations 
caused by local smaller-scale fields. We can now 
- confidently speak of the existence of such fields in 
! | 0’ = | Constellation low latitudes in addition to the general magnetic 
field of the Galaxy. For example, such local fields 
Sa can be observed in the "Coal Sack” (Fig. 4) or at 
180—170°)  108°2 25 | Mon | 20h 5™ + 35° (Fig. 5), because in both instances we 
enc Sen 7. cones observe a definite correspondence between the orien- 
110—90 89.3 274 | Pers-Cass tation of the field and the structure or size of the 
9 ¢ = = 
gong | $5.9 | 908 | Pom ceph-coms neta 
ae oy MT Wy (anes) Although the average values of 6" in Figure 3 
347-—327 89.3 112 | Sagittar-Sct depend to a certain extent on the individual groups 
301279 88.0 29 ld tle (represented in the table) any reasonable division into 
274—270 | 86.6 + 1°7 48 cus groups results in a picture which is more or less like 
64.8+3.1) 12 If Figure 3. A glance at Figure 1 (except for sparse 
300—270 87.4 40 | excluding the 


stellar regions) shows that Figure 3 correctly repres- 
ents the data. Comparison with photographs of the 
sky enables us in some instances to distinguish indi- 
vidual local fields superimposed on the general galac- 
tic field. For example, in the region of the Coal 
Sack we can clearly observe both the general galactic 
field and the local field which is in the direction, 
roughly speaking, of the apex of the "bay" to the SE of a Crucis (Fig. 4). For the galactic and local fields the 
values of @' are, respectively, 86.6° + 1.7° (18 stars) and 64.8° + 3.1° (12 stars). 


local field 
around & Crucis 


There is a certain indeterminacy at the other end of the observed region, around 190-210° (Fig. 6). On 
the one hand, we have here the close group of stars NGC 2323, for which P ~ 1.2% and @ “avo 60° (top of 
Fig. 6). On the other hand, the presence in the photograph of the very elongated emission nebula $172 and of 
a few dark nebulae elongated in approximately the same direction (at a very large angle ~ 50° to the galactic | 
equator) apparently indicate local fields. This hypothesis is supported by the orientation in the same direction 
of the plane of light vibration for two stars a little to the west of the emission nebula. It is also possible that in 
the region of longitudes 190-210* the line of sight is approximately along a spiral arm (opposite to the direction 
of the arm in Cygnus). When all of these circumstances and especially the small number of stars are taken into 
account, it is best not to consider this region in connection with the study of the general galactic field. 


It can be seen from Figure 4 that the plane of the general galactic field intersects the plane of the galactic 
equator at 18-20°; the points of intersection (the maximum and minimum of 6 ")have JZ *15° andl 195°, 
so that the galactic coordinates of the pole associated with the magnetic field are approximately Z = 105° and 
b = 72°, or in equatorial coordinates ag 12h 59 * + 46°. It is of interest that these coordinates do not coin- 


cide with the pole of the local system (the Gould belt) of B stars, nor with other systems such as the system of 
helium stars. 
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Fig. 4. General field (18 stars) and local field (12 stars) in the Coal Sack. 


Fig. 5. The region around 205m 4 35°, 


Fig. 6. The region around? = 190-200". 


The coordinates given above for the pole of the magnetic field do not, of course, pertain to the pole of 
the entire Galaxy. If the general regular field of the Galaxy were to agree fully with the indicated plane the 
center of the Galaxy would be located a few thousand parsecs from the plane of the field, which is very improbable. 
If the regular field of the Galaxy is associated with gas clouds it should coincide approximately with the piane of 
the Galaxy, as shown by 21-centimeter observations. It is most likely that the inclination of the plane of the 
field shows the existence of fluctuations of the general field, on the scale of about 1000 parsecs, because the 
separations of most clouds, which through absorption cause the polarization of light from the observed stars, do 
not exceed a few hundred parsecs. The local fields are still smaller and evidently do not exceed a few tens of 


parsecs. 
The investigations of Dutch astronomers have shown that a spiral arm determined by 21 cm radiation does 
not always lie in the plane of the Galaxy but sometimes forms an angle of as much as 40°. Fluctuations in the 
structure of the spiral arms is also observed in extra-galactic nebulae. In the present instance we may be deal- 
ing with a phenomenon of the same type. It would be interesting to obtain the distribution of gas clouds in the 
vicinity of the sun for comparison with the plane of the magnetic field which was determined above. The cloud 


distribution can be obtained by comparing data on interstellar lines, ionized hydrogen regions and 21-cm 


intensities. 


Crimean Astronomical Observatory, 


Academy of Sciences, USSR Received August 27, 1955 
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THE MORPHOLOGY OF GALAXIES, III. DIFFUSE MATTER IN 
SPHERICAL STELLAR SYSTEMS 


B. A. Vorontsov-Vel'iaminov 


Analysis of published data shows that diffuse matter is not at all uncommon 
in spherical stellar systems consisting of type II populations. It often occurs in the 
form of dust and, particularly in elliptical galaxies,in gaseous form, and is always 
present in nuclei of spirals. Observations do not exclude the possibility that diffuse 
matter may also be present in globular clusters. It is shown that double galaxies 
similar to M 51, in which the spiral arm of one of the galaxies forms a “bridge” 
joining it to a smaller, elliptical galaxy, occur fairly often. Elliptical compan- 
ions of M 51— NGC 5195 and M 31 — NGC 205 contain diffuse matter in rougly 
the same proportions as the spiral arms of our own galaxy. Evidently, both of 
these galaxies are in the initial stages of transformation into spirals. It is prob- 
able that spirals evolve from elliptical galaxies as a result of outflow of diffuse 
matter from their nuclei, if there is a sufficient amount of it. In spherical stellar 
systems this matterdoes not appear to be due to accumulation or accretion over 
a long period of time, but is present from the very start. 


1. It is a common belief that spherical stellar systems do not contain diffuse matter, the latter being 
characteristic of plane systems only. In this connection it is often held that the evolution of plane systems from 
diffuse matter is probable, but nothing is said about the origin of spherical systems. 


This belief is widely held in spite of indications that diffuse matter does occur in a numter of spherical 
systems, 
Below, we shall try to summarize and discuss what is known about the presence of diffuse matter in spheri- 


cal systems such as globular’ stellar clusters, elliptical galaxies, and nuclei of spirals. Diffuse matter is taken to 
be present when its density is appreciably higher than that in the surrounding metagalactic space. 
) 


2. Globular stellar clusters. There are no indications, so far, of the presence of dust in globular clusters. 
It would be difficult to detect dust having a density of the order of 10-* g/cc (as in the solar vicinity) in clusters 
whose dimensions are a few tens of parsecs, but denser clouds, similar to Coal Sacks, could be detected. Thus, 
although one cannot be completely sure that dust is entirely absent in globular clusters, one can be sure that if 
it is present it must be distributed just as regularly as the different types of stars in these clusters. 


The presence of gas of density of 10g g/cc, or thereabouts, cannot be detected from spectral absorption 
lines either, because of the small extent of the clusters. 
It would be easier to detect gas by its luminescence, but for this to occur the presence of hot stars is neces- 


sary. A star of type 08 has been found in M3 [1], and a star of type B5V in M13. 


However hydrogen luminescence can just barely be detected in clusters, It would be brighter if denser 
accumulations were present in clusters, e. g., gaseous nebulae, which should be close to the single hot stars found 


in M3, The probability of this is smaJ1 and hot stars are apparently rare exceptions in clusters, Thus the absence 
of observations of optical absorption and emission by (ionized) gases in globular clusters cannot be considered 


as evidence for the absence of gases in them. 


If hydrogen were present in clusters, it would be in the neutral state and hence could not be detected from 
continuous radiospectra which are due to ionized hydrogen, It is true that continuous radio emission is observed 
over spherical regions of such galaxies as our own and M31. However, phenomenologically, this emission, em- 
bracing a region of space comparable with that occupied by the spherical stellar component of our galaxy, does 
not have the same kind of concentration towards the center as the spherical component itself. Thus the matter 
which gives rise to the spherical component of radio emission does not have a direct connection with the stars 
of the spherical component. Consequently, this type of radio emission may be absent in globular clusters because 
it is not due to stars of the spherical component, but to something else. 


If the spherical component of radio emission is due to the bremsstrahlung from relativistic electrons in a 
magnetic field then in stellar clusters there may not be sufficient sources of such electrons and, in any case, 
the presence of a sufficiently strong magnetic field is fairly improbable. For this reason (and quite apart from 
the small dimensions of clusters and the probably small total emission even in most favorable conditions) one 
would not expect to observe appreciable continuous radio emission from clusters. 


The negative result obtained by Mills [2] for globular clusters at a wavelength of 3.5 m must be considered 
in the light of the above considerations. . 


There remains the problem of the possibility of detection of neutral hydrogen at 21 cm. I. S. Shklovskii 
[3] shows that, if the diameter of a globular cluster is 100 parsecs, and the density of gas in it is 10cm Si Cas 
then at 100° K the radiation from it would be the same as that along the direction towards the galactic center, 
However, such data are not very useful because it is unlikely that the density of gas, right up to the cluster bound- 
aries, is higher by one order than the density in the spiral arms of the galaxy. Furthemore, the diameters of the 
largest known clusters are only 60 parsecs, and the stars near their boundaries are very rarefied. If the gas den- 
sity in the central region of a cluster having a diameter of about 20 parsecs is much less than in the spiral arms 
of the galaxy, the expected radiation at 21 cm should be weaker by four or more orders of magnitude than that 
in the direction towards the galactic center and cannot, so far, be detected. 


Summing up all the above, we conclude that the existing and possible kinds of observations do not, and 
cannot, indicate the presence of diffuse matter in globular clusters, but, at the same time, they cannot exclude 
the possibility that in the dense parts of these clusters there may be gas and dust having densities comparable with 


those in the galactic plane. It may be assumed that if such matter is present there are no pronounced localiza- 
tions in it, 


3. Elliptical galaxies. Rarefied spherical galaxies of the Sculptor type have not as yet been extensively 
studied, and we shall not discuss them here. The spherical galaxies known earlier differ from globular clusters 
in dimensions and mass only, have the same degree of concentration towards the center, and consist of the same 
stars (including some known bright stars), Elliptical galaxies are redder than globular clusters, and their ratio 
of mass to luminosity is higher. It is probable that this difference is due to differences in luminosity functions 
for weak stars, It may, however, be admitted that the increase in the mass — luminosity ratio in elliptical galaxies 


may be partly due to a greater density of nonluminous diffuse matter. In fact, there are definite indications that 
diffuse matter is present in elliptical galaxies. 


To begin with, according to Mayall [4], the emission line \ 3727 A (OII) was observed in central parts of 
20% galaxies of type E. Actually, the percentage of galaxies of type E which contain gas is much higher. It 
is probable that this spectral region was not sufficiently well explored in all the galaxies which have been studied. 
According to Page's data [5], among 16 elliptical galaxies which he studied, the line A 3727 A was observed in 
five cases, and, in addition, in three of these, bright Hq and intense A 6584 A (NII) were also noted. The same 
but weaker lines were observed in galaxies which did not give rise to the line 43727 A. Thus among the sixteen 
type E galaxies, the line A 3727 A was observed in 31%, and the lines Hq and NII alone in a further 31% of these 
galaxies. Altogether gas emission lines were discovered in 62% of galaxies of type E. Although the data are not 
pey extensive, the large dispersion employed by Page (which weakened the continuous spectrum but made it 
easier to see the emission lines), and the fact that he covered the red part of the spectrum as well, leads us to 


10 


suppose that the figure of 62% is much closer to reality than the 20% obtained earlier from observations on the line 
A3727 A alone. Thus one may say that gas emission lines are observed in the nuclei of the majority of elliptical 
galaxies of the usual type. 


In our own galaxy the reason for the luminescence of interstellar gas is its irradiation by hot stars, except 
for rare exceptions. According to Baade [6], about a dozen blue giants are known in elliptical companions of 
M 31: NGC 205 and NGC 185, 


Nothing is known about the presence of hot giants in other elliptical galaxies, since these only appear on 
special UV photographs. By analogy with globular clusters which only rarely contain hot stars, and the small 
number of these stars in NGC 205 and NGC 185, it may be supposed that in 38% of type E galaxies without bright 
lines, the latter are not visible not because the galaxies do not contain gas, but because there is not a sufficient 
number of hot stars to produce luminescence of the gas. Thus it is probable that gas is present in all elliptical 
galaxies in appreciable amounts. Gas emission is very considerable in the spherical radio- galaxies NGC 4486 
and NGC 5128, 


Unfortunately, it is impossible to estimate the amount of gas in type E galaxies since no information is 
available on the intensity of lines in their spectra, the optical thickness of the gas within them, or whether the 
gas is luminous or not. 


Less is known about the presence of dust in elliptical galaxies, If the dust does not show pronounced locali- 
zations, it is very difficult to establish its presence. Dust localizations may even be masked by a thick layer of 
stars in front of them. Dust clouds come out better on UV photographs, but these have only been taken in the 
case of the companions of M31, which is rather convenient since of the four elliptical colmanion: of this gigan- 
tic spiral two contain large dust clouds. 


Secondly, due to the brightness of the companion of M31 — NGC 205 of type E5p, on large scale UV photo- 
graphs [6] of its nucleus, a circular dust nebula can be clearly seen due to its strong absorption, its diameter 
being 6 parsecs and its distance from the center 100 parsecs on the new distance scale (450,000 parsecs to the 
group of galaxies M31), The form of galaxies strongly depends on exposure. Thus in the second photograph of 
[7] the above dark nebula does not appear in the overexposed middle part of the galaxy. Instead, patches of less 
dense nebulae on the other (NE) side of the galaxy have merged into a large dark area of 40 x 170 parsecs in 
size. The above photograph gives the impression that a large number of adjoining clouds, a few parsecs in dia~- 
meter, surround the nucleus. They are mainly on the periphery of the central part of this galaxy, mainly in the 
SW part. The major axis of this part we estimate to be 250 to 400 parsecs. 


As can be seen from the long exposure photograph (four hours) given in [7], one end of this central ellipse 
of NGC 205 is cut off along a chord by a dark nebula 40 x 100 parsecs in size. Dark nebulae in NGC 205 are 
quite visible in the photograph given in the nebular atlas of V. G. Fesenkov and D. A. Rozhkovskii. 


The absence of photometric data makes it difficult to make an estimate of the absorption in these clouds 
or of their mass. Furthermore, lack of knowledge of the exact position of these clouds within the galaxy would 
still make this estimate difficult even if the photometric data were available. 


Assuming that the thickness of the clouds is equal to one half of their mean visible dimensions, and the 
density is 10°™ g/cc, we obtain 10° Mo for their total mass, without taking into account the more rarefied dust 
in the remaining parts of the galaxy. If, as in the case of the clouds in T Taurus, the amount of gas is 100 times 
that of dust [8], then the total mass of diffuse matter in NGC 205 is found to be of the order of 107 Mo, while 
10°: ‘Mo is quite possible. Since the mass of this galaxy is probably of the order of 108 M3 we see that the frac- 
tion of the total mass taken up by diffuse matter and the stars is similar to that in the spiral arms of our own galaxy. 


Baade [6] has discovered about a dozen blue stars of type B and absolute magnitude of 4 (new scale) in 
UV photographs of NGC 205, and he ascribes the exceptional nature of this fact to the rarity of dust in this area. 
An examination of this photograph shows that the hot stars are on the periphery of the central region and group 
themselves near the main dark nebula, a fact which has not thus far been noted. 


According to Baade, NGC 205 differs from the usual dense elliptical galaxies by a smaller concentration of 
stars towards the center, and appears to be a system in transition to a highly rarefied elliptical system of 


the Sculptor type. On red sensitive photographs it decomposes into stars right into the center. From counts of 


bright stars Baade [9] has found a major axis of 9' in NGC 205, Reynolds [10], has obtained a major axis equal 
to 12" on the basis of photometric studies of photographs taken with an exposure of 30, 


In this connection, the photograph of NGC 205 taken with long exposure and published in [7] is of parti- 
cular interest. 


In this photogrphs the major axis of NGC 205 also has a length of about 12", i. e., about 1500 parsecs. 
The central ellipse, rapidly weakening towards the edges and visible on other photographs, appears to be only the 
nucleus.of the galaxy. The ratio of its semiaxes is 1, and hence NGC 205 was classified as E5. However, in 
the photograph taken with a long exposure, the nuclear ellipse lies within a weak ellipse in such a way that the 
major axes do not coincide but form an angle of 30°, The ratio of the axes of the outer ellipse is 1:3.4 which 
would correspond to E7. However, in a high contrast reproduction made by us of the photograph given in [7] 
it is clear that the ends of the outer nebula bend in opposite directions, very much like wide spiral arms. It 
turned out that the author of the photograph (Pease) assumed that NGC 205 is a spiral. Thus, in the classification 
of spiral galaxies, NGC 205 appears to be in the initial stages of evolution of a spiral galaxy, with the character- 
istic bright ellipsoidal nucleus relatively larger than in galaxies of type Sa, At the same time, this is the 
smallest known spiral although on reproductions of photographs taken with a 48" Schmidt objective its length 
reaches 18’, i. e., 2200 parsecs. The structure of NGC 205 is not visible in these photographs. Instead, one gets 
the impression that it is joined to M31 by a luminous bridge similar to the bridges discovered by Zwicky between 
some multiple galaxies. The bridge is directed towards the nucleus of M31 and has a width which is somewhat 
smaller than the minor axis of NGC 205, This impression must be verified with the original negative. It is 
remarkable that, as was pointed out to us by F. A. Tsitsin, counts of giants carried out by Seifert and Nassau [11] 
show an offshoot curving outwards, and towards NGC 205, in the outer equal-density curves surrounding this galaxy. 
This offshoot is much wider than the bridge we expected. The nature of the offshoot and its well defined connec- 
tion with NGC 205 excludes the possibility that one is dealing here with an accidental bend in the equal density 
curves of stars belonging to M31. Judging from large scale photographs, spiral arms of NGC 205 do not contain 
giants, but the inner parts of the spiral arms of M31, and other galaxies of type Sb, do not contain them either. 


The lack of coincidence between axes of the visible nucleus and the spiral system of NGC 205 is character- 
istic of a spiral seen at an angle, The following facts bring NGC 205 close to a spiral: 1) spiral arms, 2) the 
presence of a nucleus, 3) presence of dust on the nuclear periphery, 4) presence of hot giants although small 
in number. The presence of these giants within the nucleus is anomalous. 


At a distance of 7° from M31 there is a pair of weak elliptical galaxies, NGC 147 and NGC 185, which form 
a physical system at a visible distance of 58' and also appear to be companions of M31. They both contain a 
number of globular clusters and resolve into stars. NGC 147 has an extremely small concentration towards the 
center,and, according to Baade [6], it is free from dust since distant galaxies are visible through it. NGC 185 
also has a small concentration towards the center, the latter being intermediate between that in NGC 205 and 
NGC 147. For this reason they were classified as Ep by Hubble, 


In NGC 185, at a distance of 50 parsecs from the center, and still in the dense part, one can see a very 
dense dust cloud stretched out along the periphery and 8 x 30 parsecs in size. According to Baade [6], there 
are about a doxen blue giants here, with absolute magnitudes between — 3 and —4. 


Whether there are any emission lines in the spectra of NGC 185 and NGC 205 we do not know. 


After the radio- galaxies NGC 5128 and NGC 1316, and among systems consisting of type II population, the 
elliptical galaxy NGC 5195, a companion of M51 (Fig. 1), is the richest in dark matter. We note, by the way, 
that NGC 5195 cannot be an isolated galaxy on which a spiral arm of M51 is accidentally projected as was sugges- 
ted by Holmberg [12]. It is true that one can see, beyond the edge of the galaxy-companion, a weak chain of 
hot stars A, which is a continuation of the spiral arm of M51, but this is very short. The spiral arm practically 
ends where it visibly "sticks into" the companion. The spiral arm of M51 is, at the same time, of the type ob- 
served by Zwicky in some multiple galaxies. These "bridges" were shown by Zwicky (from their color only) to 
consist of hot stars just as in the case of the spiral arm joining M51 and NGC 5195. We have discovered an even 
more remarkable fact, namely, that such pairs of galaxies are not rare. We found an essentially similar picture 
for NGC 4485-90. An elliptical companion (also somewhat deformed) is joined by a spiral arm to the main 
component. Here the main component is of the late Sc type and its equatorial plane is turned towards us much 


more than M51. It is possible that NGC 17678 is a similar system judging by the drawing given by Rosse [13]. 

In the drawing of this spiral galaxy one wide arm ends in a bright circular localization having a diameter equal 
to the width of the arm. That this is the case is shown in a reproduction of a photograph of NGC 5278-9, found 
by us later, which has an essentially similar form, NGC 5195, considered by some regular and by others irre- 
gular, has the form of an ellipse distorted by two, almost diametrically opposite, regions of large masses of 
dark matter. In addition, on one side one can see a weak, 
and reducing to zero, wide continuation of the luminescence 
of this galaxy. “Analytical photographs" by Zwicky [14] 
show fairly well that the main mass of stars of NGC 5195 is 
of the same nature as the nucleus of M51 itself, i. e., it 
consists of stars of type Il population. Dense areas character- 
istic of spiral arms and formed by concentrations of giants 
and gaseous nebulae are not visible in the galaxy-companion. 


~-o-ccrrn - 


On the usual "blue" photographs one can see dark matter 
on the background of the companion which, since it is placed 
in opposite regions, weakens the radiation emitted by the 
companion so that the latter appears to be stretched out along 
the line joining its center to the nucleus of M51. If the 
"blue" negative is superimposed on the "yellow" positive of 
Zwicky [14], the yellowish areas disappear and the dust is 
even more Clearly visible than in the "blue" photograph. 
Here the companion looks almost like a long rectangle (with 
torn edges because of the inclusion of dark matter) stretched 
along the radius of M51. Both from this picture and the 
"blue" picture it is quite clear that the end of the spiral arm 
which leans against the companion consists mainly of dust, 
with a small number of stars with moderate liminosity, and 

a very small number of supergiants and their concentrations. The width of this arm is just equal to the diameter 
of the companion. In the opposite part of the companion one can see a dark patchy nebula a with a protruding 
area and fibers b, c and d (Fig. 1). The sharp ends of these fibers point away from the nucleus of NGC 5195 
and almost at right angles to the spiral arm. For this reason these fibers could hardly be a continuation of the 
spiral arm. Also, on the opposite side of the companion, dust is more dense than in the spiral arm and hence a 
part of it may belong to the companion itself and not to the spiral arm as such. 


Fig. 1. Schematic drawing of NGC 5195. 
Dust is shown by the dark patches. 


An unusally interesting feature of reproductions of M51 given in [15] and [16] is a long, thin, and curved 
canal, 15 parsecs thick (assuming the presently accepted distance of 2.5 million parsecs) which curves round the 
central dense region of the nucleus of NGC 5195. It passes the center of NGC 5195 at a distance of only 50 par- 
secs (the diameter of this galaxy is not less than 1500 parsecs), Its length is 250 parsecs and it is perpendicular 
to the spiral arm. The fact that it curves round the dense patch in the nucleus clearly shows that it belongs to 
the elliptical galaxy and not the spiral arm of the main component. 


On the "yellow" photograph, and the "yellow" negative combined with the "blue" positive, the compan- 
ion has the form of an ellipse, the major axis of which is perpendicular to the spiral arm and the extension of the 
companion in " blue" photographs. In the compound photograph one can see, in addition, a wide continuation E 
of the luminescence of the companion which bends in a way which suggests the beginning of a wide spiral arm. 
The wedge-shaped fibers cut into.the beginning of this would-be spiral arm, At the diametrically opposite end 
of the compound photograph one can see a similar continuation of luminescence but much shorter and bending 


in the opposite direction. 


The elliptical galaxy rich in dust, NGC 5195, which consists of stars of type II population, and which can- 
not be called irregular only because its outline is broken up by dark nebulae, contains gas. In its spectrum [5], 
the emission line H, can be seen but the lines NII and 3727 A are absent, and the lines Hg — Hg are seen in 
absorption. It is probable that in this integral spectrum of the galaxy, stellar absorption lines mask the always 
weaker gas emission lines Hg — Hs. One may expect that these emission lines are due to a small number of 
hot stars in the central part of NGC 5195. NGC 5195 is evidently similar to NGC 205 in that it contains a con- 
siderable amount of dust and shows signs of formation of spiral arms, Both of these galaxies are interesting in 
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w type of galaxy which is transitory between elliptical and spiral, and is related aD the 
The presence of dust in them is very significant. One may assume that if ABE is 
present in an elliptical galaxy (undoubtedly together with cool gas), or when they appear, the formation : eye 
arms begins. The degree of development (length and strength of arms) probably depends on the amount of diffuse 
matter within the nucleus of the galaxy. The amount of diffuse matter in the nucleus may be epee e to be 
proportional to the volume of the nucleus. This does not explain the striking fact that in spite of the pliterenke 
in absolute dimensions of galaxies, the length of arms in degrees is roughly the same. For example, M31 is 

eight times as large as M33, but the length of arms in degrees is not very different. 


that they represent a ne 
nuclei of spiral galaxies. 


The presence of dust in elliptical galaxies and the nuclei of spirals should be shown in UV photographs 
(in which absorption is more marked), and with suitable exposure, In the usual photographs these objects in the 
central parts are usually overexposed, and this masks the dust, which even without this is masked by a thick 


layer of stars in front. 


In the absence of special observations it is difficult to judge the degree of diffusion of dust in elliptical 
galaxies. Apart from the three type E galaxies described above which contain dust, two of which appear to be 
revolving spirals, two radio- galaxies of type E with strong bands of dust are known. They are NGC 5128 and 
NGC 1316. We consider that these radio-galaxies are not colliding galaxies but irregular galaxies of type E, 
rich in dust and gas, and we shall discuss them in another paper. 


4, Nuclei of spiral galaxies. It is known that these do not differ in form from globular clusters and 
elliptical galaxies, and are comparable to the latter in size. The brighter stars in the systems of all three types 
are evidently the same both in luminosity and spectrum. 


However, nuclei of spiral galaxies differ from globular clusters and elliptical galaxies in having a central 
plane component rich in gas and dust. Gas and dust concentrate in the spiral arms which in late galaxies is 
rich in hot giants which cause the gas surrounding them to become luminous and visible. 


In addition, in the nearer spirals of later type, supergiants and their concentrations appear in the arms only 
at some large distance from the nucleus. On the other hand, the dark threads of dust stretch up to the nucleus 
and often even enter it. For example, according to Baade, in M31 supergiants begin at a distance of 2 kparsecs 

from the center. According to our estimate, luminous spiral arms begin at a distance of about 550 parsecs from 
the center, while the sharp ends of the dark fibers and separate dark nebulae approach the center to a distance 
of 80 parsecs, and many of them are within the nucleus. j 


Thus concentrations of dust occur, in fact, not only in those regions where there are many giants, but also 
between stars of smaller luminosity and unknown nature which form the inner and the most dense parts of spiral 
arms. These concentrations exist also in peripheral areas of nuclei which consist entirely of stars typical of spheri- 
cal components. Thus elliptical systems which are the nuclei of spiral galaxies contain dust clouds and fibres 


on their peripheries. Undoubtedly, they also contain gas since the formation of dust fibres can only be thought 
of in terms of gas currents. 


As already pointed out by us in [17], in M33 one finds dust up to a distance of 30 parsecs from the nuclear 
center. However, supergiants are not found within a distance of 100 parsecs from the center, and its visible 
nucleus, as we have shown, may not contain (or may contain only a small number) of stars of type II population 
according to Baade, For this reason, nuclei of spirals similar to M33 will not be considered here. If we included 
them among the elliptical systems consisting of type II population we would find in them a considerable gas and 
dust saturation. In all the above ares gas will not be luminous in the absence of hot stars. 


The gas-dust currents usually curve in the direction of the spiral arms, are radially placed relative to the 
nucleus and sharpen on the side of the nucleus, while on the outer side they branch out (and not just widen). 
Currents and "canals" of dust cannot occur without motion along them, and the question arises what is the 
direction of motion of matter in them, It is difficult to imagine that gas~dust currents flow into the nucleus from 
the periphery from all directions. If this were the case, and could somehow continuously occur during an inter- 


val of time comparable to the age of the galaxies, the latter would be saturated with gas and dust. This appar- 
ently is not the case. 
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Fig. 2. a) NGC 4485-90; b) NGC 7678; c) NGC 5278-79; d) NGC 205; e) Mol and 


NGC 5195. 


It is more probable that gas and dust flow out of the nucleus and down the magnetic lines of force in the 


plane of the system. In that case, however, the nucleus should have a store of diffuse matter which should be 


reproduced continually within it. 
hat there is a spiral of an early type, namely NGC 4826, with 


In this connection we wish to point out t 
nally strong band (wide and thick) is formed in this 


arms tightly joined to each other and twisted. An exceptio 
spiral and surrounds the periphery of the nucleus giving a strongly asymmetrical form to this large nucleus. 


Furthermore dark matter cannot be seen in continuous spiral arms. 
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Thus in some spirals diffuse matter is concentrated in the immediate neighborhood of the nucleus, half 
of it in its periphery, and not in the spiral arms. 


We must also consider the galaxies NGC 2537 and NGC 4038-9. The first of these has the form of a horse- 
shoe consisting of two grains touching each other, i. e., well defined stellar clouds, which in such cases always 
consist of supergiants, A band consisting of clouds of supergiants comes out of the center of the horseshoe in 
the outward direction. All this is on a weak stellar background. In this galaxy, similar in structure, and partly 
also in form, to spirals of type Sc, there is a bright periphery and the place of a nucleus is taken by weak lumin- 
escence, It is possible that the saturation of nucleus with dust is taking place here. 


Judging from the drawing of Rosse, NGC 3434 appears to be a fairly regular spiral with very bright inner 
arms having the form of a ring in the drawing. The central area within the ring, where the nucleus should be 
situated, is weaker than the outer arms. 


One may ask whether in the latter cases one is dealing with a particularly large accumulation of dust in 
nuclei of spiral galaxies, dust being always found in the vicinity of nuclei of spiral galaxies and sometimes 
approaching their center. 


According to Mayall [4], gas emission lines in SO and Sa galaxies are observed only fairly near the nucleus, 
and, evidently, this is also often the case in Sb and Sc galaxies. The percentage of galaxies which emit gas 
increases in the direction SO — Sc up to 88%. Thus gas is concentrated in, or in extreme cases luminesces, in the 
central areas of nuclei just as in elliptical galaxies (when found), while in both cases dust is observed at large 
distances from the center; in the case of spirals it is usually found on the periphery. Since there is no reason to 
assume that conditions leading to the emission of the line 4 3727A (special observations probably show other 
lines too) arise only in the central parts of these elliptical systems, we may conclude that gas and dust are differ- 
ently distributed in them: gas in the center, dust nearer the periphery. 


A series of spiral galaxies is known in which bright lines in the spectrum of the nucleus are very wide, 
showing gas motion with velocities up to 4000 km/sec. An inevitable consequence of this should be the outflow 
of gas from the center to the periphery. Comparing this with the above considerations it appears to be a new and 
strong argument for the hypothesis that gas can condense into dust, and that for some reason favorable conditions 

for this happening appear to occur in the periphery, in spite of the fact that one would think that when the central 
concentrations of gas move outwards towards the periphery their density decreases. The problem of the origin 
of the luminosity of gas in such systems consisting of type II population still remains unresolved. It is possible 
that the luminescence of gas is due to hot stars of low luminosity, Nuclei of planetary nebulae in our galaxy 
belong to this type. Furthermore, there are increasing indications that stars of type 0 of low luminosity are mem- 
bers of type II populations. In the absence of necessary data the amount of gas in the described systems cannot 
be estimated. The absence of appreciable radio emission from elliptical galaxies does not provide much infor- 
mation on the presence of gas in them since even in the irregular and spiral galaxies which are known to be rich 
in gas, radio emission has so far been discovered only in single systems. 


In conclusion one may say that in major spherical stellar systems, such as nuclei of spirals and elliptical 
galaxies, the presence of diffuse matter is not at all infrequent. The problem of its presence in globular clusters 
remains open. Indications of the presence of sometimes considerable amounts of diffuse matter in elliptical 
galaxies, and of radial motion of dust in NGC 5195 and nuclei of spirals, support the hypothesis that when there 


is a sufficient store of diffuse matter the formation of spiral arms is due to the outflow of this matter down mag- 
netic lines of force, 


5. The origin of diffuse matter in spherical stellar systems. There are three theories of the origin of 
diffuse matter,or, more precisely, ofits origin in the neighborhood of the sun, i. e., near the plane of the galaxy. 
The earliest theory was that the observed diffuse matter is the remainder of the material from which stars were 
formed. We showed in 1931 that the loss of gas by stars is sufficient for the formation of stars from diffuse 


matter thus produced, Recently, some authors have openly, or otherwise, suggested that both stars and the diffuse 
matter are produced together from some other apparently superdense substance. 


However, what can be said about diffuse matter in spherical stellar systems? 
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By the presence of diffuse matter in a stellar system we mean the case where the density of diffuse matter 
in the system is appreciably higher than in the surrounding metagalactic space. Every year there are more indi- 
cations that the density of metagalactic diffuse matter is very nonuniform, However, the passage of galaxies 
which do not contain diffuse matter (in the above sense) through clouds of metagalactic matter, and the passage 
of globular clusters through the galactic plane, cannot lead to an enrichment of such spherical systems with 
diffuse matter since these passages take place at large relative velocities. If such systems contain diffuse matter 
having a density less (say, by a factor of 1-2) than that in the galactic plane, the take-up of diffuse matter at 
large relative velocities will evidently not take place either. On the contrary, the denser medium will on 


collision remove gas from the system which itself will pass freely through the medium and will emerge from 
it free of gas. 


However, one cannot exclude the possibility that in the presence of orbital motion of diffuse masses and 
spherical systems their encounters will occur at small relative velocities, and with the presence of matter in 
the system, even with a small excess density, the capture of an additional amount of diffuse matter will be 


possible. Such cases are probably rare but as a result some spherical systems may occur which are rich in foreign 
diffuse matter. 


A certain ammount of gas in spherical systems may be produced from internal resources. Nonstationary 
objects are present in small numbers in globular clusters. Thus a new star of T Scorpius was observed in M80 
in 1860. A planetary nebula enters into M5. Ordinary stars will give certain amounts of gas. There is no reason 
to suppose that the loss of gas by stellar atmospheres is smaller than in the stars of the main sequence. The latter 
is, however, unknown. According to the hypothesis of V. G. Fesenkov on the evoluatio of stars, this loss is con- 
siderable. 


In natural conditions, gas will concentrate towards the center of spherical systems, as it would concentrate 
in the planes of spirals because of their higher density. If one assumes that gas has a density of 10°” g/cc within 
the central part of a globular cluster having a diameter of 20 parsecs, where stars are still sufficiently closely 
packed, one obtains a figure of 10°° g for the mass of the gas. If the number of stars in a cluster is estimated to 
be 10° then the mean loss per star is 10 g. If in such clusters stars lose their mass in 1014 years (which is smaller 

_ by one order than the usual assumed age of these clusters), then the mean mass-loss via corpuscular emission is 
10 g per year for each star. Figures between 10 and 10” g per year have been found by various authors for 
the corpuscular emission from the sun. For this reason corpuscular emission of even stationary stars of average 
luminosity may, over the assumed length of time, produce with the globular clusters a diffuse medium with a 
density comparable with that in the neighborhood of the sun. 


Similar results are obtained for elliptical galaxies. 


In spite of all this, indications of the outflow of diffuse matter from the center to the periphery suggest 
that this matter is present in the central parts of spherical systems from the very start and is not accumulated 


in them subsequently. 
Shternberg State Astronomical Institute Received July 26, 1956 
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THE DETERMINATION OF THE ABSOLUTE MAGNITUDES OF CEPHEIDS 
AND SUPERGIANTS AND THE ZERO POINT OF THE PERIOD- 


LUMINOSITY RELATION 


Iu. P. Pskovskii 


The radial velocities of 182 stars confirm the calibration of supergiants found 
by Morgan and Keenan from the effect of galactic rotation. The distance modulus 
for « Per, derived by different methods is discussed and the most probable value is 
found, My = — 47.0 + 01.3 (p.e). Using o Per as a mark, the spectral types, 
absolute magnitudes and luminosity classes of 6 supergiants and 6 subgiants of spec- 
tral types F5-G5 and also 6 Cep and 7 Aql were determined from equivalent line 
widths. A correction to Kukarkin’s zero point of the period — luminosity relation 
was found, — 0™.76 + 0" 40, from derived median photographic absolute magni- 
tudes. This value is in good agreement with those derived by other investigators 
using different methods, with corrections, wherever needed, to a unified mean 
system of stellar astronomical parameters. From all the determinations — 0.83 + 
+ 0104 was found as the most probable correction to Kukarkin's zero point. 


. Calibration of the Absolute Magnitudes of Supergiants from the Effect of 


1 
Galactic Rotation on Radial Velocities 


The mean absolute magnitudes of supergiants are determined from the effect of galactic rotation by esti- 
mating their mean distance according to Oort's formula 


V,=K-+ Arsin 2 (l — ly) cos? b — Vy cosi, (1) 


with the usual notation. 

This method was applied to supergiants of spectral types B— K by J. Greenstein [1] andR. E. Wilson [2] 
in 1940; Wilson considered 205 stars. In 1951 Morgan and Keenan [3] used the method to determine the mean 
absolute magnitudes of supergiants of luminosity class Ia (8 stars) and class Ib (8 stars). 

In order to determine the parameters of Oort's formula we used the radial velocities of 308 c-stars of spec- 


tral types cB6-cM5 which have been classified by different authors as stars of luminosity classes Ia, Ib and II or 
as c-stars. The data on radial velocities were taken from Wilson's catalogue [4] and from the list of Feast, 


Thackeray and Wesselink [5]. 
We excluded stars with large radial velocities and proper motions and stars in high galactic latitudes. 
The stars were divided into groups according to their positions in the sky and the radial velocities and 


galactic coordinates of each group were averaged. Thirty-one equations of the type of (1) were obtained, the 
solution os which was obtained by least squares, At the same time the probable errors of the parameters of (1) 


were obtained. 


The results are given in Table 1. 


The table shows that our value for Vo agrees with that obtained by Parenago [6]. Our absolute value for 
K is greater than that of Wilson because the latter did not drop early B stars with a positive K effect. Our value 
for the longitude of the galactic center (19) is close to the value obtained from infra-red observations (330°). 


TABLE 1 
coil ea 
No. 0 
Author Spectral | super |Voxm/sec | K,km/sec . i I; 
EPS giants es 

R. E. Wilson, 1940 cB—cK 205 49.7 —2.0+0.8] 16.3 | 324°4 
Our work, 1956 cBG6—cM5 | 308 | 14.2+41.4 |—3.4 $1.0] 41.7 |] 334.2 
P, P. Parengo, 1951 cO — cM 174 19.7 — — = 


Wilson [7] and Mel'nikov [8] have shown for a number of objects and Ogorodnikov [9] has shown theoreti- 
cally that K depends linearly on distance: 


K = Kg + Kj. (2) 


We have obtained Kp and kK, for 182 stars of luminosity classes Ia, Ib and Il in nine distance groups, with 
distances according to the calibration by Morgan and Keenan [3]. K and Ar for each group were obtained by 
least squares. Taking A = 19.5 km/sec per kiloparsec and solving a set of 9 equations such as (2) by least squares, 
we found that Kg is practically zero (Ko = 0.5 + 2.7 km/sec) and Ky =— 6.5 + 0.3 (probable error). This value 
of Ky is close to that obtained by Mel'nikov [8] for long-period Cepheids if the difference in the assumed value 
of A is taken into consideration. 


After eliminating the K term from (1) we determined the mean weighted distances of 17 star groups dis- 
tinguished by similarity of spectral type and luminosity class. The following formula was used for the mean 
weighted distance: 


>} (Ki + sin 2 (1 —1,) cos? b) (V,, — Vg cos 2) 
A}! (Ky + sin 2 (1 — ly) cos? b)? 


In order to determine the mean absolute magnitudes of the supergiants we used Parenago's method [10] to study 
light absorption (Greenstein and Wilson used cruder methods) and investigated the dispersion of the distance 
modulus for stars in each group. The results are given in Table 2. 


The probable errors of the mean distances were calculated in the usual way and the probable errors of the 


mean absolute magnitudes were computed from Equation (30) of Parenago’s article on interstellar absorption [10] 
with certain simplifications. 


The last column of the table gives the mean absolute magnitude of the stars in a group when the absolute 
magnitudes are in accordance with Morgan and Keenan's calibration [3]. The last three lines give the mean 


absolute magnitudes of Ia, Ib and II luminosity classes with probable error from the meanconsistency of values 
of the mean absolute magnitudes for groups of the table. 


This is a confirmation of Morgan and Keenan's calibration of supergiants of luminosity classes Ia and Ib. 
Stars of class II have a large dispersion of radial velocities and the calibration method cannot be applied to them. 


TABIE 2 


Spectral |Meanspec- |No.of| _ | As a 
type ; m Yr M My 
eS tral hee Neos te kps is from [3] 
B7—AO BMa IO |) PO M27 Od8. 1b 7b ONG —7m() 
DOSING A3Bta 7 5.5 1540629 26.65. O18 ee pe) 
rO—GO Fla 4) 5.0 1.20 420.16 TO SENG aeRO) 
G8—M2 MOlah 8 6.4 (). 76 45 0.36 Sh FEI) —6.3 
B8—B9 BOTb 9 Orel 0.66 + 0.16 Hh) FENG hy 
AN—F0 AAIb 9 ee ().68 + 0.24 alti Of EMIS sale '5) 
F2—F8 5b 8 yee ().36 4-20.45 —4.6 + 0.9 (hat) 
GO GOIb 9) 5.6 0.40 +4 0.44 =A) 0.) ee 
eo, G2Ib 8 nee Leal 25 Oyars sxeify PSX Set 
G3—K0 GG6Ib 13 5.3 0.49 4.0.20 ps eral) —h 5 
K1—M2 K4Ib 12 4.6 OAD AE Oe ath oh NY Sas 
A5—F6 r2U1 1A DES Oise One sah) Se SO) = Dill) 
F8—GA4 GOI 16 Re 0.39 + 0.14 = OD SER ==) 
C567 GOT 13 6.2 0.38 40.13 aad 7 SOG KO) 
G8 G8Il (1 5.9 0.35 2.0.46 SD ee OnG eA 
C9=K2 KAI 1A 5.9 O29) 2.0.30 SP TE Nas zi) 
K3—M5 KG6II 12 ee 0.28 4.0.24 SY =, 
B7—M2 Ta 32 = Ee ee One ST) 
B8—M2, Ib 70 a ay Be AeA Ss 
A5—M5 Il 80 a= =e ET (ns i 


No conclusion can be drawn from these results concerning the variation, from one spectral type to another, 
of the mean absolute magnitude within a luminosity class. 


2. Absolute Magnitudes of Five Pseudo-Cepheids and Two Cepheids from Luminosity 


The luminosity relation which Hertzsprung discovered in 1905 was thoroughly investigated in 1914 by 
Adams and Kohlschutter and was used to determine the absolute magnitudes of stars. However, the numerous 
determinations of absolute spectral magnitudes at the Mount Wilson, Victoria and other observatories from 1923 
to 1935 are now obsolete because they did not take interstellar absorption into account. 


An important advance was made by Morgan and Keenan's [11] two-dimensional spectral classification 
and calibration of supergiants with respect to absolute magnitude. At present, however, we have no determina~ 
tion of individual absolute magnitudes of supergiants from luminosity. 


The basic difficulty in determining individual magnitudes of supergiants is the almost complete absence 
of reference stars, that is, supergiants whose absolute magnitudes are known exactly. 


Among the pseudo-Cepheids only  Persei belongs to a moving cluster in Perseus and has a cluster parallax 
(Table 3). However, the determination of the parallax of  Persei from the spectrum-apparent magnitude and 
color-apparent magnitude diagrams diverges from the group parallax. The discrepancy originates in the in- 
accuracy of the old calibration of the main sequence [3]. Using the 1956 calibration of the main sequence by 
Johnson and Hiltner [22], we corrected the distance moduli for @ Persei and obtained an average value which 


is in good agreement with the average obtained from group parallaxes. 


The apparent magnitude of a Persei corrected for absorption is 117.4 according to Eggen, Harris and others, 
so that the absolute magnitude of & Persei is taken as 


Vie 0 ORD 


The agreement of the absolute magnitude as determined froin the cluster and the calibration of the 
spectral types F5Ib serves as a proof that @ Persei actually belongs to the cluster. 


TABLE 38 


a rrr 


Method 


Group method of [12] 
Group method of [13] 
Group method of [14] 
Group method of [15] 


Average for o Persei from group parallax: pl 4A OF 


Spectrum-magnitude 

diagram [16] 0"006 ee 5M 4 
Spectrum~magnitude 

diagram [17] 0.0055 6.3 5.0 
Spectrum-magnitude 

diagram [18] 0,004 6.85 5.85 
Color-magnitude 

diagram [19] 0,007 5.85 5.7 
Color-magnitude 

diagram [20] 0.010 50 5.0 
v-component of 

proper motion [17] 0.008 5.4 5.4 


Average from distance modulus of cluster 5™,.5 4 0.1, 


The trigonometrical parallax of a Persei in Jenkins’ catalogue [23] is considerable, being + 0". 029 + 
+ 0%005. But this is an error: the parallax in the catalogue of the Allegheny Observatory [24] is + 0"007 + 
+ 0" 005 and not 0."038 + 0." 006 as Jenkins assumed in connection with the absolute parallax. Our corrected 
trigonometrical parallax for « Persei is + 0." 006 + 0." 004 (mean square error) in Jenkins’ catalogue system. 
Accordingly Myrig = — baal we a Pals 

For the purpose of determining the absolute magnitudes of five pseudo-Cepheids and two Cepheids we 
obtained 40 spectrograms of six supergiants, six subgiants and two Cepheids using the 120-cm reflector of the 
Crimean Astrophysical Observatory of the USSR Academy of Sciences. The spectral dispersion was 35.6 A/mm 
for Hy. Ilford Zenith plates were used. An MF-4 microphotometer was used to process this material. The 
continuous spectruin in the spectrograms was taken from the traces for & Persei in the atlas of Hiltner and 
Williams [25], for the sun in Minnaert's atlas and for y Cyg in one of Iwanowska’s articles [26]. 


Table 4 gives the equivalent breadths of 18 lines and of the G band for stars of the program. The asterisks 
denote lines used in determining the absolute magnitudes and spectral types. 


The stellar spectral types and absolute magnitudes were determined from diagrams whose ordinates were 


logarithms of the ratio of equivalent breadths of pairs of lines and whose abscissas were the absolute magnitudes 
(visual). 


In determining the spectral types it was assumed that the absolute magnitudes of the supergiants are 
— 4.5 (from Morgan and Keenan's calibration), and the absolute magnitudes of the subgiants were calculated 
from the trigonometrical parallaxes in Jenkins' catalog. The criteria of spectral types were the ratios of equi- 


valent breadths of the hydrogen lines Hs and Hy to the lines of ionized metals, as well as the ratio of the 
equivalent breadths of the lines 4 4077 Sr* and A 4383 Fe. 
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The reference stars were 9 Cyg and a Per (with 
the spectral type taken as exactly F5) and p Her and 


log es 9 Peg (with the spectral type taken as exactly G5). 
; Thirteen pairs of lines were selected for which the 
100 6 Cyg change of the ratio of equivalent breadths of a pair of 
ue be FS ot Per lines in going from a Abe of one ea type to a star 
of another spectral type is proportional to the change of 
260 3 spectral type from one star to the other. This require- 
£8 Pos ment corresponds to parallelism of the lines which connect 
070 Se points that correspond to ratios of equivalent line breadths 
Sac ee in spectra of the same type and to equality of the dis- 


050 age Es tances between parallel lines through points represent- 
SS Gg seas ing values of the ratios for the types F5, F6,. ..G5 
G50 ts 


i eee (Fig. 1). By means of this grid constructed according to 

a ry Le ar equivalent line breadths of the reference stars we deter- 
gE ee ee mined the spectral types of the stars in Table 5. For 

MIO rie determination of the absolute magnitudes the reference 


© 9Peg supergiant was @ Per (My = — 4.™0 + 0,"3). The grid 

20 ig ee was analogous to that constructed for spectral types. 
See ee The differences between the logarithms of the ratios of 

aro ae equivalent breadths for line pairs in the spectra of @ 
| ‘Pi ae Cyg and Her provided a basis for ten intervals of the 
Tia are 2" -9" <8", grid; the line connecting points of the drawing corres- 


ponding to the logarithm of the ratio of equivalent breadths 
of line pairs for stars of a single spectral type (F5): 0 

Cyg and @ Per provided the direction parallel to which 
spectral types from the logarithm of the ratio the grid was plotted. Fifteen pairs of lines were used: 

of equivalent breadths of A 4102 (H) and the ratios of the equivalent breadths of A 4077 Srt and 


X 4045 (Fe) in the spectra of these stars. The A 4444 Ti* to the equivalent breadths for unionized 
symbol © denotes log (Wajo2 /Waq4s) in the 


spectra of the-reference stars and the dots 
denote the same in the spectra of other stars. Table 5 gives our values of the absolute magnitudes 


with the probable errors from internal consistency. The 
table shows that @ Per, 8 Aqr and 9 Peg must be assigned 
to luminosity class Ib-II. Comparison with the spectral types obtained by Gathier[28] confirms out conclusion. 


Fig. 1. Grid. for the determination of stellar 


elements. 


The spectral types of Cepheids were also determined by a different method, which is the relation between 
the equivalent line breadth and the spectral type of supergiants. This relation was plotted from the data in 
Tables 4 and 5 for Hs Hy, A 4383 Fe, A 4227 Ca, A 4325 Fe and the G-band. 


From the first three lines the spectral types were earlier than by the other criteria (the Adams-Russell 
phenomenon) but the spectral type obtained from all criteria is close to the average by the first method. 


The median absolute magnitudes of 5 Cepand 7 Aql were obtained from the luminosity curve and Eggen's 
color index [29]. The average of four values ior the median visual absolute magnitude of. o& Per is —3.™41 3 
+ 0.™04 (probable error from internal consistency of the four values) and the median visual absolute magnitude 


%, n Aql (from a single phase) is — 4.1735 4 0.740, 


The median photographic absolute magnitudes are as follows: 


8 Cep Mpg = —2770 +0704 


qAd Mpeg = 3™67 + 0740 
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Zero-Point of the Period-Luminosity Relation of Long-Period Cepheids 


The values obtained for the median absolute magnitudes of Cepheids enable us to determine the zero- 
point of the period-luminosity relation, which determines the scale of distances in the universe; this accounts 
for the great interest in the zero point and the great amount of work on its determination. 


In 1952-1954 it was shown that between the period-luminosity relations of short-period and long-period 
Cepheids there is a discontinuity at 1.%5-2™M,_ This observed fact has been interpreted in different ways. 


Baade [30] states that the mean absolute magnitude of short-period Cepheids is 0.!0 and that the correc- 
tion to the zero point of Shapley's period-Juminosity relation [34] is —1.™5 for long-period Cepheids. 


Parenago [31] states that the mean absolute magnitude of short-period Cepheids is + 0.5. This order of 
magnitude is obtained both from the proper motions and from the parameters of a subsystem of short~period 
Cepheids [28]. The correction to Shapley's zero-point is therefore the order of — 1™, 


Kholopov [32] states that Shapley's zero point is correct and that the mean absolute magnitude of short- 
period Cepheids is of the order of + 2™, 


Kukarkin [33] showed in 1949 that the period-luminosity relation can be represented by three formulas, 
one of which, for Cepheids with periods from 1 to 9 days, is the following: 


Mg med= —0"74 —1™67logP (3) 


Kukarkin obtained the correction — 0.!"5 to Shapley's zero point for long-period Cepheids. 


Using Formula (3) we obtained from the photographic median absolute magnitudes of 6 Cep and 7 Aql 
the correction — 0.™76 + 0.40 to Kukarkin's zero point (the error includes the entire dispersion of absolute 
Cepheid magnitudes, which is + 0.30). 


There are many recent determinations of the correction to the zero point, some of which have been con- 
sidered by Parenago [31]. He determined the correction to the zero point from the mass-luminosity relation for 
supergiants and from the trigonometric parallaxes of Cepheids. He converted to his system the results obtained 
by Blaauw and Morgan for the zero point from the proper motions of 18 Cepheids, and by Kukarkin from the effect 
of galactic rotation on the radial velocities of Cepheids, and also used the investigation by Thackeray and 
Wesselink of short-period Cepheids in the Magellanic clouds. 


TABLE 6 


Method Correction to Kukarkin's Correction to Shapley's 
zero point zero point 


Distance to center of Galaxy from radial velocities 
of Cepheids [39] 

Proper motions of Cepheids [40] 

Proper motions of Cepheids [41] 

Short-period Cepheids in Magellanic Clouds [42 ] 

Bright stars in globular clusters of Magellanic 
Clouds [43] 

Integr. mag. of globular clusters of Magellanic 
Clouds [43] 

Mass-luminosity relation for Magellanic Clouds [31] 

Trig. parallaxes of Cepheids [31] 

Abs, mag. of Cepheids from luminosity effect 

Rel. bolom. mag. — radius and temp. effect 


This summary of determinations of the zero point can be supplemented BY a few addimiogal ee om 
verted to the system employed by Parenago (mean photographic absolute magnitude of shore“ verge Cep r s 
+ 0.5, light absorption by Parenago's method, velocity of the sun with respect to the Cepheids 18.6 km/sec 
and distance to the center of the Galaxy from radio observations 8.2 kiloparsecs [35] ). 


Besides Blaauw and Morgan the zero point was determined from the proper motions of 18 Cepheids Bye 
A. Ja. Filin in 1953. But in calculating light absorption he exaggerated the dispersion of the absolute bedi ogee: 
of his Cepheids (+ 0.63 instead of + 0.44 when calculated as indicated by Kukarkin [33] Ghe aa all 
of this error and conversion to Parenago's system (according to Parenago [21] the reduction for the catalogues is 
+ 0.195) gives a correction to the zero point of — 0.™M9 4 0,5 (mean error), which agrees within ye limits of 
error with Blaauw and Morgan in the unit system (Table 6). Filin's correction to Kukarkin's zero point was 


originally — 1.™gi + 0,144 (rms error). 


In 1954 H. Weaver determined the correction to Shapley's zero point from the discrepancy between the 
estimates of the sun's distance from the center of the Galaxy as determined from the radial velocities of globular 
clusters and Cepheids. In its original form the correction was ~ 1.™56 + 0.™30 (mean error). When converted 
to Parenago's system the correction is of the same magnitude as by the preceeding methods. 


A correction to the zero point was also obtained from the distance modulus to the Magellanic Clouds. 
Considering the new work of Thackeray and Wesselink on the apparent magnitude of short-period Cepheids in 
the Magellanic Clouds, we obtained for the distance modulus to the large Magellanic Cloud 18.11 + 0.™4 and 
to the small Magellanic Cloud 18.12 4 0.™4 (probable error), considering the mean photographic absolute 
magnitude of short-period Cepheids to be + 0.15. Two additional methods for determining the distance modulus 
have been suggested by Shapley. By the first of these methods (comparison of the apparent magnitudes of bright 
stars in globular clusters of the Galaxy and Magellanic clouds) the distance modulus to the large Magellanic 
cloud is 18.™5 + 0.4 and to the small Magellanic cloud 18.™6 + 0.4 (in our system). By the second method 
(comparison of the integral stellar magnitudes of globular clusters in the Galaxy and Magellanic Clouds) the 
distance moduli are 18.™2 + 0,™5 (large Magellanic cloud) and 18.™2 4 0.™5 (small Magellanic cloud). The 
results are given with mean errors. The average distance moduli are identical for the large and small Magellanic 
clouds and for each method can be taken as the average for the large and small clouds. 


Using the period-luminosity relation in apparent magnitudes for the Magellanic clouds as given by 
Kukarkin [33], we determined the zero point by each of the three methods. Opolski and Krawiecka [44] by deter- 
mining the radius and estimating the effective temperatures of 5 Cep and » Aql obtained their mean bolometric 
magnitudes and a correction of — 1.4 to Shapley's zero point. Conversion to the more accurate temperature 
scale of Parenago [31] reduces the correction to 0.12, 


Table 6 is a list of ten determinations of the correction to the period-luminosity zero points of Kukarkin 
and of Shapley (1930), which were obtained from 1953 to 1956. This list does not include the estimates of the 
zero point from novae (Kopylov [36], Buscombe and Vaucouleurs [37]), uncertain consequences of the large 
real dispersion of maximum absolute magnitudes of novae and the estimate of the zero point from parameters 


of the subsystems of long-period Cepheids (Kopylov and Kumaigorodskaia [38]) which is apparently not free of 
systematic errors. 


There is remarkably good agreement among the majority of corrections to Kukarkin's zero point (assuming 
the mean photographic absolute magnitude of short-period Cepheids to be + 0.5). The same agreement is not 
sound among corrections to Shapley’s zero point (mean photographic absolute magnitude of short-period Cepheids 
0.°'0); the corrections to the zero point according to the short=period Cepheids in the Magellanic clouds are 
greater by 0.76 in absolute value than by other methods (Table 6). Thus the most acceptable mean photographic 
absolute magnitude of short-period Cepheids is that of Parenago, and Parenago's view regarding the discontinuity 
between the mean absolute magnitudes of short~period and long-period Cepheids is closest to the truth. 


The most probable correction to Kukarkin's zero point for the period=luminosity relation is — 0.Mgg + 0.™M04 
(rms error), which corresponds to multiplying Shapley's distance scale by the factor 1.8. 


P. K. Shternberg State Astronomical Institute Received October 1, 1956 


28 


LITERATURE CITED 


{1} J. L. Greenstein, Proc. Natl. Acad. Sci. 26, 259 (1940). 

[2] R. E. Wilson, Astrophys. J. 93, 212 (1941), 

[3] P. C. Keenan and W. W. Morgan, Astrophysics (ed. J. A. Hynek),Chap. I (New York, 1951). 

[4] R. E. Wilson, General Catalogue of Stars’ Radial Velocities (Wash. Obs. Publ., 1953). 

[5] M. W. Feast, A D. Thackeray and A. J. Wesselink, Mem. R. A. S., 67, part 2 (1955). 

{6] P. P. Parenago, Trudy GAISh 20, 26 (1951). 

[7] 

[8] 

[9] 
[10] 


E. Wilson, Astrophys. J. 92, 170 (1940), 
. A. Mel'nikov, Trudy GAO 64 (1950), 


P. Parenago, Astron. J. (USSR) 22, 129 (1945). 
[11] W. W, Morgan, P. C. Keenan and E. Kellman, An Atlas of Stellar Spectra (Yerkes Obs., Chicago, 1943). 
[12] 


[13] 


P 
W 
K. Bottlinger. Astron. Nachr. 212, 142 (1921), 
. Rassmusson, Lund. Medd., Ser. II, 26, 28 (1921). 
H 


O 
K. F. Ogorodnikov, Astron. J. (USSR) 21, 43, 46 (1944). 
N 
[14] H. Schloss,Z. Astrophys. 11, 117 (1935). 

[15] W. Smart, and A. Ali, Monthly Notes 100, 560 (1940), 


[16] R. Trumpler, Lick Obs. Bull. 14, 154 (1930). 


[17] N. G. Roman and W. W. Morgan, Astrophys. J. 111, 426 (1950). 
[18] 


R 

N 

M. K. V. Bappu, Monthly Notices 114, 687 (1954). 
[19] O. J. Eggen, Astron. J. 60, 65 (1955). 

D 


[20] D.L. Harris, Astrophys. J. 123, 371 (1956). 
[21] W. W. Morgan and J. L. Johnson, Astrophys. J. 117, 313 (1953). 
[22] H. L.Johnson and W. A. Hiltner, Astrophys. A. 123, 276 (1956). 


[23] L. F. Jenkins, General Catalogue of Trigonometric Stellar Parallaxes (Yale Univ. Obs., New Haven, 
Conn., 1952). 


[24] Publications Allegheny Obs., 10 (1956). 


[25] W. A. Hiltner and R. C. Williams, Photometric Atlas of Stellar Spectra (Univ. of Michig. Press, Ann 
Arbor, 1946). 


[26] W. J. Iwanowska, Wilno Obs. Bull., 18 (1937). 

[27] N. G. Roman, Astrophys. J. 116, 122 (1952). 

[28] P. J. Gathier, Rech. Astr. Obs. Utrecht., 13, part 2 (1955). 
[29] O. J. Eggen, Astrophys. J. 113, 367; 114, 144 (1951). 

[30] W. Baade, Trans. J. A. U. 8, 397 (1954). 

[31] P. P. Parenago, Variable Stars (USSR) 10, 193 (1954). 
[32] P. N. Kholopov, Astron. Tsirkuliar 148 (1954). 

[33] B. V. Kukarkin, Parerner: Zvezd. 7, 57, 69 (1949). 


99 


30 


H. Shapley, Star Clusters, Harv. Mon., 2 (1930). 

P. P. Parenago, Astron. J. (USSR) 32, 226 (1955). 

I. M. Kopylov, Izv. Krim. Astr. Obs. 13, 76 (1955). 

W. Buscombe and G. Vaucouleurs, Astron. J. 60, 155 (1955). 

I. M. Kopylov and R. N. Kumaigorodskaia, Izv. Krim. Astr. Obs. 15, 169 (1955). 
Weaver, Astron. J. 59, 375 (1954). 


Blaauw and H. Morgan, Berlin. Astron. Nachr. 12, 95 (1954). 


Thackeray and A. Wesselink, Monthly Notices A. S. S, A. 13, 99 (1954). 


H. 

A. 

A la. Filin, Biul. Stalinab. Obs. 10, 14 (1954). 

A 

H. Shapley and V. McKibben, Proc. Nat. Acad. Sci. 40, 1 (1954). 
A. 


Opolski and J. Krawiecka, Contr. Wroclaw Obs. No. 11 (1956). 


THE DISTRIBUTION OF ABSORBING MATTER AND THE SPIRAL 
STRUC TURE OF THE GALAXY 


N. E. Kurochkin 


The absorption of light in a heterogeneous medium is given in stellar magni- 
T2 ‘ 
tudes by A(r)=m—my=1.086x \ e(r)dr, where x is the absorption coefficient 
; 
per unit density and P(r) is the density, For small values of rg 1 the density 
P(r) may be taken as constant. The distribution of density in the galaxy may then 
be studied, without any additional assumptions, to within x. Absorption studies 
have been carried out in this way by Kharadze [1] and van Rhijn[2]. The density 
distribution as a function of the z coordinate and the radius (the latter was obtained 
from the absorption data of van Rhijn for six Kapteyn areas with b =<+ 3° assuming 
spherical symmetry) may be represented by 


D,= 24.38. 10-4¢e—2/216-5. Dry = 19.28.40-4 e—(R—7200)/1010 


where r and z are expressed in parsecs. 


The deviation from the logarithmic law of distribution of density for low 
latitude directions agrees with the spiral structure observed in the solar vicinity, 
determined from the distribution of maxima of B-stars (according to Bergedorf 
Spektral Durchmusterung [11] ) and the maxima in the distribution of neutral 
hydrogen (radio observations [77] ). The spiral structure in the solar vicinity is 

very irregular but its general form can be deduced from the sum total of obser- 
vations. More attention should be given to theoretical and experimental studies 
of physical properties of interstellar medium in order to determine x and the 
zero point of the derived regularities. 


Studies of galactic structure meet with very great difficulties because of the presence in the galaxy of 
lightabsorbing medium. It is difficult to say whether it is the absorbing matter or the stars which are the more 
irregularly distributed over small regions of space. Undoubtedly, the visible structure of the stellar system, 
which is mainly determined by weak stars, is markedly distorted by absorbing matter. A detailed study of the 
general laws of distribution of absorbing matter should therefore precede studies of the distribution of stars. 
Systematic studies of absorption of light which were begun in the thirties have produced a number of papers on 
the subject. Kharadze [1] has given a sufficiently detailed review of this work. Among these papers only those 
are fundamental which give a sufficiently accurate and extensive series of observations of the color indices 
from which, knowing the spectral class of the stars, one can obtain the color excesses and full absorption in the 
photographic and visual regions. Undoubtedly, an accurate determination of the system of color indices is o 
major importance, and in considering the data of different authors differences in the systems must be taken into 
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account, We shall show that the same treatment applied to results of different authors leads to strongly differ- 
ing conclusions, This is especially the case in determinations of light absorption for weak objects, which corres- 
pond to those parts of the galaxy which are remote from the sun, The determination of color indices and spectra 
of weak stars is of importance if we wish to study general regularities in the distribution of dark matter. The 
distribution of absorbing matter fluctuates so much that absorption studies over small parts of the sky give no 
clue as to the general distribution laws. At the same time, and for the same reason, studies of weak absorbing 
objects over the whole of the sky is also pointless since it would only strengthen our ideas on the chaotic way in 
which dark matter is distributed. Detailed studies should only be carried out along those directions in which 

one would expect to discover general and most characteristic properties of the distribution of matter, and, parti- 
cularly, where one can penetrate to large distances from the sun where near clouds of dark matter are absent. 
The existing data on the absorption of light and the galactic structure in the region 1-2 kiloparsecs from the 

sun may be used in a preliminary study of the distribution of absorbing matter, and this is done in the present 
paper. Accurate photoelectric data on color indices of weak stars in chosen areas, where one would expect struc~ 
tural properties of the galaxy to show up, and on visibility at large distances (the work on this is not very volum- 
inous and requires powerful instrumentation) could be used to solve the problem of the distribution of absorbing 
matter sufficiently accurately and completely. The choice of such directions is bound up with the elucidation 
of the position of spiral arms in the solar region, and preliminary studies of the distribution of dust relative to 
spirals, stellar clouds and other stellar formations. 


Method. The absorption of matter in a heterogeneous medium is given by 


dE = —k(r) Edr, © (1) 


where E is the intensity of light, k(r) the absorption coefficient depending on the distance r and, in general on 
the wavelength A, although usually one finds the total photographic absorption on the international scale of 
Stellar magnitudes and color indices which gives k(r) averaged over the given spectral region. Let us express 
k(r) in the form k(r) = kp (r), where x is the absorption coefficient per unit density and depends only on the 
physical properties of the medium. If we study the distribution of matter over a large region of space then 

this assumes the independence of k of 1, which is true in the galaxy only to a first approximation. Differences 
are possible in the composition of absorbing matter at different distances from the galactic nucleus, at different 
distances from the plane of the galaxy and even in different dust clouds. As the various methods of studying 
interstellar matter develop it will be possible to bring out these differences. 


Integrating the equation dE =x Pp (r)Edr over a range of r we have 


ties eee \ o(r)dr. (2) 


Hence, multiplying by — 2.5, and converting to logarithms to the base 10 we have 


A(r)=m—m, = 1.086. | o(r)dr. @) 


TY, 


If we have values of A(r) for different distances along the line of sight, the differences between them will 
be equal to absorption over the path difference rz—1. For small differences rg 1y the density may be taken as 
constant and equal to some mean value / (r) corresponding to the difference rz—14. In that case 


~ yy — “A (re) — A(r4) 
“0 ) ae a 
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We thus have the linear or volume density of the absorbing material, depending on the method of calcula- 


tion of the absorption coefficient k which can, for example, be expressed in stellar magnitudes per parsec or 
per unit density. 


We thus have the possibility of measuring the density of absorbing matter to within the constant multiplier 
(or on the logarithmic scale, to within the additive constant), By determining x by some method we can get 
an idea of the zero point of P(r). Random density fluctuations noticeable even over large distances rg— 14, of 
the order of 100 parsecs or more, allow us to average the values of p (r) over distances of 200-300 parsecs, which 
is not reflected on the general regularities in the distribution where distances of a few hundreds of parsecs or 
kiloparsecs are significant. Naturally, the smaller the interval the more accurate are the observations, and the 
finer the structure shown, but it is debatable whether one can expect more than a confirmation of randomness in 
the density distribution froin such observations, Only a comparison of the density of absorbing matter with the 
distribution of other matter, and studies of general regularities in the distribution of matter in the galaxy, can 
have any meaning. The above method frees us, to a considerable extent, from various preliminary assumptions 
as to the nature of the distribution of matter in space, and allows us to study densities along different directions 
in the galaxy. The determination of x thus becomes an important problem. Here it is possible to carry out, in 
addition to a theoretical treatment, an experimental determination of the absorption coefficient of meteoritic 
matter in the form of dust consisting of particles 107* cm in diameter and suspended in some transparent liquid 
or gas. A comparison of such experimental data with theoretical deductions, and with observations, would enable 
us to establish the physical properties of absorbing matter quite reliably. 


Distribution of absorbing matter along the z-coordinate. To determine the density of absorbing matter 
by the method described above we used the data of Kharadze [1] and van Rhijn[2]. In the latter work the data 


of Kharadze were used to start with, but were corrected using photoelectric and photometric observations of 
various authors [3-6]. The work of van Rhijn must therefore be given precedence. Observations of Kharadze 

and other authors contain determinations of color excesses for stars of early spectral classes in chosen Kapteyn 
areas betweenll= 0° andl = 143°, Using the conversion factor y = 5.2 (from Kharadze's data; this value was 
confirmed by van Rhijn for the international system of color indices in a summary of Leyden observations) selec- 
tive absorption may be converted into full photographic absorption A(r). Both Kharadze and van Rhijn give tables 
of A(r) for different r. These tables were used in our calculations using Formula (4). 


All SA areas were split into four groups with respect to the galactic latitude: I(/ = 0-59°), I1(2 = 59- 
81°), Il(2 = 81-118°), IV (2 =118-143°). This division was carried out because of a correspondence with 
the possible distribution of spiral arms in the solar neighborhood to be discussed below. We must, however, admit 
that no real dependence on galactic longitude was found within the limits of the above groupings, nor from one 
group to another, since random fluctuations along different directions in the solar region are apparently much 
larger. All the same, some general regularities can be observed in regions which are more distant from the sun. 
Table 1 shows the grouping of the data. The magnitudes in high latitude areas 14, 15, 32 and 12, 13, 31 were 
averaged by van Rhijn. Mean values of / and b were taken for them (1 =175°, b = + 59.7 in the first case, and 

= 119°,b=+55.1° in the second case). 


VAR Ee 


Kharadze Van Rhijn 


| 
34. 35, 36, 37, 38,39, 40 | 38, 39, 40, 87, 110 
‘ SA ea i7, 48, 419, 41, 42, 43 |.6, 7, 14 15—82, 17, 18, 19, 44, 42, 43 
ili SA | 1.2, 3, 4, 5, 8,9, 13, 20, 24, 22, 32) 1, 2, 3, 4, 5, 8, 9, 20, 24, 22 
IV SA | 40, 41, 12, 23, 24; 25; 26; 27, 28, | 10, 11, 1213-31, 23, 24, 25, 26 
29, 30, 34 
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Table 2 and Figure 1 show statistical results along the z-coordinate, based on Kharadze's data. In the first 
column are shown the values of z over which the average of kp was taken, in the second column — the mean 
value of z, next — the corresponding values of KP and their logarithms, the number of points n, and the mean 
distance from the sun in the galactic plane Txy- Figure 1 shows log KP asa function of z. As can be seen, at 
small values of z the fall in the density may be quite well represented by the barometric law. At z = 100-200 
parsecs a pronounced slow-down in the fall of the density occurs, and from this region onwards the points may be 
fitted approximately with another exponential having a smaller slope. In the absence of extensive and reliable 
data corresponding to larger values of z, one observes a strong scatter in the values of log KP but, all the same, 
the discontinuity at z = 100-200 parsecs may be inferred from the general behavior of the distribution based on 
Kharadze's data. The following Formulas (5) represent the data given in Table 2 and Figure 1 (D = KP). 


I. D = 22.39-1074e-#1725 (g < 110 ps); I. D = 54.90-10~4e—2155-6 (2 < 95 ps.): 
D = 4.90-10~4e—#1-700 (z ~- 110 ps ); D = 12.30-1074e-#768 (2 > 95 ps); 


(5) 
I], Dis 29.17 - 10 4e— 20805 (2 250 ps); IV. D=30.00-104e-2!85-5 (¢ < 150 ps ); 
D=s 6.93> 10 4e—711685 (2 > 250 ps); D = 9.50.10 4e—*!8125, (z => 150’ ps) 
lo Eze (p 500 1000 1500 
20 


Q0 10 20 
aT 
20 10 
1 U0 o- 
0 500 1000 500 Z ( parsecs) 


Fig. 1. Distribution of log KP with z (from the data of Kharadze). 
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It is clear from the above that an application of the method of Parenago [7, 8], which assumes a single 
exponential law for all values of z, gives an average value of the gradient of the distribution over all z, while 
with the corresponding grouping with respect to the galactic latitude one should obtain the dependence of the 
gradient on b. Such a dependence was found by Kochlashvili and Toronzhadze [9],who used a somewhat modified 
form of Parenago's method. To verify this, we split the material in group IV into three parts with respect to the 
latitude (b < 20°, 20°< b <= 40°, and b> 40°) and obtained the values of KP. and z for each subgroup. Because 
of the subsivision of the data, one now has a larger scatter of points, but as is seen from Figure 2, no connection 
between the density gradient as a function of z, and the latitude b, exists. In what follows, no such grouping of 
the material was found necessary. An attempt to discover the dependence of the characteristics of the distribu- 
tion of density on the distance from the sun by suitable grouping of the data was also unsuccessful, 


The problem of the existance of changes in the density gradient at large values of z had to be examined. 
Van Rhijn's data, which are undoubtedly the more reliable, were used in this. Table 3 and Figure 3 show the 
corresponding results of calculations. As can be verified, the dependence of log KP on Zz at different longitudes 
is sufficiently well represented by straight lines. The distributions can be represented by the following equations: 


I. D = 18.41.107~4e—71249-5; JIT. D = 21.88. 10~4e—z)170.3; 
Il. D = 27.54-10-4e-2/165-8- TV. D = 20.89. 10 4e—z/221.-6, 


(6) 


Fig. 2, Distribution of log KP with z for different galactic latitudes. @b < 20°; 
O20 < b< 40°: (p> 240": 


The Comesponding formulas using values of log KP corrected for the gradient of the decrease of density 
towards the anticenter of the galaxy (see below) are given below. 


I. D = 18.84. 10~4e-21205-8, TI]. D = 23.17. 10~4e—212238 


II. D= 29,51. 10° 4e-2/154.9. TV. D = 26.00. 1074e—-21278-4, 7) 


It is very diffi i i 
ee a ca to GISEONeE a systematic dependence of the distribution of density on direction. For this 
i nl ifferent directions were averaged and the magnitudes obtained (Table 4) were treated using 
st squares, The final results in Table 4 (last two columns) contain a correction for the mean 
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Fig. 3. Distribution of absorbing matter with z using van Rhijn's data, O) IZ = 0-59"; 
®) Wl =59-81°; O) Wil = 81-118"; @) IVZ =118-143°; +) mean. 


distance from the sun in the plane of the galaxy, which depends on the density gradient along the direction 
towards the center of the galaxy. The following formula represents the distribution of absorbing matter as a 
function of z in a final form: 


D = 24.38 .1074e—21216-5+15.1 


8 
+0.68 2 


From a consideration of diagrams based on van Rhijn's data, the starting point of the theory of absorption 
developed by P. P. Parenago [1, 8] is confirmed, provided that the observational material is not subject to 
systematic errors. The numerical value of the height of the "homogeneous atmosphere” turns out to be twice 
as big as that assumed by Parenago. This means that the decrease in density takes place at a slower rate than 
was assumed up to now. Using Formula (4) we can consider the mean absorption coefficient in the photographic 
region in the galactic plane (a9), With rg—1ry = 1, A(t) A (ry) = a9 and ap = 1.086KP. According to van 
Rhijn's data KPo = 24.4-10~ stellar Magnitudes per parsec, which corresponds to 2''65 per kiloparsec. Accord= 
ing to Kharadze K Py = 30.2> 1074, i. e., 398 per kiloparsec. Kochlashvili and Toronzhadze obtained 395 for 
all areas, and 31.4 with the exclusion of areas in the zone of the galactic equator. 
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TABLE 4 


Distribution of Densities with Z According to Van Rhijn's Data (Average d for 
all Directions} i 


Corrected 
Z xo - 10! log xo + 10! ie Tey 
xo - 104 hoger. 10 
2.4 22.58 4.354 426 264 24.95 4.397 
16.3 21.29 4.328 56 184 23.61 1.373 
35.9 17,27 4.237 69 222 19.39 1.288 
72.4 13.13 4.418 72 280 1524 4.182 
142.7 10.58 4.025 65 454 13.26 4°492 
240.8 6.09 0.785 27 612 9.08 0.958 
341.2 4.83 0.684 16 714 6.68 0.825 
446.6 1°35 0.130 8 678 2.18 0.338 
543.8 0.82 1.914 6 763 4,53 0.185 
679.3 0.73 1.863 a 758 4.A9 0.076 
865.2 0.18 1. 255 4 837 0.29 1.462 


Distribution of absorbing matter in the plane of the galaxy. Six areas(SA 110, 87, 19, 8, 9, 24) with 
galactic latitudes 0, 17, 80,92, 107, 128° are correspondingly disposed within + 3.1° of the galactic latitude. 
The values of KP according to van Rhijn's data were reduced to the galactic plane using Equations (6). Knowing 
the distance of the sun from the center of the galaxy Ro (assumed equal to 7.2 kiloparsecs), one can calculate 
for different Ixy = (tz + 13)/2, corresponding to the mean densities, the distance of each area from the center of 
the galaxy using 


cy == Ty + RG — 2reyRy cos (Ll — 1p). (9) 


Here 1 — 1 is the difference of longitudes of the given direction and the galactic center (1 = — 32°), Assuming 
circular symmetry in the distribution of matter one can obtain mean values of KP for equal distances from the 
center, Table 5 and Figure 4 give results of these calculations. The density of absorbing matter as a function 
of radius may be represented by 


Dr = 19.28- 4Q~4e—(Ray—7200)/1010 | ne 


The final distribution of density as a function of z was corrected using this formula and is given in 
Table 4. 


Using the law found above, one can calculate the density of absorbing matter at the center of the galaxy. 
With Rp = 7.2 kiloparsecs the absorption per parsec is a9 = 2 83 which is obviously an overestimate. It is possi- 
ble that the distance of the sun from the galactic center is much too large and it is also possible that our law is 
not sufficiently accurate for extended extrapolation. If one assumes Ro = 5 kiloparsecs then ao is il and with 
Ro = 4 kiloparsec ap = 0.11. These magnitudes are more reasonable although for the full UE in the plane 
of the galaxy one still obtains values which are somewhat high: 5717, 549 and 222'™ respectively. 


Undoubtedly in any determination of the total mass of absorbing matter in the galaxy, one must take into 
account the decrease in density both with the z-coordinate and the radius. An integration over an infinite 


ellipsoid of rotation using Expressions (8) and (10) gives 
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Fig. 4. Distribution of absorbing matter along the radius vector 
in the plane of the galaxy (according to van Rhijn's data). Rxy is 
the distance from the galactic center. 


TABLE 5 


Distribution of density as a Function Distance from the Galactic Center 


Range of Range of ae bl, 

Re, Re, Ryy n |logke - 104 
6.5—6.7 39.00 6.602 2} 1.594 7.3—7.4 21.86 7.347 112] 4.339 
6.7—6.8 32.85 6.760 PAL abel?) 7.4—7.5 47.95 7.455 TI Ay254 
6.8—6.9 PANU 6.864 Siete 7.5—7.6 13.60 esi 4| 4.41384 
6.9—7.0 23.70 6.966 Sion V6—1.7 41.22 7.646 4} 1.050 
7.0-7.4 21.18 7.055 Ge io2e 7 1—7,.8 8.77 7.779 3] 0.943. 
7.1—7.2 48.20 7.139 8] 41.260 7.8—8.0 8.27 7.934 3} 0.918 
7.2—7.3 25.20 7.256 {19} 4.4014 8.0—8.2 412.00 8.078 2 leeOLo 

M = 28,8 Rp, (11) 


where Po is the density at the galactic center. If one assumes that the smallest possible value of the density of 
absorbing matter in the solar vicinity is p = 1007 g/cc, the density at the galactic center turns out to be p = 
Seow 0: g/cc, using Formula (10) in which p is taken instead of kK Po and the gradient Bp remains con- 
stant. Substituting the values of P9 and Bp, Bz from Equations (8) and (10) into (11) we find that the total 
mass of absorbing matter in the galaxy is 


M = 6.28-216-10102.1.32- 10722. 
-(3.08- 108)9g~5.10%g~5- 108mg 


40 


‘ 
ne Oe, 
BA OS. Beant 


—— 


500 DS 7000 DS Tzy 1800s 


Fig, 5. Densities «p for various direction in the | 


Galactic plane [ the excess Alogkp above the 
values predicted using Formula (10)]. 


A larger value is obtained if dispersion is taken into 
account. Asa final value 10° m@ May be taken as the 
estimate of the mass of dust in the galaxy. This is 
higher by one order than the value accepted so far 
(Fesenkov's estimate is 10° me{10]). 


It is possible that at some distance from the 
galactic center the rate of increase in the density 
decreases and as a result our estimate is too high. 


Spiral structure of the galaxy in the solar vicinity. 
It is best to begin the study of the disposition of spiral 
arms with a study of the distribution of stars of early 
spectral classes, especially B and A-stars which are 
sufficiently numerous, can be observed at large distan- 
ces, and have a tendency to concentrate in spiral arms 
and thus bring out the spiral structure of other galaxies. 
The distribution of other objects can be compared with 
the distribution of B-stars in order to establish their 
relation to the spiral arms. At the same time, studies 
of spiral arms in our galaxy are made rather difficult 
by the presence of major fluctuations in the distribution 
of B-stars. 


A preliminary study of the distribution of B-stars 
was made using the Bergedorf Spektral Durchmusterung 
[11] which gives star counts for a given stellar magni- 
tude for different Kapteyn areas. The simple formula 


A(m) = 1.015 -10~eD (r) (73 — 73) (12) 


can then be used to study stellar density to within a 
constant coefficient (which depends on the dispersion 
of absolute magnitudes and luminosity function), at 
least for areas in low galactic latitudes. Here w is the 
solid angle; ty, fg are the distances corresponding to 
stars with some mean absolute magnitude {in our case 
it was assumed that for BO — B5 stars M = —3' 4. and 


absorption was taken into account using Parenago's tables [8]); A(m) is the number otf stars of stellar mag- 
nitude between m and m+1; D(n) is the stellar density, In order to remove the effect of factors not taken 
into account in (12), it is best to study notthe densities themselves, which for the above reason may be 


subject to fluctuations but the maxima in the density distribution along the line of sight. It is known that 
maxima in the distribution of interstellar hydrogen obtained from radio observations [12] show a tendency to dis- 
tribute themselves along certain curves which mark the spiral structure of the galaxy. By plotting the maxima in 
the distribution of B-stars it is possible to compare the distribution of B-stars and interstellar hydrogen. In Figure 
6 the dotted line shows the positions of maxima in the hydrogen densities, the open circles indicate the positions 
of maxima in the distribution of B-stars, and open squares give the positions of regions of of O-associations which 
are also the clearest maxima of distribution of O- and B-stars. The number of areas employed was found to be 
too small for a study of the distribution of absorbing matter by this method. However, a eompa spe of bilgi 

6 and 5, showing the distribution of absorbing matter along the line of sight for the a which Were inves- 
tigated (in low galactic latitudes), allows us to draw certain conclusions as to the distribution of absorbing matter 


relative to the spirals in the galaxy. 


_-_-— 


-——--—~— 


150 hike 
@ 


galactic 
center 


es 
-—- ee 


= 


Fig. 6. Distribution of maxima of stellar densities of B-stars and the densities of inter- 
stellar hydrogen. Notation: dotted line — distribution lines of maxima in the density of 
hydrogen (after Parenago), @ and ©) maxima of B-stars (open circles are less reliable), 


QO) areas of 0-associations according to Markarian [13]; continuous line — possible posi- 
tion of spiral arms. 


It follows from the above graphs that two spiral arms are brought out in the neighborhood of the sun by the 
distribution of B-stars: one of these approaches the sun to a distance of about 1000 parsecs and is directed towards 
the anticenter of the galaxy, and the other approaches the sun very closely (200-300 parsecs) from the direction 
towards the galactic center. Maxima in hydrogen density confirm the tendency of B-stars to distribute them- 
selves along these lines and to coincide with them in direction, although absolute coincidence is not observed. 

In spite of the errors involved in the determination of distances to B-stars and to hydrogen maxima (in the first 
case absorption was taken into account using Parenago's data and with B> = 100 parsecs, ap = 3.5 per kiloparsec 
which could lead to systematic errors in the distances, and in the second case, possible errors may be due to 
differential rotation of interstellar gas) the distribution of interstellar hydrogen may be much more "diluted" than 
the distribution of B-stars. If the spiral structure does not exist in the neighborhood of the sun, the hydrogen may 
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"fill in” vacancies in the distribution of B-stars and 
penetrate into areas where B-stars are not observed. 
Undoubtedly, only a simultaneous study of B-stars and 
I the interstellar medium can finally provide a solution 
| 30 to the problem of the location of spiral arms in the 
galaxy. 


Certain regularities are noticeable in the distri- 
bution of absorbing matter relative to the outlined 
450 spiral arins. The distribution of absorbing matter along 
some directions from the sun is plotted in Figure 5 
(values of A log kp). Anticentral directions (2 = 128, 


0 107°) have a general tendency, which is at once noti- 
ceable, of giving large absorption in the neighborhood 
foe. a of the sun, and then a slow fall-off reaching a minimum 
+as oe aonn - (at distances of 500-700 parsecs) at places where one 
2 OSS would expect a gap between two arms. It follows that 


a0 the absorbing matter which accompanies spirals is more 


"smeared out" in space. Dark matter wraps itself round 
the spirals. This is confirmed by the distribution of 
along Z = 0° where the line of sight enters into the arm 
nearer to us. Here the density increases as the line of 
sight enters the spiral (the data are free of the effect 

of the general density gradient along the radius). The 


+25 


a0 


Fig. 7. Distribution of absorption (continuous data for other directions confirm this. Whenever the 
and dotted curves) and luminous matter (photo- line of sight passes through the gap between the arms 
graphic luminosity) (circles) in the spiral arms (Z = 17°) the density slowly decreases (depending on 
of NGC 5194 (for different sections) (after the distance from the inner arm), and whenever it 
Holmberg). approaches the inner arm to a distance of 600-1000 


Parsecs an increase in the density is observed. Thus 

the distribution of the excess deusity of absorbing matter 
above the regular decrease along the radius vector confirms the general picture of the disposition of spirals in the 
neighborhood of the sun. In this connection one must mention the work of Holmberg [14] who studied absorption 
of light in the arm of the spiral galaxy NGC 5194 whose spiral screens the irregular galaxy NGC 5195. As can 
be seen from Figure 7, which was taken from [14], Holmberg found that the distributions of luminous and absorb- 
ing matter do not coincide with each other, the maximum of the distribution of dark matter is shifted relative 
to the arm towards the center of NGC 5194, and, finally, the absorbing matter occupies a larger region of space 
than the spiral which consists of stars. An analogous picture is apparently observed in our galaxy. Such a dis- 
tribution of absorbing matter may be explained by radiation pressure which acts on the absorbing dust particles, 
the pressure being due to gigantic hot stars which fill the spiral arms of the galaxy. A random distribution of 
stars and interstellar matter is superimposed on the spiral structure and to a considerable extent masks it. For 
this reason, we cannot expect a very clear picture of the location of spiral arms in the solar neighborhood but 
should try to extract it from a chaos of randomness, One should not reject the possibility of the existence of 
spiral arms in the solar neighborhood solely because it is only brought out by the general tendency in the distri- 
bution of Bestars, interstellar hydrogen, and dark matter. 


Spirals may turn out to be connected by "bridges" consisting of random clouds of stars and interstellar matter 
and located between the arms. Such a "bridge" may possibly be destroying the regular behavior of maxima of 
B-stars and interstellar hydrogen in directions from the sun towards 7 = 60°. In order to obtain a detailed picture 
of the location of spiral arms in the neighborhood of the sun one requires, undoubtedly, not only fuller data on 
the distribution of B-stars, but also the distributions of stars of other types, above all stars of spectral class A, 
which are more numerous and,with sufficiently powerful instruments, can be seen up to distances of 2-3 kilo- 


parsecs. 
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ON THE DISTRIBUTION OF MASS IN THE GALAXY 
A. M. Mikisha and F, A, Tsitsin 


Lindblad’s method for determining the relative mass of the galactic 
nucleus is analyzed. It is shown that a representation of the distribution 
of mass in the galaxy by a uniform model leads to more than a tenfold 
overestimation of the relative mass of the nucleus. A model of the 
galaxy is set up which takes into account the variation of flattening of 
the surfaces of equal density with the distance from the galactic center. 


1, A method of studying the mass distribution in the galaxy has been given by Lindblad [1]. Themethod 
was used by a number of authors (Oort, Titus, Chandrasekhar and others), and is described, in its most highly 
developed form, by Parenago [2] in his "Course of Stellar Astronomy." 


Since neither "solid rotation” (uniform spheroid, force in the galactic plane directly proportional to the 
distance from the center) nor "Keplerian rotation" (all mass concentrated at the center, force in the galactic 
plane inversely proportional to the square of the distance from the center) can be associated with the galaxy, 
- Parenago [2] concludes that "in reality, part of the mass is concentrated at the galactic center and the rest is 
distributed throughout the galaxy (in the first approximation, uniformly).” 


In that case ". . .qne can try (our italics, A. M., F. Ts.) the following form of the force of attraction as 


” 


a function of distance 


Fa Fy, 4+ Fee oh? + BR. (1) 


(cf. [2] Section 34). 

It must, however, be pointed out that this form cannot be taken as empirical. Since, in the first approxi- 
mation, a uniform spheroid with a nucleus is taken as the model, the force on a point mass in the equatorial 
plane can be obtained quite accurately, and is given by Equation (1). Thus Equation (1) is introduced not because 
one can associate with the galaxy some sort of intermediate state of rotation between the "solid" and the "Kep- 
lerian” but because of the particular choice of a model, 


In later calculations Lindblad and, following him, Parenago too, depart from the adopted model and,in 
fact,calculate not the mass of a uniform spheroid the the mass of a uniform flat disc, 


A correct formulation of this problem (i. e., the determination of the mass of a uniform spheroid) is given 
by us below. 


Equation (1) has often been criticized (Bottlinger [3], Parenago [2]). The main point of this criticism is 
the following. According to Equation (1) the force of attraction is inversely proportional to R? for small values 
of R .and directly proportional to R for large values. This is contradicted by observational data which show that 
at small eigtenees from the center the force is proportional to R. At large distances the force should, of course, 


fall off as R72. For this reason Equation (1) appears to be unsuitable. 


This criticism is unsound since it is based on an erroneous interpretation of Equation (1). 


Actually, Equation (1) gives an expression for the force acting on a point mass in the equatorial plane of 
a uniform spheroid with a nucleus. However, in its derivation, it was assumed that the nucleus attracts the point 
mass as a Newtonian center and, consequently, the point is situated outside the nucleus. On the other hand, 
the attraction due to the uniform spheroid was calculated for a point inside it. Thus Equation (1) is only true 
for points lying outside the nucleus but inside the spheroid and R cannot be made to tend to zero or infinity. 


Rejecting Equation (1) Bottlinger [3] proposed the following expression for the force of attraction in the 
galactic plane: 


(2) 


P, P. Parenago criticizes this expression and writes: ".. .this ingenious expression has no mechanical 
Dasis: (2). 


However, it is easy to show the mechanical significance of Equation (2). For example,* in the case of a 
spherically symmetrical system, Equation (2) results in consequence of the following density distribution: 


_ 3a { 
0) = ree BRA 


It would be interesting to consider a nonspherical distribution of mass giving rise to the law of force given 
by (2). 


In spite of the fact that Bottlinger's criticism of Equation (1) is based on a misunderstanding, this form of 
the law of force must still be rejected since, as will be shown below, it corresponds to a model whose adoption 
leads to an overestimation of the relative mass of the galactic nucleus, 


The immediate cause of our investigation was P. P. Parenago’s statement in which he formulates the main 
result brought out by the above method of determining the mass of the nucleus and the whole galaxy: “Essenti- 
ally, by determining J2t 1 (mass of nucleus, A. M., F. Ts.) one determines the mass distributed through galaxy 
within the distance between the sun and the center of the galaxy (this would be exact in a spherically symmetri- 
cal galaxy)” ([2] Section 34). Again, in another work, Parenago states: ". . .using Equation (12) (which corres- 
ponds to our Equation (1), A. M., F. Ts.) we determine not the mass of the central nucleus but the approximate 
mass contained within the solar distance (this would be true exactly in the case of a spherical galaxy)" [4]. 


However, it is easy to see that what is actually determined is the whole mass within the solar distance, but 
only if the law of force F = aR”? holds (this is exact in the case of a spherical distribution of matter). Conse- 
quently, an analogous statement about the law of force given by Equation (1) is incorrect. Hence it follows that 


the large magnitude of the relative mass of the nucleus is largely due to the adoption of a uniform model of the 
galaxy. 


To substantiate our point of view we shall consider the consequences of the adoption of a uniform model 
for nonuniform initial distribution with spherically symmetrical and spheroidally symmetrical distribution of 
matter, This method has the advantage that the initial distribution in each case is known, which is not the case 
when a model approximating a real galaxy is set up. 


Before proceeding any further we must note the following. The initial distribution is fully determined by 
a finite number of parameters, all other parameters being expressible in terms of them. We shall take the 
numerical values of the parameters of the initial distribution to be equal to the numerical values of the corres- 
ponding parameters of the real galaxy. A comparison of the remaining parameters with the corresponding 


*We are indebted to I. L. Genkin for this idea, 


46 


parameters of the real galaxy will, in addition to the main result, allow us to estimate how closely the initial 
distribution chosen by us corresponds to the real galaxy. 


2. Let 6 be the density of the nucleus, r the radius of the nucleus, Ro the distance of the sun from the 
galactic center, Rg the galactic radius, G the gravitational constant, IM ;the mass of the nucleus, M the mass 


contained in a sphere of radius Ry JM mass of the galaxy excluding the nucleus, and IR'the mass of the galaxy 
(including the nucleus), 


We shall assume that the values of Oort's constants A and B at a distance Ry from the center of the galaxy 
are known. This is equivalent to a knowledge of the force Fy at a distance Rg, and its derivative F'p. 


We recall that expressions such as "distance of the sun from the center of the galaxy”, “radius of the galaxy,” 


"galaxy" etc., mean only that we are assuming numerical values for the parameters of the initial distribution 
which are equal to those accepted for the reai galaxy. 


The initial distribution, A spherically-symmetric, nonuniform system is fully defined when the density 
of matter p(R) is given. Without assuming the form of Pp (R) we can set up a system of equations which, using 


the given parameters of the initial distribution (r, Ro, Fo and F'g), will determine the remaining parameters. 
Let 


| Pp ae o(R, 1, Xp), (3) 


where ky and Kg are the parameters to be determined. 


We can then write down the following equations 


0 (Ro, %1, %2) = Po; 


P(r, 4, %2) = 3, (4) 
where Po (density of matter at the distance Rg) and 6 (constant density of the nucleus) are also unknown. 
The force of attraction at the distance Rg is given by 
GM, GM 
Py = R? ap R? (5) 
0 
where 
Vi en 
My, = aa Treo 
Ry (6) 


M =4n \ (R, 4, %) dR. 
: (7) 


Taking into account (6) and (7) we can find expressions for Fo and F"p in terms of the required parameters 


Ro 
4 r8 4nG 9, l 
sapere teas eae aay ap = reel 
Kocoy Gnas yt Fe o (RR, *1,. *) 
Meas (8) 
8 Oy ryine ee etl Re kis ka) IE 
Fy = — $ 7G (z-) 8+ RGPo R p ( 1 2) 
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The set of Equations (4) and (8) are a set of four equations involving four unknowns (P9,6,K4 and K 2) which 
with a given form of the function Pp(R, Ky, K2) can in principle be solved without difficulty.’ 


The interpretation of observational data given by B. V. Kukarkin and P. P. Parenagoshows that the expon- 
ential distribution of density is in good agreement with observation (cf., for example, [2] Section 50), We shall 
therefore take 


p(L?, ¥1, %) = xe, (9) 


Let r = 1.5 kparsec; R = 7.2 kparsec; RG = 13 kparsec; 
A = 19.5 km/sec/kparsec; B =— 13 km/sec/parsec. 


On solving (4) and (8) we find 


Mm, = 0.61-10!° sun masses 
NY = 10.7-10!° sun masses 


i. e., the mass of the nucleus of the initial distribution is ~ 5.7% of the whole mass of the distribution, 


The nonuniform model. We shall now consider what will be the effect on the relative mass of the nucleus 
Me of the adoption of a uniform model for the nonuniform initial distribution under consideration. In this 
case instead of (6) and (7) we have 


4 
M = = xrs; (10) 
M = > «(R3—r?) p, 
(11) 
where P is the density of the uniform mass outside the nucleus. 
Equation (8) now takes the form 
Rs eo 3 Sagi pe 
9 3 2 Ouse a) P; 
-f 8 r\3 4 4 oro 
Fo=—5G(EV8 + S06 (t 4% a 
3 R, at 3 bE re )P 


Solving Equation (12) with respect to 6 and P we find finally the following values for the mass of the nucleus 
JI, and the total mass of the modek 


M, = 0.72-101 sun masses 
Mt = 1.79-10" sun masses 


i. e., the mass of the nucleus is 40.5% of the total mass. 
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Comparing the magnitude of the relative mass of the nucleus INW,/N - for the initial distribution with 
that calculated for an initial distribution using a uniform model , we see that the adoption of a uniform model 


for a spherically symmetrical nonuniform distribution leads to an overestimation of the relative nuclear mass 
by a factor greater than seven, 


It is easy to see using Equation (8) that in the case of a spherically syrnmetrical distribution of density 
we have 


dF (R 
2F (R)-+R 
BASH Saint SiaeRione. fe Sa. 
The latter may be transformed into the following form: 
B— A)(A+ 3B 
(hh Se as (13a) 


(Equation (13a) may be convenient in structural studies and the dynamics of globular clusters and galaxies of 
type EO.) 


Using Equation (13a) with R = Ro we have 


p (Ro) = po = 7.96- 1075 g/cm? 


This value of Po is smaller by a factor of 7-8 than the observed density in the vicinity of the sun: Po = 5.5— 
—6.6:10°" g/m? [5]. 


This discrepancy shows that the density distribution in a real galaxy is far from spherical. 


din p 


aR . In the initial distribution it is equal to 0.63 kparsec7* 


The quantity k2 in (9) is, clearly, — 


92 IORI be 0.27 kparsec™*. 


which leads to the following value for the logarithmic gradient of density: aR 


; | dl 
Observations give the following values for — es ([2] Section 50): 


0,25-0.27 kparsec™! (spherical component), 


0.11-0.16 kparsec™* (plane component). 


The agreement between the numerical value of the logarithmic density gradient in the galactic plane and 
~ that observed for stars of the spherical component leads us to the conclusion that the dynamic effect of the spheri- 


cal component is far greater than is usually supposed. 


3. Let the initial distribution be a nonuniform spheroid. Let us find the magnitude of the relative mass of 
the nucleus MN, /M. '.. We shall consider a nonuniform spheroid with an exponential density distribution in 
the equatorial plane in accordance with Equation (9). If one assumes that the surfaces of equal density are similar 
spheroids then it is easy to find the value of the density at any point in the body, The nucleus will be assumed 
to be uniform. Let r = 1.5 kparsec; Ro = 8 kparsec; A = 19.5 km/sec/kparsec; B = ~ 13 km/sec/kparsec. 
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The following notation is employed: r is the radius of the equatorial section through the nucleus; Ris . 
the radius vector in the equatorial plane of a spheroid; M is the mass contained within a spheroid with semi- 
major axis Rg; e is the eccentricity of the meridian section through the spheroid; ¢ is the compression 
(e= Vi-e). 

We emphasize once more that we are not attempting to model the galaxy, although the choice of initial 
data may suggest this. We are considering a fully defined initial distribution of matter with the object of dis- 
covering what errors are introduced by assuming a uniform model for it. 


Undoubtedly, the initial spheroidal distribution under consideracion is nearer to the structure of the real 
galaxy than the spherical distribution, although it does not appear to be in the best agreement with observational 
data on the galactic structure. For example, it does not take into account the observed compression of the 
nucleus [7], and it excludes the possibility of taking into account changes in compression in the galaxy. 


We can get an idea of how near the given scheme is to a real galaxy from the numerical values of the cal- 
culated parameters (e. g., compression € and the logarthmic density gradient in the galactic plane), which are 
sufficiently well known for the galaxy and are not used in the setting up of the model. 


Initial distribution, We shall now proceed to determine the main parameters of the initial distribution: 
Ky, Kg, 6 ande. Having obtained these parameters we shall find the relative mass of the nucleus M,/M. 
Using the calculated value of € we shall obtain the uniform model of our initial distribution and the relative 
mass of the nucleus for this model, and hence we shall be able to find by direct comparison the errors intro- 
duced into the nuclear mass by the choice of the uniform model. 


Since the compression € is taken as unknown in the initial distribution, we shall assume that the density 
at a distance Rg is known. We shall take it (using observational data on a real galaxy) to be 6.0- 10°" g/cm), 


We shall assume that r, Ro, Po, Fo and F'y are known. The system of equations from which we shall deter- 
mine ky, Kg, 6 ande may be written in the form 


= Mess) 
Po = eres 
Fy = F, (Ry) + Fa(fo); (14) 


Fo = Fy (Ry) + Fa (Ro); 


where F;(R) and F,(R) are the forces due to the nucleus and the spheroid,respectively. 


The attraction due to the nonuniform spheroid, with the nucleus removed, can be expressed (assuming 
exponential distribution of density in the equatorial plane and spheroidal similarity of layers of equal density) 
in the form 


one 
As eV rit 


Fy (Ry) = 4Gx, \ 


o=0 I=0 2=eV ri]? 


ol chil sh 


(Ry — 1 cos @) e~*# VY P+2* tay de dz 
(2? 4+ R’ + 1? — 2Rol cos 9)" 


rR, V eacztet 
beets \ \ \ (Ro — Lcos 9) e727" Tal de dz 
; co ey (2? -|- R --+ 1? — 2R,l cos 9) 2 } 


This integral cannot be evaluated completely by analytical methods and numerical methods must therefore 
be Used. ee shall evaluate this integral approximately by replacing it by a sum of integrals, Let us divide th 
spheroid into seven similar spheroidal layers (not counting the nucleus). We shall assume i density in each : 
layer to be constant, We shall retain the overall exponential law of density distribution (Fig. 1). : 


To simplify the computations we calculated the full 

per) attraction by summing up the attractions due to a series of super- 
imposed uniform spheroids, instead of calculating the sum of 
attractions due to the separate layers. The density of each of the 
spheroids is equal to the height of a "step" on the density curve, 
i. e., the difference in the density of successive spheroidal layers. 
For this reason, the density of the nucleus 6, for example, is 
represented as a sum of densities of eight superimposed spheroids. 


Gh aa, 


, ard 


A uniform compressed spheroid attracts a point mass lying 


Fig. 1. The continuous line represents Fy : : 
in its equatorial plane with a force given by 


the true “integral” density distribution; 
the dotted line represents the density 
distribution assumed in the evaluation 


of the sum, F(R) = 2nGp ae [= arc sin — yt : oes eek , 8) 


where a is the semimajor axis of the equator section through the spheroid, e is the eccentricity of a meridian 


section through the spheroid, p is the density of the spheroid. The force is assumed positive in the direction 
towards the center of the spheroid. 


ae 
Putting pois we obtain 


= e4 a 
F(R) = 2nGo era [arcsinu —uY1 — uJ, 
le f 2 
F’ (R) = 2xGop —.— Ese == |are sine ee te | : 
Vi-w 
Summing up over all the superimposed spheroids we finally have: 
an RR vise’ 
Qn Geo s (%— bi)ei, 2nGeo i > (a: — di) e; (16) 
i=] i=1 


where aj, bj, cj and dj are functions of Kz, e and other known quantities. 


The set of Equations (16) is solved by a graphiconumerical method in which the right-hand sides of (16) 
are tabulated for different kz and e. 


Using the values of ky, Kz, 6 and e we find the following values for the initial distribution: 
ee lL / TR; 


M@, = 3.47-10° sun masses, 
mM = 1.26-1011 sun masses, 


Thus the mass of the nucleus JJt, is 2.75% of the total mass of the initial nonuniform spheroidal distribution. 


Uniform model. If we represent the initial distribution by a uniform model, i. e., a uniform spheroid 
ee ee . . . . 
with the same compression and a nucleus, we obtain the following values for the masses (using the method dis- 


cribed above): 


bl 


y2, = 0.70-10" sun masses, 
= 1,24-10!1 sun masses, 


Thus the mass of the nucleus is 56.5% of the total mass. An attempt to represent a nonunoform distribution by 
a uniform model leads to a twentyfold (!) overestimation of the mass of the nucleus. 


We have shown above that in his account of the determination of the galaxy, Lindblad departs from the 
uniform spheroid by assuming € << 1 and actually using € = 0, and therefore goes over to the case of a uniform 
disc. 


It is possible, without writing down long integrals, and without assuming € << 1, nore = 0, to obtain a 
solution to thé problem in its general form. 


In fact, it is easy to show that (using the notation of [2]) 


Figg = 226001 +9 (asotant — 


3 tn) (17) 


where a (cf. [16]). 
Let 
3 (142 | 
9 (1) os a ( a} (18) 
then 
ers , | 
eS ae oe 


Substituting Fo, = 0,20Fo, which was obtained by P. P. Parenago [2] for a model of a uniform spheroid with a 
nucleus, we have 


__ 3.62-1088 
i= src nts (20) 


The total mass of the spheroid can be found from the expression 
Rg\3 
Mm =, + M (2) =M, + 5.93 M. eet) 


Let us construct a table of M(€) and IM, (¢). 


bel The case considered ee P. P, Parenagg, i. e., the case of a plane disc, is given in the last line of Table 1, 
: the latter case we have ~ Jt,/MN = 0.60, i. e., the mass of the nucleus ‘My is 60 % of the total mass of 
the galaxy. We note that the effect of the compression € on the mass of the spheroid is Rey! small, 


qo 


TABLE 1 


€ | l | @ (1) | M | Me Remarks 

1 0 4.000 3.62- 1049 gy =1.84- 101° m 410.7-1019 m Yo is found from 
® c) ; 

Ig => | 4.418 [2.55 1.28 7,6 SNe 

Wy loyo | 41.644 [2.24 1.40 6.5 

1, | Yas | 4.778 |2.04 1.02 6.0 

1, | Vex | 2.029 11.78 0.89 5.3 

0 ©o 2.356 |41.54-1043 g =0.77.101 m2 3.6-101 m ¢ (oo) = 3/4 is 

found from 


I | L"Hospital's rule 


4. The method which was used in studying the model of a nonuniform spheroid with constant compression 
can be used to construct a model which will approximate to the real galaxy. 


The galaxy is looked upon as a system consisting of a uniform spheroidal nucleus and a variable density 
spheroid. This spheroid is conditionally divided into two spheroidal parts, one inside the other, witha dividing spher- 
oidal layer passing through the sun (Rp = 8 kparsec) between them. The compression € is assumed to be variable 
in the inner spheroid and constant in the outer one. This assumption means that we need not consider the gravi- 
tational effect of masses which lie beyond the limits of the spheroidal layer through the sun. Thus € = € (R) 

(R is the radius vector in the equatorial plane) where 


o 


== @, == const for R <r, 
=aR+b fr r<ck<hy, 
== € <= const for «h > Ro: 


oo @ 


The quantities a and b are found from boundary conditions to be equal to 


= & Ro — Eq? 


Ry—r 


where €, is assumed equal to 1/2 (cf. [4] or [7)). 


Such a model of the galaxy allows us to limit ourselves to a consideration of the densities in a “cut” 


through the galaxy along the equatorial plane. Suppose that the density at some point (x, y, 2) in the galaxy is 
given. It is clear that the surfaces of equal density in our model are surfaces of revolution. We have to find 
that value of R which is equal to the radius of an equatorial section through a surface of equal density passing 


through the point (x, y, Z): 


x2 + y? Ze a 
a on ap ene oe) 


Equation (22) is of the fourth degree in R. Using the theorem of Descartes we know that there is at least 
one positive root of this equation and we shall denote it by u. It is then possible to write 


P(X, Y, 2) = me", (23) 


which can be used to determine the density at any point. 


We recall that the previous paragraphs of this paper were mainly devoted to the problem of when it is 
allowable to represent a nonuniform distribution of matter by a uniform model. It we take the initial charac- 
teristics of a nonuniform distribution to be equal to those obtained from observations we can obtain as an 
additional result an estimate of how closely the nonuniform distribution under consideration approaches the 
real galaxy. 


The aim of the present paragraph is to discuss our ideas on the structure of the galaxy and as an additional 
result obtain an estimate of the accuracy obtained when this model is approximated by a uniform spheroid. 


Thus, the quantities r, Ro, Po, €;, Fo and F'y are considered to be known. The quantities Ky, Kz, 6 and 
€9 are determined. 


Let us construct a table of the quantities M(eé ) and Ms, (e). 


Using a method analogous to that of the previous paragraph we have 


, 


8 8 
F 0 ’ 
OR ue enn ete = ee (24) 
i=1 


i=1 


where aj, bj, cj, dj, and yj are functions of K2, €9 and known quantities. We find 


%q == 0.36kps "i. e., ae — 0.16 kps~* 


€9 = 1/9,hence e€ (R) = — 0.0598 R + 0.5897where r = RX Rp); 
xy = 1.07-10°%g/em®; 8 = 0.62- 107? gem’. 
Finally we find masses of the nucleus and the galaxy to be 


IM, = 0.64-101° sun masses, 
Mt = 1.82-101! sun masses, 


i. e., the mass of the nucleus is 3.5% of the total mass of the galaxy. 


If we recall that the adoption of a uniform spheroidal model gave a nuclear mass which was 60% of the 


total mass we find that, by assuming a uniform model for the given distribution, the relative mass of the nucleus 
is overestimated by a factor of 17. 


5. No model can pretend to cover the entire range of observed facts. For this reason any approximation 
to the real galaxy must involve the assumption of a finite number of parameters, best known from observations 
The number of the given parameters must be such that the model is fully defined by them. The remaining 
parameters of the model are then uniquely determined. These parameters give new information on the structure 
of the galaxy which could not be deduced with confidence directly from observations. 


Table 2 summarizes the main results of the present work. 


aes ER ne Table 2 shows that the uniform model cannot be used, even in the first approximation 
or studying the distribution of mass and, consequently, the dynamic characteristics of the galaxy. 
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TABLE 2 


No. Model Log. density | Mass ofnucleus | Total mass of |Relative mass Overesti- Compress- 
gradient MM, (sun masses)| galaxy Mt  lofnucleus mationof |ion « 
(sun masses) | a = (Nty/M) | the relative’ 
» mass of the 

nucleus by 
assuming a 
uniform 
model 

F “nif /a) 

1 Nonuniform spheres = ara) 0.61- 10% 1.07: 10% 5.7% 7 times 1 
with spherical uni- 
form nucleus 

2 | Nonuniform spheroid — 0.20 0.35-10"° 1.26- 1044 2.8% 120.5 times Wey 
with a spheroidal 
uniform nucleus 
(constant com- 
pression) 

3 Nonuniform spheroid — 0.16 0.64-10" 1.82.10 3.5% 17 times € (R) = 
with a spheroidal =aR +b 
uniform nucleus forr = 
(variable com- =R= Ro 
pression) €(1)="4 

€ (Ro) = 
=, 


Model No. 3,in which the initial parameters represent most of the observed characteristics of the galaxy, 
allows us to estimate the masses of the galaxy and its nucleus, the fall-off in the density in the plane of the 
galaxy and changes in the concentration of matter in the galactic plane as a function of distance from the 
galactic center. 


Shternberg State Astronomical Institute Received July 9, 1956 
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PERMANENT CONFIGURATIONS IN THE PROBLEM OF FOUR 
BODIES AND THEIR STABILITY 


V. A. Brumberg 


The stability of permanent configurations in the problem of four bodies is inves- 
tigated for some particular cases. Section 1 contains an analysis of particular solutions 
in the problem of four bodies based on the work of Andoyer, Meyer, MacMillan and 
Bartky, Subbotin and Moulton. A solution of the problem of stability of circular orbits 
in the problem of N bodies is given in Section 2. Section 3 deals with the stability of 
particular solutions of the four-body problem in some special cases. Finally, an attempt 
is made to apply the theory of permanent configurations to the Trapezium of Orion. 


We know that in the problem of three mass points which interact according to Newton's law there exist 
only two premanent configurations, the straight line and the equilateral triangle. In circular motion the first 
of these is unstable (in Liapunov's sense) for all mass. values but can possess some conditional stability. The 

necessary condition for the stability of the second configuration is the inequality 


MyM, + Mgmg-+ Mam; __ 4 


(my -- my + mg)? ~ 27 (1) 


The stability of permanent configurations in the problem of four bodies is of interest and is the subject 
of the present paper. Despite the simple character of the investigation in principle, it could be carried through 
for only a few particular cases. 


In Section 1 Lagrange's solutions in the problem of four points are analyzed on the basis of papers by 
Andoyer [1], Meyer [2], MacMillan and Bartky [3, 4], Subbotin [5] and Moulton [6]. Section 2 discusses princi- 
pally Andoyer's results [7] concerning the stability of circular orbits in the problem of N bodies. In Section 3 
stability in some very special cases is investigated. Finally, an unsuccessful attempt is made to apply the theory 
of permanent configurations to the Trapezium of Orion. 


1. Permanent Configurations in the Problem of Four Bodies 


The equations of motion of the problem of four bodies in a barycentric system of rectangular coordinates 
with the rotation components p, q, r_along the x, y, 2 axes is as follows: 


of 


z,—2ryy + 2qe — (g? +72) + (pg —T) yi + (Pr + 9) A = 


yi— 2pas + Ara; — (r? + p*) yi + (Gr — pyat(qp+r)u = 
U4 


Sm Y; 
r] 


ij 


t — 298; +279, —(P EP) a+ (rP— Qt (rd +P) HK = 


' j 
= m; ’ 
pr wert (1.1) 


where rj; = (11 — a3)? + (yi — yi? + (21 —2))",, Uae ty. a, ae ee t=47 (the gravita- 
tional constant is taken as unity). 


A permanent configuration is one in which the ratio of the distances between bodies is constant during the 
motion, i. e., the following condition is satisfied: 


L; Yi 24 0.0 (0 
So eS SERS CF Xi, Yi, 2; = const. 
af y? 29 e( ); vy Yi, v (1.2) 


Equations (1.1) reduce to the following: 


9 (p — (g? + 7?) p] + y8 [— 279 + (pq —r) 0] + 24 [290 + (pr +9) 9] = 


0 0 
Le — 2; 


4 , j 
= —— a 
oD) i ae 


j _ 


at [2ro +- (gp +7) p] + 92 fe — (r? + p?) ep] + 22 [— 2pp + (ar — p)p] = 


0 ( 
J 


0 
A eceeitan Voce Bs (1.3) 
ner ail pater 
j 7 


as [— 2g6 + (rp —q)p) + y [2pp + (rg + P) p] + 23 [p — (p? + 9°) 0] = 
29 — 2° 


1 , j 
=a 703 


5] os 


There are three possible cases: 


I, All points of the system lie on a single straight line. For this line we take the x-axis and choose the 


y-axis so that q = 0, From (1.3) we obtain 


p=—-~+5 (p=0), (1.4) 


and the necessary condition for the existence of a ‘straight-line configuration is 


eg Qn 
>) my —a- = — nat, n= const. (1.5) 
j "ij 


Thus, when c # 0 all points describe similar Keplerian orbits (with equal eccentricities) whose focus is 
at the center of gravity of the system, and the points simultaneously occupy similar positions. Equations (1.5) 
for the problem of N bodies were first investigated by Moulton, who showed that there are’ + N ! solutions. 


Il. The points of the system do not lie in a single plane. Since each of the coefficients of x’, y°, z° 


in (1.3) must have the form const/p? we easily obtain 


=p—h, a, B, y= const, 


0 0 
Obese od 1 , v. — x; 
4 [pp — pir®] — yi. ( BY i, 
‘ j "iy 
od or” 1 yo — yf 
, ; Fi 
Serle Rams IA ioral idl toe Stic es oe 
j y 
0 
0° 4 , oy 25 
ria are 
raj 


i 
Consequently, 


Fe (Pr) Symi (x? + y%) = 0, 
t 


i. e., p*r = const and p°r? = const. A case with P = const, r = const is impossible because we would then have 
v4 
, j i a2 
); er 0. 
j ad 


Thus r = 0, f =—n?/p? (straight-line Keplerian motion) and the necessary conditions for the existence of a 


solution are 


0 
x x 
ald j 2,0 
> PES EE cal 
J 
0 
, Ua sud 20 
Win nyt, 
ij | 


of 


or in vector form 


Sin (h fii = Sym 


Since the vectors Hj are not coplanar it is necessary that 


2 M 


0 0 0 Co 00 a 
Tyg FS M9) 44 1 on SS oe a og Wis eae 


Thus the only permanent spatial configuration is a regular tetrahedron with arbitrary masses for the points of 
he system, and motion occurs along straight lines passing through the center of gravity in the same direction 


for all points. 


Ill. The points of the system lie in a single plane but not on a single straight line. 


Taking this plane as the xy plane in (1.3), we have 


p=", g=Bp"*, (a2 + 8%)p=0. 


which is impossible. Thus p = q = 0, and the equations of motion are again given by (1.4), i. e., all four points 
The necessary conditions for the existence of plane confi- 


move in a fixed plane along similar conic sections. 


gurations are 


0 
p el=aat 
J a 2.0 
yy ™ gS OR, 
: ‘i 1.6 
0 0 I) 
reel mess, 2,0 
>) ae oe 1) 
ij 


Bigs ices 


Let us consider, following MacMillan and Bartky, an arbitrary quadrilaterial with the sides, ry, rg, tg, T4. 


The directions of the corresponding vectors are shown in Fig. 1. 


1, 2, 3, 4) the distance along rj from the point mj to the intersection with the 


We denote by ajrj(i = 
aj is considered positive when the distance is measured in the same direction 


opposite side (at point A or B); 
as ¥j; otherwise it is negative. It is easily established that 
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1 
= aS [Bors + Baral, 5, = “. [Bats + Bors]; 


1 
Togs ay [Bore + Byry], typ = ; (Bits + Byrs]; Bi =a —1. ie 


Introducing the new constant Ro = n?/M , using (1.7) and setting 


1 4 
R=-—=, R=—~, GK=R—R, Teic='itya ls: etees Foal, se ig ee 


a 

(1.6) can be given the following vector form: 

(Syma + ByS'57703) F1 + (Bg S53 — Hy yIM4) Pg = 0, 
(agSa7mg + BgSg7m4) Po + (By SqMq — X51) r, = 0, 


(%qSg7q + Bg) r3 + (B25 ,7, — %pS gq) rg = 0, (1.8) 
(otS 4m + By Sgr) ry + (BgS p72 — 0%45'gMmMs) Py = O. 


Noting that 


Ai_y Bi,’ (1.9) 
and considering that the vectors rj are noncollinear, we obtain two equivalent systems: 


S1%M_ + S5Bamz3=0, Seam, + S58. = 0, 
Syagmg + SeBgm,=0, S5x3m 3 + SyBgm, = 0, 
Sgoqm, + S,Bym, =0, Sqaging + S,8,m, = 0, a) 
Sx, + SoBimM. = 0, Sam, + Sy Bim, = 0 


from which it follows necessarily that 
S Sg == Sq, = S56 (1.11) 
and 
(S, — Sg) (Ss — S3) (Sq — S,) = (S, — 55) (S53 — S2) (Sy — So). 


Eliminating Rp from (1.10), we have the following system, in which the coefficients of the masses depend 
only on the mutual separations: 


(S, -- Sg) aM + (S5 — S's) Baits = 0, 
(Sq — Sg) XM + (Sg — Sa) Bam, = 9, 
(Sg — Sg) %IMq + (Ss — Sy) By, = 0, (1.12) 
(Sq — 55) %mm, + (S, — Sa) Bim = 0. 
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The determinant is 


0 (S,;— Sp) a (S,;— S35) Be 0 
0 0 (Sy —S5) 43 (Sg — Sa) Bs aS 
A= is; See eee gM sgh aaetaaya 
(S4— Sg) % (S, — Se) By 0) 0 


ra (S; ae iS) (Ss, a Ss) (S, ee 9) (S, aa S'4) BiB oBsB4 at 
— (S_— 5) (S4 — S5) (S11 — Sg) (S3 — Sq) Xp %q%qXq = O 


in virtue of (1.9) and (1.11). It is easily seen that the riecessary and sufficient conditions for the vanishing of 
all third order minors of A are 


which corresponds to an equilateral triangle with three arbitrary but equal masses at its vertices and a fourth 
arbitrary mass my at its center. For any other permissible quadrilateral the matrix of (1.12) will be of rank three, 
so that all masses will be expressed uniquely in terms of one of them. 


We note in passing that according to Williams [8] even in permanent configurations of five bodies four 
masses can be expressed uniquely in terms of a fifth, except when four equal masses are located at the vertices 
of a square and the fifth, which is at the center,is arbitrary within a definite limit 


0.439 < i < 1.414. 


This is the essential difference of the present configurations from a permanent three-body configuration, where 
the masses of the system are arbitrary. 


The analysis of the permissible quadrilaterals in the problem investigated by MacMillan and Bartky begins 
with the skew-symmetrical quadrilateral 


Here 
a, Rs— Rs 2G) Ree Re 
mN,= — — uy 5 1 
: Bia eS Me ee 
2 
2 oot ee Va Poe ioe 2 
75. a4 may pee rg = (1 — a) (7; — 73). 
Introducing x = 13/ry and noting that __ j a Ret ; = ” bee these authors obtain the following 
critical values (Fig. 2): 
(a) 1, = Te, ye 1, , m=); 
V 3— 30,42 (1.13) 
V1 —a,+ a? (continued) 
(b) 7s=n7, (ee aa i p10 


2— a4, 
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(c) oy = 0, MN, = 0; 


(d) ify yen7,, topes ay ; (1.13) 
ay ey Ak ae (continued) 
2 
(e) re=N1, os wera Cicer 


7-6 -5 ~~ 


Ae ee Sees f GO 714, 


Fig. 2. Pigs -3. 


Each point of the shaded region gives a real solution of the problem. The intersection point of the two 
subregions where a, > 0 has the coordinates a, =1/2, xy =1/V83 and corresponds to an equilateral triangle 
- with its center at the point my. 4 < 0 corresponds to a convex quadrilateral and 4 > 0 corresponds to a con- 
cave quadrilateral (this means that when cy > 0 the point my is inside the triangle my Mz m3). We note that 


n= MRy can be written as 


n? = 2m,R, + mR, + mR3. (1.14) 


Special cases of skew-symmetrical quadrilaterals are rhombuses and squares. In the first instance my = Mg, 
Ma = my; in the second instance my = Mg = M3 = my. 
The second class of permissible quadrilaterals with axial symmetry consists of equilateral trapezoids (Fig. 3), 


Here 
= peat 
Ee : Vi8/Se 


When an equilateral trapezoid changes from an equilateral triangle to a square the ratio mg/m, changes 
continuously from 0 to 1. For all given values of my and mg there exists one and only one equilateral trapezoid. 
In general all permissible quadrilaterals can be divided into those which are convex and those which are concave. 


For any convex quadrilateral we have the inequalities 


Ty, To, T3) Ta ST0 NBs M6 


with the two extreme sides lying opposite each other. 
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Let us consider Figure 4, where ByOA2, AyOB, and by ay ag be are semicircles of radius rg, drawn with centers 
at the points my, mz and O. These determine the areas A and B. For each point m3c B there exists one and 
only one point mgc= A (and no point outside A) such that the points my, Mg, mg, mg form a permissible convex 


quadrilateral. 


Fig. 4. 


We now give without proof the following important properties of a permissible convex quadrilateral: the 
ratio of the diagonals lies between 1/V3 and v3; each interior angle lies between 60° and 120°; finally, the 
diagonals divide each of the interior angles into two angles each of which is less than 60°. 


In the case of a convex quadrilateral we have 


M1) To) 731013) Ta Mos 


when only ry is a maximal side. In a construction similar to Figure 4 with some limitation of the regions (Fig. 5) 
the following theorem can be proved: 


Fig. 5. 


For each point mg T, + T, there exists one and only one point mg; + Sy such that the points my, mp, 
m3, Mg form a permissible concave quadrilateral. When Mg Sy, we have mgoT}. 


MacMillan and Bartky prove, finally, that for any four masses in a given order there exists at least one 
permissible convex quadrilateral. 


These are the principal results from the theory of permanent configurations in the problem of four bodies 
which will be needed hereinafter. 


2. The Stability of Circular Solutions in the Problem of Four Bodies 


We shall first give the results which i : 
rare ich were obtained by Andoyer and Meyer and which apply to the problem 
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Everywhere in what follows we shall consider p(t) = const, i. e., only circular motions are considered. 
The plane of relative equilibrium of the given bodies will be the xy plane; 
the center of gravity of the system; 


the origin of coordinates will be 
2 axis willben — 


the constant angular velocity of rotation of the X and y axes aroung the 
the angular velocity of circular motion of the points. The equations of motion are 


1 oF 
Cit Sar ra 
. 4 OF 
2nx, = 
Ys + s m, Oy, (2.1) 
Pe Wat) 2 ss 
eer ae 19 J 


Mm, m; 


Fazm Ya+y)+ > 


r 
¢ 
8, J j 


~] RM mts. , ' ; 
In addition, >) m,x, = 0, >y mM, ys = 0, > m,2Z,=On Dy or (x2 + y? + 22) = F + const (Jacobian integral). 
8 s s 8 


With the notation 


pha ets Veer 
Vs == Li — Lys. 
we obtain a : 2 OF 
us + 2niu, = ™m, Ov, 
coli nage ruts Bat OF 
0, == NV; — m, du, 3 
iy, Jeng pO (2.2) 
Zs = m, 02, ’ 
ad Ee 2 ; — v,;) +- (3, — 2;)*. 
4 at Ee ra, = Us — Uj) (Vs vj) + (2s j 
sla aap thers ines I sj ’ 8) ( s # 
8 


erat: ye ig 
= rae : one : en ae =x. ti = const, be. = 
When ug = ag, Vs = bg, Zs = 0 are the positions of relative equilibrium (when ag = x, + iy, » Dg 
= iy? = const in the chosen coordinate system) we have 


OF - .« Ta ‘ eee. et) 
Ba, = na; Dy (% Aj) as; , 
j 
(2.3) 
OF ey 
PF nt by — hm; (bs —b,) ae) = 0, 
j 


sj = [(4 — 4;) (bs — 6;)“8 = = , re = (a — at)? + (y, — y')% 


sj 
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- 
Introducing into (2.2) us=@;-+%s, Us = be+ Ns, 25== 581, where §s» “Nes Cs are assumed 


to be small, and neglecting terms of the second and higher orders, we obtain the following set of linear differ- 
ential equations with constant coefficients: 


. ° . y 3 , 
bs + Qnit, = nbs + > DE Pa 125 Si Py 
j j 


a Py 3 y 1 
Ns — 2niqe =n? Ne + | DEP Dy Psi» 
j I 


(2.4) 
i + Doi Psi = 0, 
4} 
where 
P= Sy Mj; Asj, P53 = — Mj 4sj, 2; P,; = 0, 
j j 
, f Ces ‘ ey 4 
;= M;A5; ———-, Psi = —Mjas3 —>—— » P=) 
Pa 2 j “8) b,— 4, 87 J a b, — 6; > i] ’ 
, bib. ; Dest) i 
Ps = >) Mj Qej ——— nica PF 2 eo > Psi = 0. 
s Ht} : 
j j 
(P’ and P” are conjugate quantities). 
The equilibrium conditions of (2.3) can be rewritten as 
n? ds => > 4; Psi = Dr ae 
: : ; (2.5) 
n2b, = >s b; Py = >| a; P3}: 
j j j 
I. The determination of ¢s. 
Let us set ¢, = C,eht; then 
Cy h? + > CP =) (2.6) 
j 
and the secular equation will be 
Py th Pr. Pr Py, 
A (h?) Pe Poy Poo + h? Po Po 0 
P3y Pp Py3 + h? Py, ial oy) 
Pay Pr Py; Pay + h? 


Multiplying the rows by Vm, Vm, V ms, Vm respectively, and dividing the columns by the same 


quantities, we obtain the determinant of a real symmetric matrix, from which it follows that all roots of (250) 
are real. In expanded form 
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Sy nae tnt a H,h? + H,| =0; 
Hy = my (Ay + Ag + Ay4) + My (gq -E Gag + yy) + my (431 + 52 + 34) + 
+ 1M4 (dq, + gp + O43); 

Hy, = ME (42 G13 + As Ay4 - yg Qy) + M5 (2g G24 + Gog Gay + gy Gog) + 
aie m: (454 Ag + Ay Agq + (gg A34) 4- mt (Q14 Qag + Ago Ay3 + O45 A) + 
My Ms (Ayo 23 + Ay3 Ay4 + Qyq4 Ay + Aye Qy3 + Ay3 Agq + Ayg dy) + 
+ My Ms (413 A3q + Agq A144 + Ay Qy3 + 43 Ay + My2A34 + G34 G3) + 
+ My Mg (Ay4 Ayg + Ogg Ayo + Ayo Ag + 14 G13 + G13 Ayg + Ago A414) + 
+ Mg Mz (493 Ag1 + 431 oq + Aoq Mog + Ggg Dp + Agi Agq + Aq Gy) + 
+ Mz M4 (gq Agy + Ag; eg + Gog Gog + Gog oy + Gy; Agg + O43 Go4) + 

+ Mg Mg (gq Gq + Ag Agq + Ago gq + gq M31 + 51 Age + Ayo Asa); 
H,;=M [4 Ayo Ay3 Ayq + M2 Ag) og Agq + M3 A31 Ago Agq + M? Ay Ago Ogg + 
My My Aye (414. Aq3 + yg Aqq) + My Mz Ay3 (dy4 Ugg + Ay2 Aga) + 
+ my M4 Ay4 (AyQ M43 + 213 Gaz) + Me Mz Ay (Coq Ag + Aoi Ag4) + 
+ Me M4 Ag4 (Gy Ag + Aq A41) + M3 Mg Ag4 (Aq Ag + Gay Ae3))- 


Thus the equation A(h?) = 0 has one zero and three negative roots. When h? = 0, the condition emsCg =0 


can be satisfied only for C, = 0,s =1, 2, 3, 4, i. e., we again obtain Lagrange's solution. Multiplying the columns 
of (2.7) by ay, a2, a3, a4, respectively, and then by by, bg, bs, by, adding and taking (2.5) into account, we see 

that h? = —n? is a root of A(h?) = 0. It is a simple root if all points are on a single straight line (then all ag 

are proportional to b,); otherwise it is a double root. In the first case we can take C, = ag, and in.the second 
case C, =Hag + vb, (H, v = const). 


Thus for any masses and initial perturbations perpendicular to the plane of unperturbed motion, the circular 
solutions are stable in first approximation and the perturbed motion with respect to the z-axis will be periodic. 


II. The determination of &s, 1s. 


Setting & = A,e®t, y,= Be, we have 


1 eS ’ 
Ay (i +n)? +> SAP t+ > DB, Pes = 0. 
: j j (2.8) 
y oo 
j i 


and there is a characteristic equation (see page 69). 


Therefore A(A”) = 0 is aneighthdegree equation in d? with real coefficients. Adding the first and last 
four columns, we find that A? = n” is a double root. Multiplying the first four columns by ay, ag, ag, a4 and the 
second set of four columns by by, bg, bs, bg, adding and taking (2.5) into account, we obtain the following two 


columns: 
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. 
ete ‘tee el ope) cel x9 16 


(sai, eae 3, 4) 


[ ; ) Il 4 z 


“ Miss ones A 2 0. 
Thus, of the eight roots of (2.9) four are known: the triple root \* = n° and the simple root 


Equation (2.9) can also ne obtained in another form. Introducing into (2.1) 
° a ee r, 
r= rt, + 4 Y=Yy, +h y= 2, +E 


, . . . . * . . ns 
x} = const, y; = const, zt = 0 are Lagrange's solutions which satisfy (1.6), we obtain the system in the variatio 


i — In; = n2q; + > Gj Pxit + >) bi9853 
i) j 

be + angi = 2 + Dy oi9 + DOP rir; (2.10) 
j 7 


Xi + DG Pi; = 0; 
j 


——_——— a 


where 


wnt wh =e 
P,yii = >?) M;Qi; = -b 3 a P ii = M;a;; {—3 02 , 
7 "ij raj | 
3 
>) P ii = 9, W= 2, Y; 
j 
0 0) /,,0 0 0 0) /,,0 0 
, (x; — 25) (yj— 93) _  (%;— x) (yj — 99) 
Vii= 3; M;Qi; pale PM its ds a : », ij=— 3m ;0i; ee a : 
j ty "aj 


D4 = 0, 
J 


Considering the equations for the determination of yp, and Vj, We set oj = Aje®" i = Bies. There 
is a characteristic equation (see page 69). . 

Equations (2.9) and (2.11) are two ways of writing the same equation. By setting g =~)? we can derive 
one » from the other. The order of these determinants can easily be reduced by two through use of the double 
root g*=—n?*; but we did not succeed in reducing the order further from the knowledge of two additional 


roots g =—n? and go" =.0: Generally speaking, a determinant of order 2N can be reduced to a determinant of 
order N by means of Schur’s formula: 
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(IT'S) 
2o— zu + hy ca i oft 
wi 7 wll 7 go—u ety 
o= nf 7 af 7 ufl 7 
suz + 5 8vF, 2h 
"hug + cp 
rh &2h Buz + 82h 
Ih eIp Z1h 
(6°3) 
e(u—\) +a = be = we = 
"gy a 2(u—V)+ Fg . wey a 
Bae eq (uv) Hg 5 
wd = aly ns aly 4 
Q= 
bo fies, i ee 
nd = ed = ed a 
OR eae 
a + a; a ae eo 


Pe 


wh 7 

wil 7 

go — gu + mg 
hh 

1eh 

Ih 
8uz + Db 


g(u+y) +8 : 
G 
837 
a J 
C 


Sie; Lee 
an 


a 
=x 


a 


A 
| ca |x oo |N co |x 


= 


a 
oO 


5 


4 
ay 


3 = 
Qy 


$ 
AY 
ain |x =|x =|N ae) 


g(u+y) +%y 


a 


| 
| 
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AB ioe 
a=(¢ Dl=|0n~ cap =|A|| D —CA™1B|=|AD— ACA™B|, 


where A, B, C, D are quadratic matrices and [Al ==; 0°; But in practice this method is ineffectual. 


If the determinant (2.9) or (2.11) is calculated the characteristic equation of the system in the variations is 


A (g?) = g? (g? + n2)3(g8 + ag® + bgt + cg? + d) =0, 
where a = 5n?—(Pyq + Poy + Pag + Pgg), but for the calculation of b, c, d one must know the sum of the principal 
minors of the second, third and fourth orders, respectively, of the determinant A(0). 


It is known that the equation 


agtt + 4a,x3 + 6agz? + 4azr + a, = 0 (2.13) 


then and only then has all real, negative and different roots when the following conditions are fulfilled: 


ay > 0, a, > 0, Gy =O, ig >; ap > 0; 
At — Ap, > 0; 12 (aj — doa)? — apf, > 0; D = fi—-27f3>0, 


2 2 3 : 
where fi, = aa, — 4a,a, + 3a; fo = Apdgd, + 2ayaQa3 — aya3 — aja, — ag. Setting ag =1, 
4a, =a, ba, = 6, 4a; = c, a,=d, we obtain the necessary conditions for the stability of the position of 
relative equilibrium. When D = 0 Equation (2.13) has multiple roots. 


We finally give the formulation of the problem of the stability of the libration points in a limited problem 
of four bodies. |The equations of motion of a point of zero mass will be 


we 2niu = 25, 
b—niv = 29", | 
u 


ae 
Of” 
yom 
F = muv + >, or m, = 0 
j 8 


If ¢= Ae, y = Bet, the characteristic equation will be 


N92 (DFP 1 2\" 9 pp 
\? — 2? (2n3 — + (gP+r)—5Pl == (), (2.14) 


, 


PPL OP Spe aa 


So SENS Stability of Some Particular Circular Solutions of the Problem of Four 


Bodies 


In the case of a straight-line configuration ee (2.9) can be completely investigated. Taking 
ye= =0,s=1, 2, 3, 4, we obtain ag = by = Keo Psj co heat oe sj° If h? is a root of (2.7), then according to (2.8) 


m,C,h? ob D)m,C;Pis =(. 
Jj 


Multiplying the first four equations and then the second four equations of (2.8) by mgC, and adding, we 
have 


[a +n)? — + n?| py m,C.4,— +h? > m.C,B, = 0: 


[@Q — ny — 5 24] SmC.B,— 2 12 Sim.C,4, = 0, 


i. e., (2.9) in this case allows the roots of the equation 


14 — (h2 + Qn) B+ (n? + A?) (n? — 2h?) = 0. (3.1) 


2 


If both roots of (3.1) are to be positive it is necessary that — n? < n< — = n°. But, as was proved by 


' Meyer, all roots of (2.7) satisfy the condition h? < —n? (except, of course, h? = 0, h? = —n?), Therefore (2.9) 
for a straight-line configuration has, in addition to Ay,2 = 0, A34.5 = n, Ag,7,3 = — 2, two pairs of real roots of 
opposite signs:: Ng = — Axo, Ag, = — Age, and four pure imaginary roots: yg = ia, Ay =iB, Ags = ia, Ayg = =iB. 


The existence of the last two roots proves that the straight-line configuration is unstable. Nevertheless, with 
suitable choice of the initial values the straight-line configuration possesses conditional stability, and there 
exist two families of periodic orbits close to the unperturbed positions of relative equilibrium. 


It is also possible to calculate the coefficients of the characteristic equation in some very simple cases of 
plane configurations. The calculations are lengthy and will be omitted. 


I. The square. With the notation 


DT amen HU gle Ie gs d= yg. = Gog Og, = Ug), 
we obtain 
2 23 | o4 4 2 mag? + aa 5 V0) ma? || g4 + n@g? + (3.2) 
Bp) = (gh) |B auees 2 
fe PD 
bere | 
and 
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Since the discriminants of the last two factors are 
= 9 Ne 
Dy = — 16 (4+ V2) ma? <0, Dz = (gy — 34V 2) mPa? <0, 


each of these fourth-degree equations has roots of the forms +- iB, «— ee me ip, —o%-> ip. There- 
fore (3.2) has four roots for g with positive real parts, from which the instability of this configuration follows. As 
in the case of the straight-line configuration, this instability is not of absolute character. 


Il. The rhombus. Considering that 


m, = Msz, My = ™4, n? = 2(m,Q42 + Molo4) = 2 (m4 Q13 + M2042); 


we obtain 
A (g3) = g?(g? +n)? [g* +g? (2n? — M) + 
+ (M2 (—2 + i risria) + Mn? +n} | [gt + mtg? + A] = 0, (3.3) 


where 
A = 36 sin? w cos? © (m, — mz)? + 36 mymz (a{y COS? © + ay4 Sin? w), 


ry © o 
® =arctan ~, 30°<o = 60°, ry. = 1. 
4 


The discriminant of the first biquadratic equation is 


| 


hed 1 212 9 
D, = 9M? (1 ae rors) — 8Mn? = M? [8 (pleceaed Ate (4 aay: )] 
But from (1.11) we derive 
yee R,R; — R2Rg & RoRy.-— RsReo = R5Ro — Ri Rs 
0 Ri + Rs — Rp — Ry Ro+Re—Rs—Re Rs +- Re—Ri— Ro" (3.4) 
Substituting in this case 
{Re ae 


Sees nS 3 v 
2ri3Toa — ("13 + P54) 


we os that Di < 0. Investigating the discriminant of the last factor in similar fashion we obtain D, = nt 4A < 0 
Thus in this case also the characteristic Equation (3.3) for g has four roots with positive real parts. 


We choose the units of length and mass so that 


Pog = T34 = Tan = |, n* = 3(m, 4+ V3 m) dog = 4; 


and let 
3V3 = 
M, = —™, M, = (8V'3 — 1) mg. 
Then 
A (g?) = g?(g? + 1)° [(g4 + pg? + 9)? + r4g?] =0 (3.5) 
or 
A (i2) = 22 (0? — 1)8 [4 — ph? + 72 +g] PA — ph? — rh + gq] =O, (3.6) 
where 


i 3 9 
p=1—+-M,, r?=3M,, g=-}(1+M,)—->-M, (1+ 3M). 


Let us consider the equation 


f(a) = (2? + pe +qy+rz=0. 


y 
Yj (2) 


Ye) 
Fig. 6. 


= (x? + px + q)? to the straight line yp(x) =— ix (Fig. 6), we 


Analyzing the relation of the curve y4(x) 
for the reality (and consequently the negative character) of the 


find that the necessary and sufficient conditions 
roots of f(x) = 0 are 


A) p>0; B) q>9; jy ea pl 


This means 
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2 ape ee 
A) Maxi B) Mi<aapamy ° 


Seen 
3M, V 2—3M,— ~,- M3 +6M2+4 
C)M,> 9(1 +3 M2) 


But conditions B) and C) lead to 
27M3 + 522M3 + 36M, —8 <0, 


SY=5. «9 55015, 


from which Mz < 0.093860. At the same time it follows from B) that Mg < 


Thus for any finite masses this configuration is unstable. 


IV. Isosceles trapezoid. In this case my = Mg, Mg. = My, nee (mq + M3) (ayo + ay3) (Fig. 3). The 
characteristic equation could not be investigated; we note only that in the limiting cases where the trapezoid 
degenerates into an equilateral triangle (mg = 0) or a square (Img = mj) the configuration is unstable; in the 
first instance we have 


27 
A (g?)=9%(e4+1)° (6? +2¢°-+7) = 0; (rin = 1, m = 1), (3.7) 


and the second case is discussed above. 


V._Rhombohedron, The characteristic equation for this configuration could also not be investigated. 
We only note some critical cases corresponding to a), b) and c) of (1.13). 


Let ty2 = Tyg, T23 = 134, Mg = My (Fig. 1) and a point is found on the straight line mg, mg. Then with mg, = 0 
the configuration is unstable, since otherwise the collinear configuration of the three points my, my, my would 
also be stable. This follows from the analysis of the characteristic equation: 


A (g?) = g7(g? + n*)® {[g? — n? — 2my,a,3 + 


+ 2meA23 (1 — 3 cos? B)] [g? — n? + myay3 + 2mza23 (1 — 3 sin? B)] + (3.8) 
3 a 47 
+ 4n?g?} [e + g? (m, ay ms) a2 —Z (3m, ts ms) mac, | = 0, 
Tie 


B =arctan -, 
18 


the last factor of which has two pure imaginary roots for & and two roots which are equal in absolute magnitude 
but opposite in sign. The second biquadratic equation, which according to (2.14) determines the stability of the 
liberation point of ms, was not investigated, because, according to Huttenhain [9], the libration point in this 


case is unstable if my/mz > 11.72 138 (libration points lying on the axis which passes through m2 and my are 
unstable), 


Let us consider the equilateral triangle my, Mg, mg where mg = mg, (Fig57); 


The coordinate origin is taken at my. The fourth mass my, = 0 can be either outside or inside of this tri- 
angle. 


Calculating (2.9), we have 
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A(X?) = 22 (02 — ni)? fs + 22 (Py, — 2n?) + 


+ (G = & Paul ae - Pu Pil} [2 — nn? + # (mys =} oy ‘at, | = os 


The problem of the stability of such configurations also divides 
into two problems: the stability of an equidistant configuration of 
three points where two of the three masses are equal, and the stability 
of the corresponding libration point in the limited problem of four 
bodies. The first follows from the fact that the positive discriminant 
of the last factor leads to the inequality 


m? — d0mym, — 23m? > () (3.10) 


which also follows from (0.1)*when mg = mg. 


We denote the libration points by C, A and B. Each of the three 
systems of four mass points is stable if and only if, according to (3.10), 
my/mg > 50.4558, and the corresponding libration point is stable. 
The libration point on the axis x = 0 is determined by the following 
equation given by Meyer: 


EES 
Vora: u ee (y a z) = (3.11) 
Yate A — B) — a seny (y2+V3 y+ 4)? == 0, 


where 739 = 1, me = m3 = ——, M1 — tL. 


Meyer also showed that there always exists a libration poing C between 0 and i; the point B is between 
— V3 and—1, but when wp > 0.43242 the libration point A ceases to exist. Hence it is clear that the configu- 
ration mm ,Am," is unstable, because in our case 1 > 0.96187. 


It follows from (2.14) that for the stability of the libration point it is necessary that 


2n? — P>O, (P + 2n?)? —9P?2 >0, D=9P? — 8r?P > 0. (3.12) 
In our case 
a = AEE Mee eared ee a 
n® = (m, + 2m,) aig = mee [yt 4 V3yiiye’ 


Te 4 4 Mom) 
ee ae we mage BY. 


Giving a value to y, it is easy to find the corresponding values of # and then the values of P and P’. Cal- 
culations showed that the configuration my Mm2Bmz, is stable for my/m, 2 k = 367.0 (with y> —1.00153), and 
that the configuration my mz2Cms is unstable for any finite masses. 


* As in original — Publisher's note. 
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4, Conclusion 


We now make an attempt to apply the theory of permanent configurations of four bodies practically to 
the Trapezium of Orion, which consists of six main stars, two of them, E and F, differing considerably from the 
others A, B, C, D, by their large proper motions and radial velocities. 


The projection of the Trapezium of Orion on the celestial sphere looks approx- 
et imately like Figure 8 [10]. 


Observational data communicated by Professor P. P. Parenago is given in the 
, : table. The last measurements were made in 1920-1930 and the first measurements 
\ in 1840-1850; @ is the position angle, rj; is a distance between stars in the plane 

! of the figure; JI is the mass of a star in terms of the solar mass. 


\ 
| \ The distinct difference in the proper motions of stars E and F, on the one hand, 
ib ee PE and A, B, C, D on the other hand, suggests that E and F are only projected on the 

jet 4 general background of the trapezoid. But even if this is not the case, by using the 

smallness of the masses of E and F and their closeness to A and C, respectively, the 

Trapezium of Orion can in first approximation be represented by a system of four 

bodies A, B,C and D, where the distances CD and AC.are approximately equal and 
the masses D and A are also approximately equal. If this observed system formed a permanent configuration it 
would necessarily be a rhombohedron (or a quadrilateral very similar to it), which is not the case in reality. 
We now inquire whether the quadrilateral ABCD is only the projection of some rhombohedron lying in a plane 
that forms the angle i. with the plane of the figure. Let us now consider the following problem of descriptive 
geometry: to construct a rhombohedron from its projection. 


BG oo. 


AB 


8"770| 127913 | 24305 3884 167750 | 197082 | 137334 37743 


..,| 87748] 12"942 | 217477 | 3”927 | 167723 | 197271 | 13°489 | 37913 | A= 40 

if} s"701} 127944 | 217595 | 47054 | 167598 | 197333 | 137432 | 3”894, | B=6 
87730 137004 | 247550 | 47204 | 167759 | 197318 | 137332 | 47005 | C =24 
32°30] 131°36 | 95°55 | 352°28 | 163°29 | 299°59 | 240°58 | 122°40 

9 | 32°28] 131°28 | 95°57 | 352°02 | 163°66 | 299,43 | 240°97 | 122°70 | B=!2 
32°34] 131°55 | 95°73 | 354°38 | 163°91 | 299°69 | 241943 | 4122°59 | B=2.5 
32°20] 131°20 | 95°58 | 350°90 | 163°97 | 300°00 | 241°36 | 421°¢0 | “=2> 


Fig. 9. Fig. 10. 
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Let the projection of the segment a in plane P on plane Q be c (Fig. 9). In a parallel projection the seg- 
ments and the angles between them are related by the formulas: * 


c=aV1— sin?esin?i, a=cVY1+sin?Ptg?i, tg¥= tg ycosi’, ) 


sin B ie, sin a cosi 
V1 — sin? 9 sin? i V1 — sin? (p | a) sin? i : (4.1) 
sina = Sth eevee Care 27 ee | 
V1 + sin? ¥ tg?i V1 + sin? (¥ + B) tg?i 


When in plane Q a quadrilateral is given with the sides cy, C2, C3, C4 and angles By, Bo, B3, By, (Fig. 10), 
the position angles of its sides are related by 


Yn = 180° + V1 —B, = ¥, + (n —1) 190° — S8,, 


ets (4.2) 


which is valid for any (convex or concave) n-sided figure. 


Thus from a quadrilateral with the sides cy, C2, C3, C4 (Cy # Co, Cz # C4) We require the construction of a 
thombohedron with the sides ay = ag, a3 = ag, lying in a plane inclined at the angle i to the projection plane. 
The equalities ay = a2, a3 = ag lead to the equations 


peo es : 
‘tg? i (cj sin? VY, — c} sin? V,)= c? — c?, 


hee (4.3) 
tg? i (c3 sin? Y, — c{sin? W,) = cz —c?, 


The angle Wy is determined from 


Atg?¥, + Btg¥,+C=0; (4.4) 


A =e} (ci — ¢) — ¢3 (ca — 3) COS*By — 5 (Cg — C1) cos? (By + Ba) + 


“cy (Ch ¢;) COB’ - 
B = c2 (c2 — c2) sin 28, + cg (cz — cy) sin 2 (8, + Bs) + cf (cp — cj) sin 28,; 


C= — A (h—o2) sin, — cf (cf — ef) sin® (Pa + Ps) + (cd — ef) sin? 


A+C=0. 


The solutions of (4.4) are 


* In this paper ‘tg’ = ‘tan’, ‘ctg’ = ‘cot’ — Publisher. 


TA 


/ —B+LV B+ 4A} * .B—V B+ 4A 
Tg ee tg Y= ———_,,——__ 


~~ 2A 


from which 


tg V7, + ctg i = 0, oe oo" 


Thus (4.4) has four roots ¥"'y, ¥"4, Vi; ¥{3 separated from each other by 90°. Therefore 


, ” . 
sin? Yh j= (oe? iy a farky ane ye 


and from (4.3) we obtain 


ctg? i’ + ctg?i” = — 4. (4.5) 


ctg”i" = 0 gives the formal solution of the problem: 


(as and ag are diagonals). 
Therefore the solutions of (4.5) are 


ctg?i" = 0, ctg?i’ = — 4. (4.6) 


When i = i’ we have 
a,=c,|cos¥al, n=1,...,4, a,=0, a,=—a, + aj. 


Consequently in this case we obtain a straight line as the limit of a rhombohedron with ay = ag, ag = ag 
and the diagonal ag—* 0. Any quadrilateral (cy # Co, C3 # cg) can be regarded as the projection of some 
segment, with the angle between the respective planes i = i‘ determined from (4.6) and the position angle y = 
= 1 from (4.4). 

The solutions (4.6), which are not rigorously derived from the theory, were checked a few times using 
numerical examples. One of these examples is the Trapezium of Orion. We first note that if.we begin to con- 
struct the quadrilateral ABCD directly from observational data (using the sides rij) we will not obtain a closed 
figure because of errors. Therefore, taking definite values for CD and BC and the angles A, B, C, we determine 
all remaining elements of the quadrilateral. The final values are: 


AC = c, = 12.886; | A=B8,= 99°00’00"; 
CD = c, = 13.332; | C =B, = 110°09'36"; 
DB=cy = 19.237; |. D= 8, = 58°38'24”; 
BA= c, = 8.8878: |B =f, = -92°42'00’. 
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From these values weobtain bi = 54°34" 59", » } = —35° 36" 01" and the values of (4.6) subject to the 
accuracy of the calculations. We add that the straight-line configuration of three points which is obtained for 


i =i': 


a,=a,=7.501, a,;=a,=—7.947; a,=0, a,=15.448 


is not permanent for the given masses. 


Thus within the framework of the present model (with the stars A and D of identical mass) the Trapezium 
of Orion cannotform a plane permanent configuration. It is only possible that the stars of the Trapezium form 
a spatial permanent configuration in the form of a reguiar tetrahedron, or that, according to the Pohlke- Schwarz 
theorem, any quadrilateral in a plane can be regarded as a parallel projection of any given tetrahedron. 


In conclusion I wish to thank Professor G. N. Duboshin for directing this work. 


P. K. Shternberg State Astronomical Institute Received April 18, 1956 
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ON SOME PECULIARITIES IN THE MOTIONS OF COMETS 
P. G. Dukhnovskii 


The changes in cometary orbits, caused by planetary perturbations, are 
considered. This problem is solved approximately by the doubly averaged 
limited circular problem of three points. 


A formula for calculating the period of rotation of the line of apsides 
and formulas for the amplitude of the variation of inclination have been 
derived. 


Calculations according to these formulas, in the case of perturbations 
from Jupiter, for cometary orbits with semi-major axes within the limits 
10-100 AU and perihelion distances from 6 to 10 AU, show that the ampli- 
tudes of variations of inclinations are small and reach only 2-3° and that the 
periods of rotation of the line of apsides are measured by tens of thousands, 
hundreds of thousands and millions of years. The variation of the inclina- 
tion angle can have no special cosmological significance. The rotation 
of the line of apsides causes a periodic variation of the smallest distances 
to which comets can approach a planetary orbit. 


These distances change within the limits q~ R to p— R, where q, p 
are the perihelion distance and the parameter of the cometary orbit, R is 
the radius of the circular orbit of the planet. As a result of the variations 
of these distances, comets can evidently sometimes approach close to 
planets and be transformed into short-period comets. 


1. Introduction 


Several investigators have established [1-3] that the solar system includes a huge number of comets moving 
at great distances from the sun. Such comets are inaccessible to observation and can exist practically indefinitely 
since their great distance from the sun enables them to retain their constituent gas. 


A comet can be observed only when the perihelion distance of its orbit is reduced so that it falls into the 
region of visibility. 


The orbits of comets in the zone outside of Neptune are continuously changing due to perturbations by 
the planets. This is evidently the reason for the appearance of new short=period comets. 


In the present paper we shall endeavor to explain how cometary orbits change over a long period of time 
due to planetary perturbations and how this leads to the appearance of observable short- period comets. 


PAS Planetary Perturbations 


We obtain an approximate solution for the changes of cometary orbits under planetary perturbations by 
using the doubly averaged limit circular problem of three points. 
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The differential equations of this problem have the following integrals [4]: 


a = const, 


v (1) 
V Pcosi = const. 


(2) 


For the purpose of determining how the angle of inclination _i changes along with the orbital parameter P, let 
us consider the following pair of differential equations [5]: 


2m 


QT 
di ; 
dt aE | | rWeost ane dM’, (3) 
do 4 se 
di EVE | \t— P cosVS + (r+ P)sin VT —re cot isinUW]dM dM’. (4) 


Here r, P, e, w, i, V, M, U denote, respectively, the radius vector, parameter, eccentricity, angular distance of 


the perihelion from the node, inclination, true anomaly, mean anomaly and cometary latitude; M' denotes 


the mean anomaly of the planet; S, T, W are the components of the perturbing acceleration along the radius 


vector of the comet, in a direction perpendicular to the radius vector and lying in the plane of the osculating 
orbit, and along the normal to the orbital plane. 


We represent S, T, and W as follows: 


S = 5S, —S,; (5) 
Pie) ae (6) 
W =W,—W,, (7) 
soils CR | 
‘ee (= Me ip (Fe ie t (8) 
mia(E ys = CHE | 


Here A is the distance between the comet and planet, m is the mass of the perturbing planet and R is the 
radius of the circular planetary orbit in the plane of the ecliptic. 


From Figure 1 it follows that 


k2m (r — R cos L) , 
Spe asilla ite SR 


e : (9) 
k2mR ; ig 9 10 
T, = ——jr (sinU cos M —cosU sin M’ cos i); (10) 
2 
Wi= as sin M’ sint. Coy 


Here L is the angle between R andr. 


From simple geometrical considerations we find that 
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es | Rsinn'sint 


4 
s 
SS 


Fig. 1. 
ae 
S,= St (cos U cos M’ + sin U sin M’ cos i); (12) 
k?m |. P : , “\. 
T= — a (sin U cos M’ — cos U sin M’ cos i); (13) 
W,= — sin M’ sini. (14) 


Substituting in Equations (5), (6) and (7) instead of Sy, Ty, Wy and Sg, T2, W, their values from (9)-(14), we 
obtain 


S sx eed = (cosU cos M’ + sinU sin M cos i); (15) 
T = km R(s — gz) (sin U cos M’ — cos U sin M’ cos i); (16) 
W= Mf Ae sin M’ sin i} oe ~~ sin M’ sini. (17) 


Substituting these values of S, T, W into the differential equations (3) and (4), we have 


2m 
\(e ee sin M’ sini 
0 


ai en Bp + | cos U a’ dM; (18) 


dt 4n2k VP 


on 
dw k2m \ 


e =. = 
dt 4rn2k VP ) 


{— Poos v [745° * -— Fa (cos U cos M’ + 


ony 2. 


-b sin U sin M’ cos i) | + (r + P)sinV |— (sin U cos M’ — 


— cos U sin M’ cos i) + sa (sin U cos M’— cos U sin M’ cos i)|— 
—recot i — sin U | 5 sin M’ sini + sin M' sini}} dM dM’. (19) 


Considering that through integration the second term in (18), and the second, fourth and sixth terms in 
(19) vanish, and also transforming by means of the formula 
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2 
dM ay PB dV (20) 
from dM to dV, we obtain 
di____kmRsini tt 78 
bees km Rsini ? ( r3cosU sin M’ i 
a= Tae \ Vane. dV dM’: (21) 
00 
QT QT 
do ht km \ ’ 2 { ~~ PcosV (r — Roos L) 
es 4 AS ai 
00 
__ (r+ P) sin VR (sin U oe cos U sin M’ cost) _ Rerctgi ee U sini | dV dM’. (22) 
It is easily seen that A? = R? + r?— ORrcosL. We now find 1/3 
JA ae 4 ae 4 [1 oz 2Rr cos r= 
A’ (R? + r? — 2Rr cos 1)? (R® + r2)"!* Rear? 
1 3Rr cos L 1 3Rr cos L 
=———, [1 + 2 aa | = ose Ge ee (23) 
De ava CR a en ay ale | 


The last terms in the expansion (23) are small enough to be neglected 


Substituting into (21) and (22) instead of 1/A* its value from (23), we obtain 


2 2 
pee eer ne ia \ \. 3cosU sinM’ : ane) ee dM dM: (24) 
dt 4r2Pa l2 (R2 ) (R2 +. r2)"l2 
09 
2m 27 
g2? bs Ve 4 EHD | 5 
dt 4r2Pa | , (R? + r2y'2 (R? + r2)'l2 

0 


x [— PcosV (r— Reos L) — (r + P) sin VR (sinU cos M’ — 


—cosU sin M’ cos i) — RerctgisinU sini]dV-dM’. 


(25) 
2 
4 4 : R 
i i Se eee st and ee POW EISIOL : 
In order to simplify (24) and (25) we expand ona aoe p ay. 
1 “| R2 \“*Is 4 [ Se Kay Site | 
———~_ => —) =—|i-=—; —— ae ) 
(R? + r2)'2 inane (1 Be =) r8 : a 7 ai 8 r4 ste (26) 
4 3 R? 15 Ra ; 
eee eee ei sogmay: uel 4 
4 4 FY "=5[1-3 o FR 30 R3 «|= 
(R24 raya or (1 er een, ae es (27) 
1 oy eh eee 
gs es es a er Ea 
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Since we are considering only those comets for which we always have r > R, we can confine ourselves to 
only the first terms in these series, 


It is easily seen that the following simple relations exist: 


cos L = cosU cos M’ + sin U sin M' cosi; (28) 
cos U = cos V cosw — sin V sina; (29) 
sinU = sinV cosw—+cosV sinw. (30) 


We substitute in (24) and (25) instead of 


4 4 


See Sidi — ee fs COS 17, COD aie, 
(R? + r2)'ls (R? + r?) la 


their values from (26)-(30). After integration and a simple transformation we obtain these equations in the form 


“di 45kmR¢ “A Tr 

= ae (a — P) sin 2: sin 2u; (31) 
dw 3 kmR? a: Wagt tg 
a= B Shape (3-5 cos? i — sin? i cos? w — 4.5]. (32) 


For more convenient integration of these equations we transform them through introduction of the quanti- 
tiesa, P, r, A, f and W by means of the following relations: 


R, 
a=, (33) 
P = aP a Op. 
a (34) 
| 
dt ae iL 7 dt, 
-kma (35) 
from which 
R’l 
t—t= 
kma'ls ” 
V Pcosi= A, (36) 
from which 
— A2 
Lee cos? i ’ 
(37) 
cost =f; (38) 
20 = W. (39) 
Using these formulas, we transform (31) and (32) to 
af 45 a2 : 
a = 3 gl (PAY) (1—P) sin W; (40) 
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aw 
Ge = REL +P) — 5 (1 —P) 008 W. oe 


We shall seek the solution of this set of equations as series of powers of the small quantity a: 


f= fot af, 4 otf, 4 (42) 


W=W,+ a°W, + aW, ie (43) 


Differentiating (42) and (43) and substituting the values of these derivatives into the left-hand sides of 
(40) and (41), and also substituting into the right-hand sides of these equations instead of £ and W their values 
from (42) and (43) and equating coefficients of equal powers of a, we obtain 


d 
we a= Q, fe; “fo = const; (44) 
dw 3 f 
=F a [4(1 + 72) — 5 (1 — f2) cos Wo). 
(45) 
Transforming (45), we obtain 
dW, An — fF 
ae = a + f)[1— — cos Wo]. (46) 
We introduce the notation 
qo 0 i+f)=B 
Aa 0) = (47) 
and 
—ft 
8 3° (48) 
: 4 + fi ; 
Substituting B and n into (46), we have 
aw, 
dt 2 == B (1 — ncos W,) (49) 
For the determination of fz we obtain 
df 45 fF 0 
dt Setoea “A8 (fo ere A?) (14 — f2) sin W,. (59) 
To determine f, we perform the following transformations. 
We can write 
dfs _ dfs aWe 
dt dW, at ° (51) 
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Using (34), (37), (47), (48), (49), (50), Equation (51) can be represented as 


sin Wo 
pepe) SL 0 ees 


dfy Ldn foa 
le P » (4 — P) Tacos W ; (52) 
By integrating we obtain 
h= Pah (e*) In (4 — ncos W) + const. (53) 
Therefore, neglecting terms containing higher powers of a, f can be represented in the form 
>for ® safe Pin (4 — ncosW,). (54) 


This equation shows directly that f is a periodic function of Wo. If in the expansion (43) we neglect 
terms containing the small quantity a? as a factor, we shall have the equality 


Wo = 2w; 
from which it follows that f is a periodic function of w. Equation (54) shows that f has a minimum value at 


Wo = 0, that is, at w = 0 and a maximum at Wo =17, that is, at w = 1/2. 


From these relations we obtain the formula for calculating the amplitude (am) of f in the form 


(55) 


In calculating amf from this formula we shall assume that for t = 0, Wo = 1/2 and thus f = fg, that is 
fo is the initial value of f. 


Equation (55) was used to calculate amf from orbits with different semi-major axes, different parameters 
and different inclinations in the case of perturbations by Jupiter. These calculations showed that the amplitudes 
of f are so small that the corresponding changes of inclination i amount to only 2° and in a few instances 3°. 

It follows from (2) that with a change of i there will also be a change of the orbital parameter P, and that a 
change of the parameter is accompanied by a change of the perihelion distance of the cometary orbit. This 
change is small and does not exceed 0.2 AU. Such changes are so insignificant that they cannot play an impor- 
tant part in the transformation of cometary orbits. 


We shall now consider the result of the rotation of the line of apsides of a cometary orbit. We shall begin 
by deriving a formula for calculating the period of this rotation. 


Since f changes very little with time, we shall assume in (41) that f is constant and equal to fy. 


Then, using (47) and (48), we represent (41) as 


dw 
Bdt =: erates (56) 
from which 
Ww : 
Ba Ne ee (57) 
{ = ncos 4 
i) 
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For convenience in integrating we introduce the new variable o which is associated with changes of W 
by the relation 


1—n oe (58) 


It is easily seen that we then have the following relations: 


dha yg ’ (59) 
1-+neoso ’ 


dW = aieer do. 
1+ ncoso (60) 


1—ncosW = 


Substituting in (57) instead of 1— ncosW and dW their values from (59) and (60), we have 


1 
Br = ma dc. (61) 
0 


from which 


Vin (62) 


We shall assume that when o is increased by 27, t changes by T', and shall therefore write 


o-+ 2 


B a5 c)= SSS SSG (63) 
( . ) Vin 
From (63) we now subtract (62) and obtain 

pee Ee (64) 

BV1— n2 
We now go over from T to t. 
When T varies by the amount rT’, let t vary by the amount t’. Then, using (36), we obtain 
3} 
tees ee (65) 


kma !2 


Using (64) and (47), we obtain a formula for calculating the time interval during which w changes by 7, 
in the form 
Qra’lt Aa 


ue 3ftR2mk(1 + f2)V1—n® (66) 


Taking the solar mass as unity, the astronomical unit as unit of length and setting k equal to 27, t' is expressed 


in years. 
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By using this result as well as some of Equations (33)-(39), we can represent (66) in the form 


v= le Ne aay . (67) 
3mR? (4 4- f2) VI — 


This formula has been used to calculate the periods of variation of the angularseparation of the perihelion from 
the node which results from perturbations by Jupiter. The calculations were performed for orbits with semi- 
major axes from 10 to 100 AU at intervals of 10 AU, inclinations from 5 to 30° at intervals of 5° and fot peti- 
helion distances from 6 to 10 AU at intervals of 1 AU. 


The results are shown in Tables 1-6. 


The tables show the periods of variation of the angular distances of the perihelion from the node for 
cometary orbits with different semi-major axes, inclinations and perihelion distances. 


In order to obtain the period of variation of the angular distances of the perihelion from the node in years, 
each value in Tables 1-6 must be multiplied by 10*. These tables have been used for the plotting of curves 
showing the dependence of the periods of change of w on the semi-major axis for different angles of inclination 
and different perihelion distances. 


The curves for the dependence of the period of variation of w on a for different i and.q are shown in 
Figures 2-7. 


TABLE 1-2 
For t= 95° For i= 10° 
q q 
6 7 8 9 10 6 7 8 9 10 
a a 
10 3 3 3 4 4 — 
20 12 12 16 19 23 26 
30 20 20 30 44 50 60 
40 40 41 94 69 86 104 
50 58 09 78 100 124 AS 
60 78 79 106 136 168 204 
70 100 104 136 174 | 217 265 
80 123 124 168 216 | 270 | 329 
90 148 150 | 202 264 327 399 
100 175 177 240 | 309 | 387 474 
TABLE 3-4 
For i= 15° For i= 20° 
q q 
6 7 8 9 10 6 U 8 ) 10 
a a ‘ 
40 3 3 4 4 — 10 3 4 4 4 — 
20 13 16 20 23 27 20 43 16 20 24 28 
30 26 34 42 51 61 30 26 32 43 03 63 
40 42 50 70 87 103 40 43 o7 72 90 106 
50 60 80 | 102 126 153 30 62 82 105 130 157 
60 84 108 138 172 | 208 60 83 444 142 EL 214 
70 103 138 178 227 270 70 106 142 183 | 228 277 
80 127 174 224 PH ta Ya ree 5) 80 1314 176 PRAY 283 345 
90 154 | 206 | 266 333 | 406 90 158 242 274 342 418 
100 181 244 | 314 | 394 | 483 10 186 206 324 405 495 
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TABLE 5-6 


. For i= 25° For i = 30° 

q q 

6 7 8 9 40 6 u 8 9 10 
a a 

10 3 4 4 4 = 10 3 4 4 4 — 
20 43 17 24 20 29 20 14 18 22, 26 30 
30 27 33 45 54 65 30 28 35 46 57 67 
40 44 58 74 92 110 40 46 61 77 96 114 
50 64 850 |) 108 91 13457) 162 5 66 89 | 113 | 140 | 169 
60 Son Ada 147) =| 482 ie 4 60 89 | 449 |; 153 | 190 | -230 
70 109 | 4147 189 | 235 | 286 70 AA SDS wel ehOOm |e OS 
80 BOS) 182 234 292 356 80 140 189 244 304 371 
90 NCSe 219M 283.1 3353 a 434 90 169.) 5 2262 | 9294 4) 36317429 
400 192 |, 258 | 3384 | 448 (1/9544 100 199 | 269 | 348 | 436 | 532 
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Fig. 2. Fig. 3. 


The data in the tables and the curves show that the angular distance of the perihelion from the node of 
cometary orbits varies from 0 to 27 in tens of thousands, hundreds of thousands and even millions of years. 


It is easily seen that the rotation of the lines of apsides of cometary orbits results in a change of the mini- 
mum distance of approach of the comets to planetary orbits. These distances will be smallest in the two cases 
w = 0° and w = 180°. In these two cases the smallest distance of approach by a comet to a planetary orbit is 
q—R, where q is the perihelion distance of the cometary orbit and R is the radius of the circular planetary 
orbit. 


This distance will be greatest when w = 90° or w = 270°. It will be greatest for comets moving in orbits 
inclined at 90°. In this case the distance is P ~R, where P is the parameter of the cometary orbit and R has 
its previous value. 
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The changes of these distances can evidently result in close approach of comets and planets; this can 
result in a radical transformation of cometary orbits and the appearance of short-period comets. 


In conclusion I wish to thank A. A. Orlov for suggestions. 
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THE DISTRIBUTION OF COMETS OF JUPITER'S FAMILY 
K. A. Shteins 


A two dimensional function of the distribution of semi-major axes and 
parameters of the orbits of comets, captured by Jupiter, is determined by the 
generally accepted method. The distribution of long-period comets, disinte- 
gration and the probability of discovery of short-period comets are taken into 
account. It is assumed that the probability of discovery of short and long- 
period comets is the same and determined from observations empirically. 
According to observational data it is assumed that the distribution of the peri- 
helia of long-period comets is proportional to the distance of the perihelia 
from the sun. The influence of the disintegration of comets is investigated 
in accordance with suggestions of S. Orlov and B, Levin and the gravitational 
hypothesis. The results are in good agreement with observations. The almost 
complete absence of short-period comets with retrograde motions is a result 
of the fact that these comets, in most cases, have smal] periheliondistances 
and disintegrate rapidly or, to a lesser extent, that they are transformed 
from a comparatively small number of long-period comets, which have small 
perihelion distances. 


The Plane Problem 


§1. The existence of a group of short-period comets whose aphelia are close to Jupiter's orbit, shows 
that Jupiter has played an important part in their origin and perhaps intheorigin of other short-period comets. 
Investigations by Tisserand, Schulhoff and others [1] have shown that some of the characteristics of Jupiter's 
comet family indicate that this group of comets originated through capture by Jupiter of some long-period 
comets. H. A. Newton [1] investigated the theoretical distribution of short-period comets and showed that 
direct motion must predominate over retrograde motion among short-period comets. Newton's results cannot 
be compared exactly with observations because he considered one-dimensional distributions. The majority of 
scientists have regarded Newton's investigations as confirming the hypothesis that short-period comets originated 
through capture; but others, like S. K. Vsekhsviatskii [2], have regarded these same investigations as disproving 
this possibility. Newton's calculations show that comets whose planes have small inclinations should arise with 
direct and retrograde motion in about the same numbers. In actuality, the majority of short-period comets 
have small inclinations and are not observed to have retrograde motion. B. M. Shchigolev [3] has made a special 
study of the orbit elements of periodic comets captured by Jupiter. He considered the plane problem for the 
hypotheses of Newton and of Fesenkov. According to Newton's hypothesis the velocities of long-period (parabolic) 
comets are distributed uniformly at each point of interplanetary space, and according to Fesenkov's hypothesis 
the velocities of parabolic comets outside the sun's sphere of influence can have completely arbitrary orienta- 
tions. Fesenkov’s and Newton's hypotheses were used for ease of calculation and are, first of all, the source of 
the radical discrepancy between Shchigolev's results, which gave approximately identical probability for direct 
and retrograde motions, and our results, according to which direct motions are almost the only possible motions. 
Secondly, Shchigolev used Newton's approximate formulas, thus introducing another small difference in the 
results. Thirdly, Shchigolev did not treat the matter of comet disintegration; in 1924, when his article was 
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written, this question had not arisen. It must be emphasized that Shchigolev’s work shows that an investigation 


of direct versus retrograde motions requires great caution. His investigation shows that Fesenkov’s and Newton's 
hypotheses lead to different results. ‘ 


In the present article we again investigate the distribution of the comets in Jupiter's family, but we take 
into account the distribution of parabolic cometary orbits, the lifetimes of short-period comets and the probability 
of discovering short-period comets. We assume that we have a plane problem because the inclination of the orbits 
of Jupiter's family is small. This is not entirely legitimate. Although all the members of the family of short- 
period comets possess small inclinations, in the original family of long-period comets arbitrary inclinations are 
encountered. In the second part of the present article we shall consider the three-dimensional problem, and 
shall determine the order of the error introduced by our assumption. Since the three-dimensional problem leads 
to fourfold integrals, the modest calculational resources of the author will cause it to be treated less precisely 
than the plane problem. We have therefore considered it necessary to divide the work and to consider the plane 


problem separately. We propose subsequently to examine the problem of double capture of comets and similar 
questions. 


In its general features the scheme of this investigation is classical, We assume that a stationary flock of 
parabolic comets is revolving around the sun [4]. Stellar perturbations [5] insure the stationary condition. Per- 
turbation by Jupiter changes some of the flock of parabolic comets into short-period comets. We assume that 
Jupiter has a circular orbit of radius Rp around the sun. Jupiter does not influence the motion of a comet at a 
large distance, but inside its circle of influence only Jupiter attracts the comet. Dynamic effects of the rotation 
of the circle of influence are not taken into account. As a result of the assumed symmetry the passage of a 
parabolic comet through Jupiter's circle of influence is considered for only one definite position of Jupiter. 


We shall describe the path of a comet which undergoes capture. The orbit of the parabolic comet is 
determined by the perihelion distance q and perihelial longitude w. The longitude of the perihelion is computed 
from the straight line CJ which connects the sun (C) and Jupiter (J). The position of the comet in its orbit is 
not fixed, since it is assumed that a continuous flow of comets takes place along the orbit in both directions. 
The comet intersects the circle of influence at the point K(a, P), where @ is the angle in Jupiter's polar coor- 
dinates with the pole J and polar axis CJ. The radius of Jupiter's circle of influence is denoted by p. The angle 

“between the comet's velocity vj, and CJ is denoted by B. The comet enters the circle of influence with the 
relative velocity vy, which is determined by the coordinates a,6y. To determine the motion of the comet 
inside the circle of influence we introduce the coordinates (&, By), where By = By — @ is the angle between the 
normal of the circle of influence at the point « and the relative velocity. Inside the circle of influence the 
comet describes a hyperbola and leaves the circle at the point @ + 2g. The relative velocity of the comet 
v';(vr = v;) at this point forms the angle 1 —8, with the normal. The absolute velocity va = v'; + vy, where 
vy is the velocity of Jupiter, forms the angle Bg with CJ. The elliptic orbit of the captured comet is determined 
by the semi-major axis a and the parameter p. 


The problem of the present article is the determination of the theoretical distribution f(a, p) of captured 
comets. The problem will be solved in the following form: 


\ 14 (a, p) | dadp = 


D “ D a, 6 ra D a -+- Pa, Po 
Be RAN cosy Dee ce pw plaad wan 


» P) 


Here f,(q, w) is the distribution function of parabolic comets, f(a, p) is the probability of the lifetime 
of short-period comets and f,(a, p) is the distribution function of the discovery of comets. The Jacobian 


D (a+ Bs, Be)| is defined as the area between the curves a =cj, p=c'y in the rectangular coordinates 


D(a, p) 
a + By, By. This procedure is simplest and the limits of integration are obtained directly. 
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§2. From ACKJ we have 


(CK)? = R5 + p? + 2Kop cos a, (2) 
sin < KCJ = 


elle, 
(GK) sino. (3) 


We shall regard p/Ry * 0.06 as a first-order small quantity and shall hereinafter restrict ourselves to only 
first-order quantities. Therefore 


ify ° 
Ch = ts E + h, cos «| . (4) 


= KGl ee i sin a. (5) 


The parabola intersects Jupiter's circle of influence at two points. We shall consider each parabola twice, 
because comets move along parabolas in both the direct and retrograde directions. Then one definite pair of 
values of a , B corresponds to each parabola. Directions on the parabola are chosen so that in relative motion 
the comets travel inside of the circle of influence. 


On the basis of the familar parabolic property that the angle between the parabolic axis and the radius 
vector is twice the angle between the tangent and the radius vector we easily obtain the following equality: 


Guat Bp ae tae OT. (6) 
' Equations (5) and (6) give 
= 28 — nF sin a, | (7) 
Since 
GK vote eee Se ae 
1 +- cos (wo — Xc KCJ) ’ (8) 
we have 
Ry 
q= 2! — cos 28 + (1 —cos 28) £- cos « — 2 sin 28 & sin al (9) 
0 0 
From (7) and (9) we easily obtain 
D(q, ®) 
es (a, B) |= 2o | sin B cos (a — B) j. (10) 


The comet meets the circle of influence with the absolute velocity 


n= V205[1— 8 cos a]. 


2R (11) 
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Inside of the circle of influence the comet travels with the relative velocity v; = Vk —vy; therefore 


Yr COS By = v, cos (a + B,) = v,, cos B; (12) 


v, sin By = v, sin (« + By) = vy sin B —v,. (13) 


It then follows that 


5 / | 
met [—sin (*+ Be) + VY 2—2F- cosa — costa + fr). 


(14) 
Because of the smallness of p/Ry ~ 0.06, | — By, |< 9 fora < 5.2 we replace (14) by 
Vp = VJ E sin (a + 8.) + V2 +2 aS cos (« + 8) — cos? (a + °.)| : (15) 
In determining the Jacobian 
Pedc CA Mikes Bhagat 
“D(a + Be, Ba) Sen + O(a + Be) | 035 (16) 


we shall confine ourselves to zero-order quantities. Then in (16) the term 06/08, can be dropped, because 
By =(a +8,)—«a appears in the equation only with the factor p/Rg- Consequently 


D(a, B) | Uy 
SE | SS a ee 5 (17) 
D(a + fe, Be) Py V2 199 aaa (a + %) — cos? (a + Bs) 
0 
On the basis of (12), (13) and (9) we obtain the following approximate formula: 
Ste cet an asaya 
pee ae raed es (18) 


This a very good representation of the exact solution because if in the derivation we also included terms with 
the factor pP/Ry we would have obtained the formula for q with the factor 


[1+ Se sin ge +0], 


Substituting into (10) the values taken from (12) and (13), we obtain 


| ems =o||t+ isin +2. | sine + 5p costa || 


a) 


| =f 2 sin (« + | ee + sin (a + 2) 


=p 242 cos(a + Bs) — cost(a + Ba) V po (a9) 


When we take (17) into account we finally obtain 


_ Dee, a eee (20) 
D (# +82,B2)| " VW Ro vy 


§3. The angular coordinate of the comet which is moving around Jupiter is given by the formula for a 


conic section: 


12 


P= i+ Foose’ le, 


where P is the parameter and E is the linear eccentricity of the hyperbola. P is determined by the law of areas 
according to 


” v2 Sin? fo 


P=— (22) 


where m' is the mass of Jupiter and k? is the constant of attraction. For motion around Jupiter we also have 
conservation of energy given by 


2 2 4 
Op = k?m' [= ae 7 | ? | (23) 
where A is the semi-major axis of the orbit and 
If 
ke =1— a (24) 


From (21), (22) and (23) we easily obtain* 


[2k%m’ tg2o + 9v2] + V [2k?m’tg2e + pv]? — 4k4m'*tg2@ (1 + tg? ) 


Sill 65 25 
2ev; (1 + tg%e) a 
or for values of 2y which are very close to 360° the more convenient formula 
. 2m! 
sin B, = Oe aes SMG ae ae 5 
Vor, V2 tg? pk®m’ + pv? (26) 


*"tg' = ‘'tan' — Publisher. 
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The comet leaves the circle of influence at the point a + 29 and angle 1 — 6, with the normal. The 
transition to fixed space obeys the following formulas: 


Va COS Py = v, cos (x — B, + & + 2); 


(25°) 
Yq Sin Ps = v, Sin (x — By + a + 20) + vy; (26") 
Va = Vr + 2v,vz sin (w — By + w + 20) + v3. (27') 
Taking as the radius vector from the sun to the comet for the exit point 
EE Oe [1 - a cos (a + 29)] , (28) 
0 
we have, according to the energy integral for captured comets and (27), 
2 ; \ 
1 1 20 v7 + 2050, sin ( — Be + @ + 29) 
7 eT Re mn ase ean taped 9 aac maar agen (29) 


0 


Considering the smallness of the factor 2p /R%, we can in the first term of the right-hand side of (29) replace 
a +2 by 71—By+ +29. We then obtain from (29): 


Vp = — vy Sin (t# — By +o + 29) + 


2ek2 4 1 
+ y% sin? (r—B,+a-+4 20) — a cos (t — By + a + 20)—k? e — =| : (30) 
0 
Similarly, by the law of areas, for captured comets we have 
EVP oy 71 co8 (x — Be + @ + 29) 
4 mn 0 0 : 
2 aes sin (™ — Bz + 4 + 29) oy 


For direct motions Vp > 0, and for retrograde motions vp < 0. 


The curves a(@ + By, By) = const and Vp( a + By», Bz) = const were determined as follows. For a = 2.6, 
285052, 5.2 and vp=— 2.4, —2.2,..., 41.3, +1.4,...+ 2.0, + 2.2 with definite valuesof m= 62+ 
+ +29, v, was determined from (30) and (31). From a previously compiled table of v, with the argument 
m + By + a (15) the corresponding values of ™ + Bp + & were determined. In first approximation 2g was 
determined from the formula 


2p = (mw — By +0 -+ 2p) — (m+ Ba + a) + 282 (7 — fa +% + 29) —(* + Br +2). 


With approximate values of(2¢, 7 + a+ B,) Equation (25') was used 2 obtain psd ee values of 7 — By. 
In the second approximation the differences of 2B, from 360° were taken into account. The feu of ws calcu- 
lations are shown in Figure 1. The horizontal axis represents ™ + @ + B, = a + B" and the vertical axis represents: 
B" = 1 + By with opposite sign. The continuous curves represent a = const; the dashed curves represent Vp = const. 


The author possessed a more detailed graph. The values of the Jacobian 


D (a+ Be, Be) 
Dee eer 


‘ : i rked with crosses represent the number of comets 
are proportional to the areas in the figure. Thus the areas ma p 
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with direct motion for which 3.2 < a< 3.6; 0.6 < vp < 1.0, and the dotted area is proportional to the number 
of captured comets with retrograde motions for which 2.6 < a< 2.8; 0.4< vp < 0.8. This result was obtained 
on the basis of the assumption that a + 6" and 6" are distributed uniformly and that there is no selection of 
observations nor disintegration of comets. By taking into account the other factors in (1) we can make the solu- 
tion more precise. 


§4. The lifetime of a comet depends on several physical factors. First, the evaporation of the gases in- 
creases cometary instability, and the loss of gas can cause a comet to lose its luminosity. 


In a number of articles Vsekhviatskii stated that the brightness of a comet decreases monotonically from 
one appearance to the next, but this has been disputed. Some scientists believe that photometric data on comet 
brightness as well as the theory of cometary luminosity are too incomplete to allow such conclusions. There 
have been calculations of comet lifetimes assuming desorption but no adsorption of gases. S. V. Orlov [6] 
determined the number of Ny of molecules lost by the nucleus in unit time. The average for comets is 


N= 7, (32) 


where s is a constant and r is the distance from the sun to the comet. According to Orlov, in each passage a 
comet loses an identical quantity of its gases. Upon this hypothesis it is easy to calculate the life T of a comet. 
We obtain the following relation between T and the elements of the cometary orbit: 


T, T 
T = s,0°h \ | = skp'ra’h|\ (1 + ecosv)? dv} == Syp*a*h [2 +. e2]71, (33) 
0 


0 


Here sy and sg are constants, 


B. Iu. Levin [7] has derived a formula for the loss of gases from the nucleus of a comet, based on the assump- 
tion that the temperature of the nucleus is inversely proportional to vr. He proposed that we consider the store 
of gases N of a comet to be diminishing in accordance with the formula 


N =N,(1—hy)™, (34) 


where py is the number of revolutions of the comet, and k, is the fraction of its molecules lost in one revolution. 
Assuming ky, to be small (which leads to Orlov's hypothesis), we have, to first-order small quantities, 


T =s3(1 —e)~kel ValR-350° ‘ (35) 


where L is the heat of sorption, R is the absolute gas constant and q, is the perihelion distance of the comet. 
On the average L(R- 350° y= 10; 


Secondly, the life of a comet depends on the attractive force. Comets have been observed directly to 
disintegrate into several parts which subsequently disintegrated into a swarm of meteorites. It is reasonable to 
assume that the disintegration resulted from an attractive force proportional to dt /*. The overall disintegration 
P during one period of revolution To as unit of time is given by 


ee wid lie 
P=a a a = Sa lap la, (36) 


es S4 is a constant. If it is assumed that this "disintegration" is a continuous depletion of the comet, we 
then have 
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P= s40"ls p's, (37) 


where sg is a constant. 


Fig. 1. Isoenergetic and isosectorial lines. The continuous lines are curyes of identical 
values of the semi-major axes; the dashed lines represent identical parameters. Negative 
values represent comets with retrograde motions. 


TABLE 1 Since comets cease to emit light or disintegrate 

—— into swarms of meteorites we do not directly observe 

Ve | q | ip the distribution of Jupiter's comet family which results 
only from capture by Jupiter. In order to determine the 

observable distribution we must obtain the probability 


0.0—0.3 i 33 Ry ey ee ee of the existence of a short-period comet. It is reason- 
ae 0.73 | 4.4 | 0.12 2.0—9 5.0—¢ able to assume that this probability is proportional to 
0.8 0.58 | 1.6 | 0.074 = const the lifetime of the comet. Then in accordance with the 


argument of this section we have three hypotheses for 
the function fg(a, p) corresponding to Formulas (33), 
(35) and (37). 


§5. We have a well-developed theory on the probability of comet discovery f(a, p). But, as N. D. 
Moiseev [8] has pointed out, this theory does not agree with observations. Therefore we shall determine f(a, p) 
and f,(q, w) at the same time by an empirical method. For this purpose it is sufficient to assume that the 
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discovery functions of short- and long-period comets are of the same form and that for very small q of long- 
period comets fg(a, p) = 1. This results from most theories on the probability of discovering comets. In this 
case for small q we see from Figure 2 that fy(q, W) = sq, where sg is a constant which is of no interest to us here. 
We take f,(a, p) directly from Figure 2 and obtain’ Table 1. 

$6, From measurements of the elementary 
areas in Figure 1, from Equations (1), (14), (18), 
(20), (33), (35), (37) and from Table 1, assuming 
fy = sgq we have obtained the results in Tables 
2 and 3. The numbers within parentheses denote 
the numbers of detected short-period comets accord- 
ing to the catalogue of [9], for which i < 10°. Six 
comets were not taken into account because they 
cannot be included in Jupiter's family. 


. Tables 2 and 3 show directly that our theory, 
10 20 IO 40 50g both in the form based on Orlov's hypothesis and 
in the variation based on gravitational disintegra- 
tion, is in good agreement with observations. The 
Fig. 2. Distribution of long-period comets discovered fact that some of the discovered comets are some~ 
since 1700 as a function of perihelion distance q. what outside of the limits predicted by theory, was 
explained by Tisserand [1]. The hypothesis based 
on Levin's proposition represents observations poorly. 


As a form of analysis of our calculations we shall explain why we have obtained very small probabilities 
for comets with retrograde motion. This does not at all result from the fact that we did not use Newton's formulas 
but derived our own formulas (Section 3). The difference in the curves for B” > 0 and B" <.0 gives us an idea 
of the inaccuracy of Newton's formula, according to which the curves should be identical for B" > 0 and B" < 0, 
The existing discrepancy evidently cannot result in an essential difference in calculations of the probabilities. 


The almost complete theoretical absence of short-period comets with retrograde motions results from the 
fact that these comets in most cases have small perihelion distances and disintegrate rapidly or, to a lesser 
extent, because they originate from among the relatively small number of long-period comets which have small 
perihelion distances. | 


Astronomical Observatory, 
Latvian State University Received June 19, 1956 
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ON THE COLLIMATION OF A TRANSIT INSTRUMENT 


G, P. Pil'nik 


It is shown that the collimation of a broken transit instrument varies from 
the beginning to the end of observations and also depends on the zenith distance 
During one hour of observations collimation varies on the average by 0°. 05. 
Collimation differs by 08.2 at the zenith and at zenith distance 55°, Therefore, 
the suggestion that collimation can be excluded by a reversal of the instrument 
during observations of each star is unacceptable. 


For observations with a portable transit instrument it is assumed that by reversal of the instrument, 
collimation is totally excluded. Such an assumption is correct only if during observation the collimation remains 
constant. However, our investigations have established that the collimation changed throughout the observation 
and also with the zenith distance. Data on the variation of collimation with time were published for three transit 
instruments in the proceedings of the Shternberg State Astronomical Institute Vol. 27 (1956). Further studies of 
the instrument APM-1 No. 530001 (d = 100 mm, f = 100 cm) showed a considerable dependence of the collima- 
tion on the zenith distance. 


haces 


R Mo Sy Pade Many 4 A, fig fy gg 
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ae fx ———>" 


Fig. 1. Turns of micrometer screw of APM-1 No. 530001. 


For the calculation of the collimation, we used running observations of the error of the clock. As the star 
passages were registered with identical turns of the micrometer screw, it is possible to calculate from the obser- 


vations the mean distance of selected contacts from the noncollimated line (n. c. 1.), In Figure 1 we see the 


contacts of part of the micrometer screw of APM-1 No, 530001. During observation of the star in the order (we) 


on the tape of the chronograph there were registered the contact TIMES tog, trots» © «» CHO-g00 © ov Ldedy + + 0» 


ti1-9, t1g-4- When the star was observed in the order (ew) the registered contacts were ty7.5, ty7.6,+ + +» tit-g» 
-9y ti2-1- 


tyg-4tr- « 0 tigegs Uigeds + © o» L196 (Figure 1). 
We introduce the notation: 


T'= sates + tiga t--- + hoe + 414 +... 4+ tay + tes) 
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and 
T” = oe (tins t --- bho + tsa «++ beet toa «++ + fips): 


. , sa] 
Then for each star we get, depending on the order of observation, the average times T', and T', or T", and 
T",.. The average distances of the chosen contacts from the n, c. 1. were determined from the formulas: 
we 


, 


hoe = 5 (Te =A) cos 8, 


few = = (Tw — Te) cos 8. 


We note that the passage times were also calculated from the same initial data: 


| ped he ad De 


DO 
Lew = S (Ce + Pe 


The calculation of the errors in the half differences can obviously also throw light upon the errors in the half 
sums, i, e., on the passage times to be determined, 


z | ) | few | z | § | it 
61° 630 228642 10° 46° 188762 
12 68 22.621 5 54 18.762 
28 28 22.654 at 29 | | 18.766 

4 57 22.585 27 29 | | 18.776 
35 24 22.648 55 69,, 18.624 
39 17 22.664 17 73 18.748 
17° 39° 22.589 25° a ta 18.767 

Average 22.657 Average| 18.744 


We give the results of the calculation of collimation of f,,, and fy,¢ from observations of 14 stars on 
August 16, 1955, based on one clock error. In the considered example the observation in the lower culmination 
is marked with the letter H. The considerable deviations in the individual values of foyw and fy. can only be 
explained by the variation of collimation or by the inequality of the pivots, because in the above-mentioned 
article it is shown that it is impossible to explain these deviations by variations of declination or of azimuth. 


It is also impossible to explain them by the inequality of the pivots, which in this instrument does not 
exceed 08, 01, 


In order to show the variation of collimation with time we carried out calculations of few and fwe based 


on the observations of two clock errors in one night. Then we obtained the difference of the absolute mean 
values of collimation: 


Nes = Cena ae (Few)s and YA = (fwe)s as (/we)i: 
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The results of these calculations are given in Table 1. This table also shows the date of observation, the 
temperature of the air, the mean values of fey and fy and their sums few + fwe. As the observation of the two 
clock errors was carried out without interruption, the differences Afew and Afye show how the collimation changed 
on the average in one hour of observation, Table 1 shows that there are no seasonal variations of collimation for 
the instrument used. The change of collimation from the beginning to the end of observation is conclusive. The 
differences Afey are negative and Af, positive. 


TABLE 1 
Collimation Values for APM-1 No. 530001 
i 


1955 dates 
day. * mont te Jew A few Swe A fwe fow + fe 


20.4 8 228552 188826 4A $378 
—08027 +05068 

20.4 + 2.1 22 O20 18.894 41.419 

4.9 + 7.6 Des aye 18.836 41.363 
—0.031 +0.019 

4.5 + 6.4 22.496 18.855 41.354 

20.6 $13.4 22.6418 18.810 41.428 
—0.077 -+-0.036 

20.6 +12.4 220041 18.846 41.387 

1.8 +:22.3 22.687 18.706 44.393 
—0.059 -L0.061 

1.8 +24 .0 22.628 18.767 41.395 

8.8 +18.7 22.675 18.739 41.414 
—().042 -+-0.015 

8.8 +17.5 22.633 18.754 AA 387 

15.8 +49.3 22.699 18.7419 44.418 
—0.056 +0.057 

15.8 +18.0 22.643 18.776 41.419 

16.8 22.0 22.657 18.744 41.404 
—0.012 -+0.039 

16.8 +20.5 22.645 18.783 AA 428 

2.9 +18.1 22.693 18.704 41.397 
—0.043 +0.048 

SY) +17.2 22.650 18.752 41.402 

12.9 +18.6 22.688 18.697 41.385 
—0.029 --0.056 

12.9 +17.4 22,.659 18.753 41.442 

4.10 +10.0 22,976 18.834 44.4410 
+0.008 —0.002 

440 + 9.6 22.584 18.832 41.416 

17.410 + 5.6 22.634 18.772 441.406 
—0.022 : +0.017 ae 

7 4.9 D2non, 18.789 4 

2 us ne 22.620 18.740 41.380 
—0.144 +0.164 

0.42 —19.2 22.476 18.904 | 41.380 

€ 


It is interesting to note that on observing the same stars on different nights, Af was different. Such obser- 
vations were made on August 15 and 16, September 2 and 12, and October 4 and 17 (Table 1). The change in 
collimation was particularly marked on December 5, 1955 (Table 1), It must be admitted that SUD Se day 
the temperature was about 0°C and by nightfall a sharp frost had begun. It is possible that the large differences 
of collimation in the observation of the first and second deviation can to a certain extent be explained by the 
drop of temperature. We are inclined to suppose, however, that purely mechanical displacements a Shier 
prism, and micrometer are of no less importance here, When not in use, the objective of our transit instrument 
was at the bottom and during observation it was at the top. It is possible that transfer of the ae from the rest- 
ing to the working position caused a slow change of the relative positions of the parts of the instrument and 
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o°| 2285907 | — 66 
5 |. 99.603 |F "00 
10 | 22.609 | — 54 
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4B | 99. 916p|  eeobe 
4 50 1 22.746 | +83 


55 | 22.784 | +424 


Fig. 2. Dependence of few on the zenith distance 
for the observations on August 16, 1955. 22.663 ode 


consequently the collimation also changed. If the variations in few and fye were only to be ascribed to a tem- 
perature decrease, then a considerable seasonal variation would be observed. The results of Table 1 show that 
the variation of collimation cannot be explained by a temperature influence alone. It is also quite impossible 
to explain the dependence of collimation on zenith distance as the result of thermal influences on the instru- 
ment. 


In the present example of absolute collimation calculated from the observation of one error of the clock, 
the zenith distances are also indicated. If we plot z and fy, ona graph, we notice that with increasing z, 
few also increases, In Figure 2, few is shown as a function of z, based on the observations of two clock errors 
on August 16, 1955. Similar graphs were drawn based on the data of twelve nights of observations. After draw- 
ing a smooth curve, the values of fey for zenith distances 5° apart were taken from the graphs. The collima- 
tion values taken from the graph (for our example) are given in the small table. 


AE AIS YA 


Deviations from the Average Collimation for Different Zenith Distances (in thousandths 
of seconds) 


dates, 1955 


20.4 — 82) —70| —54) —40) —26} —13} + 4] +17 + 77) +104 
4.5 — 84) —68) —52) —37) —22) — 7] + 8) +22] +37] -+52} + 68] + 84 
20.6 —105| —92) —77) —61) —45| —26) — 5} +49} +45) +76] +415] +158 
1.8 — 78} —70) —60) —50} —37} —24) — 8] +44] +33] +60] + 92) +132 
8.8 — 82) —77| —72| —64) —53} —38} —19} +10) +46] +81] +415} +4153 
15.8 — 64) —61) —58} —50) —40) —27| —12] + 6) +25} +55) 4+ 91] +135 
16.8 — 66} —60) —54) —47) —37| —25) —10) + 8] +29] +53] + 83] +421 
2.9 — 91) —50} —47) —46} —41} —34) —23) — 8] +13] +44] + 88] +153 
12.9 — 78} —70) —61) —52) —40) —28) —13) + 6] +28] +57] 4- 94) 4154 
4.10 — 90) —49} —42) —35) —28} —24] —12) + 1] +20] +44] + 73! +102 
17.40 — 67) —59} —50) —39| —29) —19| — $] -++ 8] -+-28 03] + 77) +101 
Sy 442 — 66) —54) —43) —32} —19} — 8} + 4) +47] +30] +44] + 56] + 70 
average — 73] —65| —56} --46} —35} — 22] — s +19 +31) -- 56} + a +122 
\ 


106 


The collimation difference between the zenith and z = 55° is equal here to 08.187. Such collimation 
variations were also observed on other nights. 


Table 2 shows for each star the differences between few and the mean value of fey, which characterize 
the collimation variation with the zenith distance. 


As is seen from Table 2, f,, increases with increasing zenith distances in all cases. When we plot z and 
fwe ona graph, we notice a reciprocal relation between them, i. e., with increasing zenith distance fye de- 
creases, The collimation at z = 0° differs from that at z = 55° by an average of 08,2, This collimation varia- 


tion is most serious, because the observation of stars is carried out at different zenith distances andthus evidently 
with continuously varying collimation. 


Unfortunately, up to now it has been impossible to determine the collimation variation during the observa- 
tion of a single star. However, there cannot by any doubt that the assumption that collimation for each star can 
be eliminated by reversal of the instrument is seriously in error. 


A further increase in accuracy of astronomical observations with transit instruments is impossible without 
careful study of and allowance for the collimation errors, 


P, K. Shternberg State Astronomical Institute Received July 12, 1956 
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INCREASE OF THE SENSITIVITY OF SOME PHOTOGRAPHIC 
EMULSIONS BY MEANS OF PRELIMINARY SHORT EXPOSURE 


T. M. Muliarchik and K, I, Petrova 


The increase of sensitivity of three types of photographic emulsions by brief 
preliminary illumination before the main exposure is described. 


Films of high sensitivity are usually intended for very short exposures (t*1/1 00°), They are inefficient 
for photographing faint objects through long exposures because of strong deviations from the reciprocity law. 


A well known method for reducing these deviations consists in two successive exposures, the first of which 
is a short bright flash of light and the second a long exposure of low intensity [1, 2]. When registering brief 
processes as well as photographing very weak luminosities, such a double exposure gives an increase in sensiti- 
vity. In the first case, the emulsion must be submitted to a weak, long, uniform exposure after the main one 
[3]} in the second case it must be exposed uniformly during a short time to an intensive source before the main 
exposure [4]. Such a preliminary exposure can have a strong effect in photographing spectra of the night sky. 


A selection of the most favorable pre-exposure 
conditions was carried out for two types of films — 


D PF-3 and “type 10,” and also for Isopan ISS plates. The 
analysis of the films was done with the sensitometer 
100 FSP-4 GOI (with ordinary light). The films were devel- 
re , oped in a concentrated Chibisov developer, recommended 
H for PF-3 (8 min at 20°C). The preliminary exposure 
60 4 density varied from 0.03 to 0.40, the time from 1/10000° 
ae to 48, 
ne The nonpre-exposed parts of the film and those 
pre-exposed to a certain density were each time simul- 
7 a fanaa taneously exposed under the same conditions. Finally 
< 8 the characteristic curves were constructed; the trans- 
mission logarithms of the step wedge of the sensitometer 
Fig. 1. Characteristic curves of film "type 10" were plotted along the x-axis and the density of the 
for different preliminary-exposure densities. corresponding steps of the wedge minus the fog density 
1) D=0,40; 2) D =0.25; 3) D = 0,10; were plotted along the y-axis. For a pre-exposed part 
4) D = 0,00 (nonpre-exposed part), of a film the fog density is the sum of development fog 


and the density of pre-exposure. 


The characteristic curves of the parts with different pre-exposure densities are represented in Figures 1 


and 2 for"type 10" and PF-3, respectively. In all cases, the preliminary exposure was 1/1008 and the principal 
exposure was 1 hour. 


The analysis of the influence of pre-exposure on the sensitivity of Isopan ISS plates was carried out with 
spectograph ISP-51, because the selectivity of the sensitizing mechanism was also examined. The plates were 


100 200 200 log I 


Fig. 2, Characteristic curves of film PF-3 at different pre- 
exposure densities. 1) D=0.40; 2) D=0,.30; 3) D = 0.15; 
4) D=0,08; 5) D=0,03; 6) D = 0,00 (nonpre-exposed 
part). 


pre-exposed to ordinary light, The preliminary~exposure density varied from 0.05 to 0.40, One half of a plate 
was pre-exposed and after that, spectra of an incandescent lamp were projected on both parts. In order to extend 
the range of densities, several spectrograms for different slit widths were projected on each half of the plate. 

A common characteristic curve was obtained by shifting the individual curves along the x-axis. 
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Fig. 3. Characteristic curves of Isopan ISS at different preliminary 
exposure densities. 1) D=0.25; 2) D= 0.10; 3) D=0.05; 
4) D = 0,00 (nonpre-exposed part). 
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The preliminary exposures were of 1/508 and 1/1008 while the principal exposure lasted 108, The 
spectra were measured with a photometer in the wavelength interval dX = 4300-6300A, The characteristic curves 
with A = 6300A are presented in Figure 3 for parts of a plate with different densities. 


The common feature of the curves in Figures 1, 2 and 3 is a lower contrast of the sensitized layer and an 
extension of the underexposed region in the direction of lower intensities. Thus, an intensity which on a nonpre- 
exposed film gives no blackening at all gives an appreciable density on a pre-exposed layer, the magnitude of 
which depends on the density of the preceding exposure. When increasing the preliminary exposure density, the 
sensitization increases, It has a maximum for Dpe = 0.30-0.40. The sensitization is different for different 
densities of the main exposure. For a preliminary exposure density of about 0.30, the gain in sensitivity is: 


For D=0.4—0.5 jn PF -3—2 times,in ISS — 1.5 times 
De Ors » on eet 2. » 
Dy LSM! » By 3a 
» D=0.05 — 0.03 » 8 ou 10 
» D<1 pre-exposure gives a loss in sensitivity 


Figure 4 shows the relation between the sensitization and the length of preliminary exposure, The prelim- 
inary exposure density of 0.30 was obtained by exposures of 1/1000S, 1/1008, 1/258, 18, 48, Pre-exposure of 
1/1008 and 1/258 produces the greatest effect. An exposure after the main one gives no increase in sensitivity 
of the emulsions. 


100 200 300 log 


Fig. 4, Characteristic curves of film PF-3 for Dy, = 0.40 
for different pre-exposure times. 1) t =1/10008; 2) t = 

= 1/1008; 3) t= 1/258; 4) t =18; 5) t= 48; 6) the non- 
pre-exposed part. 


For PF-3 film with preliminary exposure density and time of 0.1 and 1/1008 respectively, the main ex- 
posure varied from 1§ to 4 hours, For all these exposures the same gain in sensitivity was pbiaified On a pre- 
exposed plate a density of 0,2 was obtained in 1/4 of the time for a nonpre-exposed one, Messiienrcne of the 
spectrograins on ISS plates did not show the effect to be dependent on wavelength, 


Thus, if the considered emulsions are used for emission spectra of the night sky and other weak objects 
’ 


a preliminary exposure of the film to a density of 0.2- é 
y of 0,2-0.3 at an exposure of 1/1008 allows a sh ing 
exposure by several times, ‘ Hite oe 
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The pre-exposure must be carried out immediately before the principal exposure, because the increased 


sensitivity is retained for only a short time, 
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THE EARLY THERMAL HISTORY OF THE EARTH 
V. G. Fesenkov 


Phenomena directly connected with the thermal history of the earth 
are considered: the flow of granite from magma, the origin and formation 
of continents, the dissipation of the primary atmosphere, the formation of 
oceans, the possible significance of underwater ridges, the motion of the 
poles on the earth's surface, the internal structure of the earth from seis- 
mological data and the earth's magnetism, the distribution of radioactive 
elements in the interior of the earth and the possible production of heat by 
radioactive processes according to data of different authors. Taking these 
phenomena into account,an attempt is made to represent the initial state 
of the earth and its formation from a primary nebula of gas and dust, 
surrounding the sun. 


Considerations are given in favor of the loss by the earth of a con- 
siderable part of its mass, mainly light elements, during its formation. 


Reasons for the differences between Mars and the earth are given. 


The problem of the early thermal history of the earth is shared equally by geology and astronomy and must 
be studied by both sciences. Geologists must provide facts concerning the observed characteristics and structure 
of the crust and the inner shells of the earth, Astronomers can relate the earth's thermal history to the existing 
data on the formation of stars and planets of our galaxy and also, in a measure, draw comparisons between the 
earth and other planets of the solar system. It is essential in these studies to adhere strictly to facts and avoid 
preconceived opinions, which are especially dangerous in cosmogony, which has not yet freed itself from a priori 
hypotheses. First of all, then, the geological and astronomical facts pertinent to the problem must be summar- 
ized. 


1. The existence of great granitic masses of which the continents are composed makes it necessary to 
discuss the origin of granite. Granite consists essentially of three minerals [1], quartz, feldspar and mica. It is 
the product of slow cooling of magma which can occur only at some depth; it has a characteristic texture and 
a low melting temperature. Many attempts have been made to duplicate this texture in laboratories, and differ- 
ent hypotheses have been advanced to explain its origin[2]. Without going into details of these investigations, 
it will suffice to say that apparently the only way in which granite can form is by crystallization from hot 
magma, Other theories of the origin of granite, by metasomatism of sediments or by migration of ions through 
solid rocks, are unacceptable, according to O. T. Tuttle. Granite, which constitutes the bulk of the outer shell 
of the earth~sial—must have formed at a temperature not lower than 1000°C. This requires explanation in the 
discussion of the early thermal history and evolution of our planet. As noted by P. N. Kropotkin [3], according 
to D. S. Beliankin granite may form as the result of a prolonged heating of siliceous rocks which, remaining in 
the solid state, exude a combination of oxides corresponding to the fraction with the lowest possible melting 
temperature. For this reason the melting temperature of granite is about 600°C lower than that of the ultra- 
basic rocks, Sometimes the opinion is expressed that granite may have formed from sediments fused at great 
pressures. Such is the view, for example, of D. I. Shcherbakov [4]. V. A. Obruchev ([1], p. 196) holds, on the 
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contrary, that this is the mode of origin of granitic gneisses. In any case, there is no doubt that granite is the 
product of slow crystallization of melts, 


Since granite is the main material of continents, the problem of its genesis is intimately associated with 
that of the origin and gradual differentiation of the continental masses. 


2. The origin of the continents must also be discussed in connection with the general problem of the 
origin of the earth. 


It has long been known that continental rocks differ in composition from those underlying ocean basins. 
The former are mainly granites, the latter are the denser basalts. The boundary surface between the crystalline 
basalts and the underlying, still denser noncrystalline material, the so-called Mohorovicic discontinuity lying 
at the base of the crust, is much farther from the surface under the continents than beneath the oceans, Accord- 
ing to E. Bullard [5], this discontinuity is due to change in composition of the material and not to a polymorphi¢ 
transition. This is known from experiments in which the temperatures and pressures corresponding to the actual 
depth of this materiai were duplicated. 


It is clear that the sialic continents, being lighter, must be isostatically elevated relative to the general 
surface of the earth. Thus, the heavier basaltic foundation must necessarily be occupied by oceans. 


The origin of the continents is by no means clear as yet, but some clues to it may be obtained from a 
study of intermediate cases, namely, island chains. The rocks of these archipelagoes are mainly granitic but 
they are underlain by relatively thin layers of basalts covered with sediments. As noted by P. N. Kropotkin [3], 
the view is widespread that in the irreversible evolution of the earth's crust a gradual growth of the relatively 
stable continental platforms has taken place at the expense of geosynclines. In the more active geosynclinal 
belts, intensive folding develops, and igneous activity brings magmas close to the surface. These processes are 
active mainly near the platform borders and lead to the gradual growth of continental masses by the addition of 
newly formed volcanic island chains. Of course, some parts of the continents may undergo local transgressions, 
but the general trend of the processes remains unchanged — the granitic continental mass increases in volume and 
area. The irreversibility of the evolution of the crust was very convincingly demonstrated by V. A. Magnitski [6] 
(see also [7]). 


Similarly, M. Kay [8] has pointed out that the island chains are very probably embryonic continents and 
illustrated this by the example of North America. Near the eastern and western margins of this continent there 
are volcanic chains which at one time were archipelagoes but later were joined to the continental block. In 
the regions of present day California and Florida there existed in the past only separate volcanic islands. 


It is believed by some that in the beginning the entire surface of the earth was covered by the ocean, 
and that from it gradually rose.separate islands which grew by addition of sialic material slowly formed in the 
depths of the earth and then brought to the surface. 

The rise of magma from the depths, the development of platforms and the gradual diminution of geosyn- 
clinal belts represent an irreversible slowing down process. It is characterized by gradual diminution of fissure 
volcanism as well as volcanism in general throughout the history of the earth. In the past, as A, P. Pavlov [9] 
so vividly described, great lava fields were formed, over millions of square kilometers; they can be seen now 
on the Yenisei and Tunguska rivers, in Greenland, at Spitzbergen, in India, in parts of North America and else- 


where on the globe. 
This process was accompanied by emanation of water vapor which contributed to the filling of the oceans. 
The development of continents must have been most intensive in the earliest epochs of the earth's history. 
That enormous masses of silicate melt rising from the depths of the earth contributed to it is evident, but the 


dynamics of the process is by no means clear. 


V. V. Belousov has pointed out that formation of granite occurred in the main in pre- geologic time [10]. 
In his opinion, the process of continental growth as described above is acceptable for the Pacific Ocean but not 


for the Atlantic and Indian oceans, which are much more recent features [11]. 


3. In order to discuss the earliest history of the earth it is necessary to consider the problems of formation 
of the atmosphere and the ocean basins. There is no doubt that the present atmosphere is secondary. 
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It is known that oxygen was produced by green plants and to a large extent by oceanic plankton; nitrogen 
is also partly due to the activity of plants but mainly it came from the depths of the planet from which, to a con- 


siderable extent, it comes now. 


Only a few traces remain of the primitive atmosphere of the earth: they are the negligible quantities of 
inert gases such as neon, krypton and xenon, which do not enter into combinations with other elements, It is 
believed that helium, which is mainly a product of radioactive decay, can not remain near the earth's surface 
but is dissipated into inter- planetary space because of the high temperature of the layer of dissipation at the 
height of 400-500 km. 


A question may be asked: What process was responsible for the dissipation of the earth's primitive atmos- 
phere which affected even such heavy gases as krypton and xenon? It is evident that if the gases were dissipated 
separately in accordance with their atomic weights, the heavy ones, very abundant in the universe, would have 
remained on the earth in their original concentration, It must be assumed, then, that the process of dissipation 
was determined essentially by such a light gas as hydrogen, the main constituent of the primitive atmosphere. 
Escaping into space, hydrogen carried with it all the minor constituents; it was responsible for the dissipation 
of the primitive atmosphere, and this process did not require high temperatures. The escape of the heavier ele- 
ments, for example of calcium, aided by hydrogen is observed on the sun, 


It may be supposed that water vapor, if it existed in any quantity in the primitive atmosphere, also escaped 
during this process of dissipation. 


The process of dissipation was not restricted to the primitive atmosphere. It began much earlier during the 
time of formation of the earth. This concept is not based on a hypothesis. The main fact which speaks convin- 
cingly of the loss of a considerable mass by our planet, mainly in the form of light and uncondensed elements, 
is the existence of the moon, as will be considered in greater detail later. 


Where did the water of the modern oceans come from? It must have been derived from the depths of the 
earth during the period of its cooling. 


Every volcanic eruption is accompanied by great quantities of water vapor which immediately falls in 
torrential rains, During the extremely intensive magmatic activity which accompanied the formation of contin- 
ental masses, the surface of the earth received vast quantities of water which gradually formed the oceans, 
Graton [12] believes that the content of gases in extrusive magma may be as high as 10% by weight, and Fried- 
lander [13] thinks that it may be much higher. 

Goranson has shown [14] that granitic magma under high pressure corresponding to the depth of 10-20 km 
may contain great amounts of water in solution, up to 10% by weight. With decrease of pressure, the capacity 
of magma to hold water in solution decreases rapidly and it separates as water vapor. Separation of water vapor 
from cooling lavais a familar phenomenon. Therefore, every rise of sialic material to the surface led inevitably 
to separation of vast amounts of water, which remained on the surface of the earth and filled the spaces between 
the relatively high continental masses, 


The average depth of the ocean uniformly spread over the entire surface of the earth would be 4-5 km, and 


this mass of water must have been derived from the outer layer or crust of the earth approximately 20 to 30 km 
thick. 


It is doubtful that this could have taken place on an initially cold earth only gradually warmed by radio- 
genic heat. 


The presence of the oceans as we now know them served, therefore, as a strong argument for an initially 
hot earth cooled by intensive convection currents. 


According to Graton, even modern volcanism points to the existence of a primitive volcanic source related 
to tectonic displacements and lying at a depth much greater than 100 km. 


It is interesting that volcanoes, both ancient and modern, are not isolated from each other but occur in 
groups covering large areas. In each group, individual volcanoes are alike in type and petrography, although the 
country rocks within the limits of the group may vary. Graton considers this to be evidence in favor of the 
existence of substratum of uniform composition which circles the globe (see p. 93 [12]). 
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4. Very interesting features which it may be supposed have a direct relation to the thermal history of 
the earth are: deep trenches, valleys: and mountain ranges which are distributed in regular sequences on ocean 
floors and apparently remain unchanged for long periods of time. 


Besides these mountain ranges, there are deep canyons on the ocean floors, sometimes having numerous 
tributaries, such, for example, as the canyon which starts near Greenland and extends south along the bottom of 
the Atlantic ocean to the 38th parallel. 


The question of their origin was long ago posed by Academician L. D. Berg. At present it has been proved 
on the basis of new data that such canyons are formed by gigantic turbidity currents initiated by slides on the 
margins of continental platforms during certain times in the earth's history, as for instance at the end of the Ice 
ages , resulting from rapid melting of ice, or during strong earthquakes. Under normal conditions the velocity of 
submarine currents diminishes with depth and becomes negligible [15]. 


The enrichment of oceans with different salts necessary to prepare a medium for the inception of life could 
have occurred only with the participation of the continents, whether we consider the great water cycle or nature 
of the submarine currents. 


More than ten deep trenches have been mapped in the Pacific and these may be intimately related to 
the origin of the continents. The general tectonic character of the oceans was described by Belousov [16]. 
Academician D. I. Shcherbakov [17] has pointed out that there are many submarine mountain chains in the region 
of the Antarctic continent forming a complex system which ties it to South America and partly to Australia. 
Huge mountain chains are separated by deep troughs, such as the South Hawaiian depression with its depth of 
8264 m, for example. 


Many of these trenches and mountain chains exceed in size those existing on the earth's surface. There 
are, for example, four great escarpments (Mendocino, Murray, Clarion and Clipperton) which start near the 
shores of North and Central America and trend westward at a considerable angle to the meridians, parallel to 
and nearly equidistant from each other [18]. These and other fracture zones characterized by sharply defined 
mountain ranges must have been formed as the result of tremendous dislocations of the earth's crust. 


According to G. Menard [19], these structures could have been formed by huge convection currents immed- 
iately below the crust. It is assumed that hot materials from the depths of the earth rise and the relatively cold 
materials just below the crust sink. This phenomenon, Menard points out, has never been observed directly, for 
it could have taken place only during the earliest history of the earth, but it is the most likely explanation for 
the existence of these structures. The drag exerted by the convection currents against the lower surface of the 
earth's crust could have produced parallel fracture zones. It may be concluded that in the remote past the outer 
shells of the earth were plastic, a condition which no longer exists. 


5. The distribution of radioactive elements within the earth has an important bearing on its history. 


It is well known that the main radioactive elements, uranium and thorium, are enriched in the uppermost 
layers of the earth some tens of kilometers in thickness and that the abundance of these elements is greater in 
the younger rocks. V. G. Khlopin [20] presented convincing data on this subject. Meteorites, the analogues of 
the deeper parts of the earth, are only weakly radioactive. 

It is natural to suppose that the primary material of the earth and other planets had a uniform distribution 
of the heavy radioactive elements simply because it was in a dispersed state. This raises the question, what 
process resulted in the present peculiar distribution of these elements. 

There is no doubt that if the entire earth had the same concentration of radioactive elements as its surface 
layers it would never have solidified, and yet seismological data show convincingly that the earth is solid except 
for the inner core. 

The cause of the selective distribution of radioactive elements, according to Goldschmidt [21], is the great 
size of the uranium and thorium ions, which does not allow them to enter into the silicate lattices during the 
process of crystallization of magma. 

Therefore, during the initial stages of solidification of the deeper parts of the earth these elements were 
gradually forced upward, and enriched its upper layers. In this connection, it is interesting to note that, in 
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1928, Pigott [22] had pointed out that in granite the deposition of radium takes place on the surfaces of mineral 
grains and not within them, Holmes [23] has postulated that the volatile matter emanating from the depths of 
the earth carries radioactive elements with it. 


This leads to the conclusion that the radioactive elements were uniformly distributed in the primitive earth 
while it was in a molten state. As it gradually solidified, the radioactive atoms migrated toward the surface and 
became concentrated in the crust. If it is assumed that the earth was cold in the beginning and never reached 
the molten state, the cause of the present distribution of radioactive elements cannot be explained. Their dis- 
tribution is easily understood, however, if it is assumed that in the beginning the earth was molten and was 
cooled by intensive convection currents causing the radioactive elements to migrate upwards as it gradually 
passed into the solid state. 


A detailed discussion of migration of radioactive elements was given by V. V. Belousov [24]. 


6. Perhaps the problem of the wandering of the poles should also be mentioned. Evidently no noticeable 
shifting of the poles occurred during the Tertiary. This was emphasized in various papers during the symposium 
at Cambridge devoted to this problem, and especially by Professor Dershem (University of California) on 
the basis of his study of fossil marine fauna. of different epochs and by Dr. Arnell (Cambridge) and S. Runcorn, 
on the basis of the paleomagnetic studies [25]. According to them, several reversals of polarity occurred during 
the Tertiary but no change in position of the poles. 


It appears probable, however, that before the Tertiary the average position of the magnetic field changed 
slowly with time. Runcorn explains this by slow shifting of the earth's axis. 


Other specialists believe that the evidence is still insufficient for them to accept the wandering of poles 
as a fact since in the past, as now, the Arctic Ocean was characterized by low temperatures, Corals did not live 
there and deposition of limestones was relatively slight. 


If, however, wandering of the poles is accepted, it must have occurred in distant geologic times when 
large movements of masses within the earth were still possible. A theoretical analysis of this latter phenomenon 
was made by the astronomer of the Greenwich Observatory, Gold, but factual evidence is scant. Undoubtedly, 
shifting of masses within the earth and changes of relief on its surface must affect the position of the earth's 
axis. 


7. Of great importance to the discussion of the early history of the earth is the knowledge of its internal 
structure obtained by seismological methods in the past few decades. The basic facts are that the earth is solid 
throughout its entire mantle, from the core to the surface, and that the mantle consists of at least seven shells 
with densities changing abruptly from one to the next. These shells, according to Bullen, differ in chemical 
composition and can be only partially explained by polymorphic transitions. The outer shell is composed of 
olivine, In the deeper shells change takes place from one modification of olivine to another. Still deeper, 
heavier compounds dominate (Si — Mg — Fe), and beneath that lie shells of a variable composition. Both longi- 
tudinal and transverse waves pass freely through the mantle. Longitudinal waves do not pass through the core, 
and this is the main proof of its liquid state. It has been known since 1936 that the core itself is not homogeneous 
but consists of at least two concentric parts. The inner and denser part may be solid. According to Bullen [26] 
(see also [27]), it consists of iron-nickel, possibly with an admixture of denser materials. 


As long ago as 1939, V. K. Lodochnikov suggested that with sufficiently great pressure, materials of the 
earth must pass into a metallic state stable in that environment. A more detailed investigation of this was made 
by Ramsey [29, 30]. The metallic state is marked by high electric conductivity because of weak bonds between 
atoms and electrons, This is accompanied by a sharp increase in density (for example, density of hydrogen in 
the metallic state increases sevenfold, and helium becomes 10-20 times denser). The metallic state involves 
destruction of electronic shells and of saturated molecular bands, and is characterized by a greatly increased 
electrical conductivity. There is no need, therefore, to think that the core is metallic in the usual sense, that 
it consists of a substance like an iron-nickel alloy which behaves like a metal under normal conditions of temper- 
ature and pressure. The change from the mantle to the core may simply be a change of state and not a chemical 
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change, although either is possible. Since theoretical calculations related to the change into the metallic state 
are extremely difficult, although not impossible in the sim 
paring the earth with other planets. 
and oblateness of Mars indicate that 


pler cases, the final solution may be reached by com- 
It is possible, according to Ramsey, to show that the observed mass, size 

it is composed of the same materials as the earth if the difference between 
the two planets arising from a different distribution of densities with depth is taken into consideration. 


On the other hand, Jeffreys has shown that the concept of similarity in chemical composition between the 
earth and other planets is invalid if the core of the earth has a composition substantially different from the rest 
of it. Admitting that all planets are composed of the same material modified by the existing pressures it is 
still possible to assume a considerable admixture of the heavy elements in the central parts of our planet. This 
is not based on observed facts but appears probable in view of the known differentiation of material in the mantle 
and is suggested also by the striking difference between the metallic and stony meteorites which are undoubtedly 
remnants of a planet or of several planets which once existed in the solar system. 


8. The magnetic field of the earth provides evidence for judging of the properties of its interior, That the 
secular variations of the magnetic field of the earth originate entirely within the earth has already been shown 
by Gauss. According to the modern view the magnetic field of the earth is generated by electric currents within 
its liquid core [31]. This theory was worked out in 1939 by El'sasser [32]. The existence of the magnetic field 
shows that the core is liquid and not merely plastic, since varied and intensive movements occur within it. The 
change in polarity repeatedly observed in "fossil" magnetism of the Tertiary suggests that these currents may 
die out momentarily and then flow again but in the opposite direction. 


All this points to the high temperature of the earth's interior, which must have remained practically un- 
changed since the beginning, partly because of low thermal conductivity of the earth's mantle and partly because 
of lack of exchange of material between the core and the mantle. 


These are conclusions concerning the internal structure of the earth based on an analysis of existing data. 
It is possible to show that they are not compatible with the concept of an initially cold earth. On the contrary, 
very high temperatures must have existed on our planet from the very beginning. 


Let us consider this in detail. 


Four basic sources of heat existed during the evolution of the earth. One is the heat generated when a 
compacted body forms from a dispersed mass. If the process of agglomeration is rapid then not less than one- 
half of the kinetic energy developed from the potential energy during the process of compaction is heat. 


The second source is radiogenic heat, which accumulates very slowly and whose effect becomes significant 
only over long periods of geologic time. When generated within the inner portions of the earth, nearly all of it 
is used up in heating, for little can escape upward. Thermal conduction becomes perceptible only at small 
depths below the surface. The effectiveness of this source of heat depends on the law of distribution of radio- 
active elements within the earth. 


The third source of energy during formation of planets derives from chemical reactions. The importance 
of this source during formation of cosmic bodies has been repeatedly emphasized by Urey [33, 34]. 


To return our planet to its initial dispersed state, it would be necessary not only to overcome the gravita- 
tional forces holding it together but also to break chemical bonds of various complex molecules in order to 
obtain simple molecules and atoms of the protoplanetary cloud. The energy required for this would exactly 
equal the energy expended during the formation of the earth, 


The birth of our planet involved not merely the aggregation of the originally dispersed matter but also 
formation of very complex molecules. 

Finally, the gravitational differentiation of materials of the earth must have had a certain effect on its 
heat budget. 


It can be easily shown that if gravitational condensation were the only process involved in the formation 
of the earth, even then sufficient amounts of energy would have been released to melt the entire globe. Ina 


heated mass convectional currents would be readily started and lead to a rather rapid cooling. During this 


process radioactive elements would be carried towards the surface and become concentrated in a relatively thin 
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layer. Thus, radioactive materials migrated toward the surface before the effects of radiogenic heat became 
perceptible in the interior and their influence on the heat budget of the earth as a whole became negligible. 
Concentration of radioactive elements in a thin surface layer affects the geothermal gradient and in case of 
unusual local concentration may possibly result in local fusion, Assuming the average geothermal gradient to 
be 30°C per kilometer, we find that the maximum temperature at the base of this layer cannot exceed 1000°C. 


Similar temperatures must obtain throughout the earth, if it is assumed to be without any other sources 
of heat. With the tremendous pressure within the earth, such temperatures obviously can have but little effect. 


In the past, when the amounts of radioactive materials were much greater, the geothermal gradient was 
much steeper and the temperature at the base of the radioactive layer was higher but still inadequate to cause 
melting of rocks at the base of the mantle where pressures of 1,3 million atmospheres exist. 


Therefore, the exact way in which the earth was formed from a dispersed cloud becomes of primary im- 
portance in tracing its thermal history. This brings us to a discussion of the mechanism of formation of the 
planets. 


According to the well-known cosmogonic theory of O. Iu. Schmidt [35], the earth and the other planets 
were formed from a cloud of meteoric dust by gradual agglomeration and welding together of individual 
particles colliding with one another, The original part of Schmidt's concept is in his postulate of capture of 
such a dust cloud by the Sun. This capture of a rather massive cloud has not been theoretically discussed and 
appears to be quite impossible. Moreover, such a mechanism is unnecessary, since the planets were formed not 
from a cloud consisting of dust alone but from an extensive and rather dense mass of gas in which dust particles 
constituted a minor admixture, There is no need to cite arguments for this concept, since no one is likely to 
deny it now. 


Essentially the process of planet formation differs from that of the formation of stars only in scale, just 
as planets differ from stars merely in having much smaller masses. 


There exist several calculations of the rate of heating of the primitive earth based on the postulate of its 
origin from a cloud of dust by accretion. Such a calculation in strict conformity to Schmidt's cosmogonic hypo- 
thesis was made by V. S. Safronov [36]. 


The rate of growth of an embryo planet depends on the frequency of its encounters with particles of dust 
and on the density of the cloud. Safronov calculated the rate of growth of the planet and as a corollary the rate 
of heating of its surface. The temperature of the surface for every given moment was determined by consider- 
ing the balance between the heat generated by the fall of particles and the radiation of heat into space accord- 
ing to Stefan's law. 


Having made probable assumptions as to the speed of the dust particles, he found that the temperature of 
the primitive earth during its formation could not have exceeded 1000° abs, at the depth of approximately 
1000-1500 km, 


The time required for the growth of the earth to its present size under these conditions was given as 20 
million years. This is a short period of time compared with the age of the solar system, but it should be noted 
that the duration of solar radiation in the process of purely gravitational contraction of the Sun from an infinitely 
large size is also approximately 20 million years, according to well-known calculations. The relatively huge 
mass of the sun falling upon itself from infinity was capable of maintaining radiation without aid from atomic 
reactions for only this extremely short period of time. 


The arguments of Safronov are not very convincing, since energy is dissipated not only during the process 
of accretion but during an indefinite period of time after the fall of a particle on the planet's surface and this 


must lead to an incomparably greater lowering of temperature. Thus, the result obtained by Safronov must be 
regarded as a certain maximum value only. 


The author of the meteoritic hypothesis, Schmidt, considered the process of accretion insufficient for heat- 
ing the planet (see Four Lectures on the Origin of the Earth[37]). In general, it was believed that in the light 
of this hypothesis the earth must have been completely cold in the beginning. Such was the opinion of P. N. 
Kropotkin [3]; E. A, Liubimova, in her calculations of radioactive heating of the earth, based on Schmidt's 
views, assumed in her principal work [38] that the initial temperature of the earth was about OC. 
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It is clear that the mechanism of formation of the earth Proposed by Schmidt does not lead to a noticeable 
increase in temperature. In any case, this mechanism could not have resulted in the melting of the whole earth 
or of a large part of it, including the core which is still liquid, and therefore intensive convection currents could 


not have been generated within the earth. Analogous calculations were made by T, Gold [5] on the assumption 
that the earth grew by accretion. 


All these calculations show that even in this process of formation of planets a certain amount of heating 
will take place, reaching perhaps 1000°C, but this is quite insufficient to melt the earth and especially its inner 
core and to insure rapid cooling. Therefore, in order to explain the present high temperature of the earth and 
especially of its interior, it was necessary to turn to a second source of heat, radioactive decay. A number of 
works by E. A. Liubimova in the Institute of Geophysics, Academy of Sciences, USSR, and also by some foreign 
authors (39, 40], are devoted to this problem. 


E. A. Liubimova, like Safronov, starts with Schmidt's hypothesis by assuming that the primitive earth was 
cold when heating by radiogenic heat began, and implies a uniform distribution of radioactive elements throu gh- 
oug the earth. This is a highly artificial approach since it makes it necessary to assume that the content of 
radioactive elements in the earth is as low as in meteorites and lower than that observed even in the deepest 
terrestrial rocks. This concentration corresponds approximately to that which the earth would have if all radio- 
active substances were concentrated in a layer 200 km thick. 


Liubimova's calculations show that the thermal field of the earth can be divided into two regions: the outer, 
where heat flows towards the surface and is partly radiated into space, and inner where temperature increases 
continuously because of generation and conservation of radiogenic heat. The depth of the layer separating these 
two fields is of the order of 6Vkt (k is the thermal conductivity; t is the temperature) and increases slowly witht. 


At present, this depth is about 1000 km (k = 0.01). 


These calculations show also that in the interior the temperature continuously increases, while in the outer 
region it decreases. The maximum temperature which could be obtained in the course of the earth's evolution 
from radiogenic heating of the lower part of the mantle at the core boundary is about 4000°C, Cooling of the 
outer part of the earth causes a pronounced contraction to the depth of 700 km, and this gives rise to stresses and 
earthquakes. From this point of view, in spite of the postulated distribution of the sources of radioactive decay, 
the earth has always remained solid. However, occasional local concentrations of radioactive substances could 
have led to local melting, to pockets which served as sources of igneous activity. 


A number of known facts concerning the internal structure of the earth make it necessary for Liubimova to 
abandon the concept of an initially cold earth and to seek new factors which would insure sufficient heating. For 
example, in order to explain the concentration of radioactive substances in a thin surface layer, it is necessary 
to assume that the entire earth was heated and brought to a plastic state so that convection causing migration 
of radioactive substances to the surface became possible. However, this convection must occur only after a 
complete fusion of the core, since the extreme slowness of radioactive heating and the rapid escape of these 
substances would not allow time for heating the inner portions of the earth to a sufficiently high temperature. 
Such artificial constructions must be introduced merely because Liubimova is committed to Schmidt's concept 
of formation of the earth by slow aggregation of dust particles, a process which does not provide the primitive 
planet with the necessary heat energy. 


Actually, it must be assumed that the earth and other planets were formed not from a quasi-meteoric cloud 
captured by the sun in passing through a hypothetical nébulosity, but from a gas and dust medium surrounding the 
primordial sun and identical with it in chemical composition. 


It is quite evident from rather elementary considerations that during formation of a star from a gas and 
dust medium only a certain part of the condensing cloud goes into the making of the star; the remainder and 
possibly the bulk of the medium must remain in a dispersed state because too Aa cee eT can- 
not be concentrated in one body without its becoming unstable. Only under exceptional conditions could the 


entire mass of a cloud be gathered into one compact body. As a rule, a new star will be surrounded mainly in 


its equatorial plane by a cloud of gas and dust with a rather high density and a relatively low temperature, 
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Because of gravitational instability this cloud provides more than enough material for secondary condensations 

or planets. The angular momentum which cannot be concentrated in one compact body must necessarily, during 
the advanced stages of condensation, be distributed among separate gravitational centers with their own orbital 
motions, Even then, much of the cloud must remain in a dispersed state. Illustrations of this are numerous 
double and multiple stars in our galaxy and also the earth-moon system, It is almost certain that the moon, 
during its formation, acquired a considerable part of the initial angular momentum and thus allowed the earth 
to complete its condensation. But if this point of view, based on a modern concept of the formation of stars, 

is accepted, an elementary calculation shows that the modern earth represents only a small fraction of its ori- 
ginal condensation. Indeed, if the initial angular momentum and consequently the mass of the earth had not 
been lost to a large extent the combined mass of the earth and moon united into one body with the total angular 
momentum which now characterizes this system would have been unstable and the period of its rotation would 
have been at most 1.4 hours. In reality because of this loss it is about 4.2 hours. 


It may be supposed that the process of formation of a planet, once it started, would go on rather rapidly. 
According to modern calculations the time of condensation of a star from a globule is comparable to time needed 
for a simple fall of an element of mass from its surface to the center. To initiate this process it is necessary 
only that through a preliminary condensation the density of individual heterogeneities which must inevitably 
form within the diffuse cloud satisfied the condition of tidal stability, which for the central parts of the solar 
system , for example, amounted to about 107° g/cm’. 


The temperature of planetary condensation can be estimated from the fact that about half of the potential 
energy is expended on heating it during the process of shrinking. 


In the case of a homogeneous sphere, the potential energy is 


In condensation of a heterogeneous sphere the expenditure of energy increases, For instance, in case of 
heterogeneity satisfying Roche's law, which formerly was widely applied to the earth, namely: 


Saabs [1 oa a(x) | » B =0.38,/ 


we have 


U = 3.6-10° ergs 


In the case of heterogeneity satisfying the exponential law 


8 = Bye A0r1R), 


we obtain 


where 


j (a) = 25 bites stoke nee a ee 
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OR 
and for ee = 0.15 we have 
0 TCM econ Ih 
In the case of modern earth: 


GM2 7 
A = 3-75-10; R= 6400 km; M =6-1027¢. 


All these values of potential energy are of the same order of magnitude and lead to heating measurable 
in tens of thousands of degrees. 


In the Roche case, 1.8-10°* ergs or 1.5.104 cal/gram are available for heating. 


It is quite evident, then, that primitive heat generated during the process of formation of the earth is of 
the greatest importance in its later history. Radioactive heating is very slow and before its effect could become 
significant all radioactive substances would have escaped toward the surface where their role is localized and 
does not affect the entire earth. 


Jeffreys, in the last edition of his well-known book [41], writes, in complete accord with the above, that 
initially the earth was gaseous and could have radiated all of its energy in about 500 years. 


The total energy of condensation of the earth to its present radius is estimated by him as 10,000 cal/g 
and the energy due to initial temperature as approximately 6000 cal/g. In the presence of convection currents, 
the earth could have cooled to a liquid state in only 5000 years after its birth. With further cooling, which took 
place from the outside inward, the total time of cooling and solidification with cessation of convection was 
much longer. Jeffreys assumes the following conditions in discussing the process of cooling ([41, p. 286): the 
initial temperature of the surface is near melting point and gradually decreases to 300°K; the radioactive layer 
is very thin and radioactivity decreases rapidly with depth. 


The formulation of the problem given above agrees well with the modern concepts of science. 


It remains to consider a number of side issues of great interest. First, why does not our neighbor planet 
Mars have water bodies and is poor in water in general? 


It may be asserted at present that Mars can easily hold water vapor without allowing its escape into space. 
As was shown by L. Spitzer [42], in order to allow the escape of water vapor, the temperature of the appropriate 
layer of the martian atmosphere should exceed 500°K, Actually the temperature of the martian atmosphere is 
much lower than this. On earth, the temperature of the high ionosphere increases considerably because of the 
presence of atmospheric oxygen produced by photosynthesis. On Mars, oxygen is absent and no such heating 
effect can take place. The loss of water vapor which comes to the surface of the planet from its depth or exists 
in its atmosphere must occur inevitably in the earliest epoch of its history when the entire mass of the planet 
together with its surface is still hot [43]. However, the presence of only local regions of fusion on the surface 
of a planet is sufficient to effect escape if they are large enough to heat the atmosphere Be them to the zone 
of dissipation, which on the earth is 400-500 km above the surface. Therefore it became possible for Mars to 


retain water only when the formation of its crust was far advanced. 


The very small amount of water existing on Mars at present is evidently the product of volcanic eruptions 
which occurred after the planet had almost completely solidified. The moon, with a mass almost a hundred 
times smaller than that of Mars, differs from it merely in that it has not only completely lost water vapor during 


the process of cooling, but continues to lose it, if (as some selenologists believe) it still issues from its interior. 
, 


The evidence presented above leads to the conclusion that because of the mechanism of its formation, 
the earth (from the very beginning) had considerable energy and therefore high temperature, and that the slowly 
accumulating and gradually diminishing radioactive heat was of secondary importance for the earth as a whole. 
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THE PHYSICAL NATURE OF THE DIFFERENCES BE TWEEN 
GEOMAGNETIC DISTURBANCES OF SUDDEN AND NONSUDDEN 
COMMENCEMENT 


E. R. Mustel' 


The differences between sudden-commencement geomagnetic storms and 
nonsudden-commencement storms are considered. The discussion shows that the 
ejection of condensations, leading to nonsudden-commencement storms is evidently 
caused by some shortlived nonstationary process of electromagnetic nature in the 
region of sunspots. This process is always accompanied by peculiar phenomena 
of radio emission with a shift of the level of radio emission in the corona. It 
can be accompanied by chromospheric flares, but does not necessarily cause them. 
Disturbances with nonsudden-commencement characteristics are caused by 
flocculi. This is confirmed by the fact that the passage of flocculi through the 
center of the sun's visible disk is accompanied by geomagnetic disturbances _ 
Liseeels 


The difference in the character of sudden-commencement and nonsudden- 
commencement disturbances must be connected with differences in physical 
properties of the streams and not with their geometrical properties. Condensa- 
tions causing sudden-commencement disturbances have larger velocities than 
those causing nonsudden-commencement disturbances; this influences the 
rapidity of the variations of the Earth's magnetic field (Formula 1). 


Besides, condensations, causing sudden-commencement disturbances and 
ejected from the region of sunspots, have their own considerable "frozen" mag- 
netic field. This should also cause additional variations of the Earth's magnetic 
field. Considerations show that the field of the condensations, causing nonsudden- 
commencement disturbances should be much smaller than in the first case. 


Geomagnetic disturbances constitute one of the main sources of information on the corpuscles emitted 
from the sun. In the first place, such disturbances are basically corpuscular in character (except for such phen- 
omena as the "hooks" observed at the time of chromospheric flares), In the second place, we have at our dis- 
posal singularly uniform data collected over many decades and pertaining to geomagnetic disturbances. 


It has been determined by now that geomagnetic disturbances can be divided into two groups. The first 
of these consists of disturbances of sudden commencement, the second — disturbances of nonsudden-commence- 
ment. It is known that the intense disturbances which usually accompany brilliant chromospheric disturbances 
are characterized by sudden-commencement. However, the weaker geomagnetic disturbances also include those 
of sudden commencement; thus, weak and medium-intensity disturbances are subdivided into nonsudden-commence- 
Ment ones (these are in the majority) and sudden-commencement ones (the minority). It is very important to 
note that these two types of disturbances differ not only in the character of their commencements but also in 


many other characteristics. Considering at first only the weak and medium disturbances, we will list the basic 
differences between the two types, 
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a) They are entirely different as to the character of their numerical distrubution in the cycle of solar 
activity. The number of the sudden-commencement disturbances follows the progress of sunspots in the cycle 
quite closely [1]. On the other hand, the curve corresponding to the number of the nonsudden disturbances is 
displaced by some three to four years in respect to the curve corresponding to the number of sunspots (the dis- 
turbances occurring later than the sunspots), and is altogether different in character (ll. 


b) The nons = i i i i i 
P . 4 onsudden PsoIumencement disturbances are characterized by their 27 day sequence, which is 
specially pronounced prior to the minimum of solar activity. On the other hand, the sudden-commencement 
disturbances only rarely occur in such short sequences; note especially the work [2]. 


c) The most intense of all the weak and medium disturbances are of the sudden-commencement type. 


d) The two types of disturbances differ, as is indicated by their very name, in the nature of their begin- 
ning phase. The same can also be said of their terminal phase [1]. 


e) The duration of the nonsudden-comimencement disturbances is sometimes very great and may extend 
up to 10 days, while the sudden-commencement disturbances last usually from 24 to 36 hours. About 10 per 
cent of the latter disturbances last over 36 hours, and some of them are shorter than 24 hours. 


f) The disturbances differ also in the character of their progress; those with the nonsudden commencements 
exhibit weak oscillations, usually with different phases of hydrogen and deuterium oscillations. The sudden- 
commencement disturbances oscillate rapidly or very rapidiy, and the character of the oscillations is reversed 
from time to time. Furthermore, the sudden-commencement ones* often exhibit this storm-time variation, 
while the nonsudden ones do not show this feature. 


g) Streams comprising the nonsudden disturbances are characterized by high-intensity radio emission. 
The following facts corroborate this: 


1) A number of nonsudden-commencement disturbances exhibit a sharply defined seasonal regularity, 
with spring and fall maxima [8, 4]. 


2) These maxima are most sharply defined in years of minimum solar activity. 


3) The curve representing the number of nonsudden-commencement disturbances is displaced (in 
the direction of retardation) in respect tothe curve representing the number of sunspots. This is indicative of the 
importance of the low solar latitudes and, at the same time, the radio emission of the streams, There are also 
a large number of other factors indicative of radio emission. 


The characteristics listed above pertaining to weak disturbances of sudden commencement are even more 
emphatic in the case of intense disturbances of sudden commencement. Thus, sudden-commencement distur- 
bances generally differ in many respects from those of nonsudden commencement. 


Everything that has been said above is indicative of the fact that the solar origins of the streams compris- 
ing nonsudden and sudden-commencement disturbances are altogether different. The following facts are expecially 
important in arriving at this conclusion: 1) the completely different character of the distribution of the two 
disturbance types in the solar cycle; 2) the difference in the character of the 27-day sequence, and consequently, 
in the persistence of streams in the two types; 3) the difference in the character of radio emission of the streams; 
i. e., geomagnetic disturbances related to bright chromospheric flares are produced even when a flare occurs at 
a distance of up to 45° from the center of the disc. This means that in this case the band of corpuscles is very 
broad. We may assume (see below) that the same is also true for the streams comprising weak disturbances of 
sudden commencement. 4) The very character of the disturbances is entirely different. 


In this article we propose to investigate thoroughly the problem dealing with the origin of disturbances of 
both types. Let us consider the sudden-commencement disturbances. It follows from the facts listed above that 
these must be closely related to the sunspots. It is obvious, however, that the sunspots themselves are not the 
source of the corresponding streams, Also, in the case of intense geomagnetic disturbances of sudden commence~ 
ment the source of the corpuscles is normally in the bright chromospheric flares, which in themselves are very 


*We have discussed here the weak disturbances. The storm-time variation is nearly always present in the case 


of intense disturbances of sudden commencement. 
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closely related to the spots [1]. Furthermore, observations of solar radio emissions indicate directly that the . 
solar corona at the time of a flare exhibits an outward-from-the-flare movement of some disturbed agent which, 
from all the indications, represents a flow of atoms away from the flare. It is possible also to conclude, on the 
basis of the same radio-emission observations, that atoms in the investigated stream have a very complex veloc- 


ity spectrum. 


All of this refers to the intense geomagnetic disturbances. The matter of the medium and weak disturbances 
of sudden commencement is not quite clear. It would be possible to assume that these, too, are produced by the 
flares but with less intensity. Thus, for example, Newton [5] claims (see, in particular, Figure 2c of his article) 
that those spot groups within which chromospheric flares are most commonly produced are geoactive, At the 
same time. a number of data indicate that sudden-commencement disturbances can be observed without any 
apparent connection to chromospheric flares. This allows us to made the following suppositions. 


The latest investigations make it possible for us to conclude that the outflow of atoms from chromospheric 
flares is not conditioned by the occurrence of the bright-line emission of the flare spectrum, but, on the contrary, 
that both these phenomena (that is, the emission of atoms and the appearance of a bright flare on the disc, or 
at its edge) represent the results of the same extremely complex process which is, apparently, of an electromagnetic 
nonstationary character,* In fact, if the emission of atoms from the flares were conditioned by the bright radia- 
tion of the flare, then the only mechanism capable of producing such an emission would be that of the light 
pressure. But the author has pointed out in one of his recent articles [6] a number of reasons making the accept- 
ance of this concept difficult. There are also other even more serious objections to it; study of the velocity 
spectrum of particles emitted from the bright flares (see above) indicates that light pressure can be of no con- 
sequence here, because the spectrum includes, besides the velocities on the order of a thousand kilometers per 
second, considerably higher velocities of the order of several thousand kilometers per second, not to mention the 
emission of particles with the energies of cosmic rays. 


But if all this is true, it is quite possible that the indicated nonstationary electromagnetic process is not 
necessarily accompanied by the appearance of a bright chromospheric flare on the disc. In other words, it is 
quite possible that in a number of cases (or in a great many cases) the activity of this process does not appear in 
all the possible forms, which is to say that this process may be accompanied by the emission of corpuscles, 

_ but may not be accompanied by the appearance of a chromospheric flare. Opposite cases are also possible. 
Indeed, observations of chromospheric flares and of radio emission confirm the complicated nature of the whole 
picture. Further proof of the fact that in all the cases of corpuscle emissions we deal with the same specific 
nonstationary process (regardless of whether a bright chromospheric flare does or does not appear) is provided by 
the identical character of the intense and weak disturbances of sudden commencement. This fact is of great 
importance. 


It follows that for the sources of the corpuscles producing the sudden commencement distrubances we must, 
apparently, accept those regions in which those originate specific nonstationary phenomena related to the inten- 
sification of radio emission and to the displacement of the level of this radio emission in the corona. These 
phenomena do not have to be accompanied by the bright chromospheric flares. We can only venture a guess 


that the indicated regions are very ephemeral in nature, and that they are related to the sunspots (see, for instance, 
the previously mentioned results obtained by Newton [5]). 


We will discuss briefly the origin of corpuscles producing geomagnetic disturbances with a non-sudden 
commencement, The author believes that the corresponding solar regions coincide with the regions of the floc- 
culi light with the rays of H, K and Ca*. This is confirmed by the fact established by the author [7], that the 
passage of nearly every flocculus (or facula) through the visible center of solar disc is accompanied for a certain 


period of time by geomagnetic disturbances, Consideration of these matters is contained in author's other articles 
[8, 9] and we will not stop to review them here, 


rere oe i ; 
In fact it is possible that the phenomenon is even more complicated; it is possible that it encompasses the 
Processes involving impact waves and others. 
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We will pass now to the main question for this article, which is the consideration of the nature of differ- 
ences between the disturbances of sudden and nonsudden commencements. It is obvious that these differences 
are related to the differences in some properties of the streams producing the two types of disturbances. Allen 
[10] believed that these differences are conditioned by the variations in the character of entry of solar corpus- 
cles into the upper layers of the earth's atmosphere; this is to say that in the case of disturbances with a sudden 
commencement, we are dealing with a relatively short period of emission of the corpuscles (lasting only as long 
as the formerly indicated nonstationary process), after which the corpuscles advance in a broad front of a defi- 
nite breadth (Figure 1, a). On the other hand, in the case of disturbances of nonsudden commencement, in which 
the regions of emission of the corpuscles are relatively stationary and last for months, the intrusion of the corpuscles 
into the earth's atmosphere takes place because the stream of the corpuscles is transposed in space in the direc- 
tion of solar rotation, and after every 27 days “overtakes” the earth. Allen believes that in the first case (Fig- 
ure 1, a) the foremost border (forward front) of a stream is sharply defined, and that this is the cause for the 
sudden commencement in corresponding disturbances. In the second case (Figure 1, b) the earth falls into the 


Te " s . . . 
side" portions of a stream where the border is less sharply defined; this fact leads to disturbances of nonsudden 
commencement. 


Fig. 1. Differences in the geometrical character of the streams producing 
geomagnetic disturbances of sudden and nonsudden commencement. 


It is difficult however, to agree with Allen's view. In the first place, it follows from the very character 
of geomagnetic disturbances, whether their commencements are sudden or nonsudden, that the stream consists, 
in both cases, of separate condensations divided from one another by the zones of considerably lesser conden- 
sation of matter. Consequently, there is no reason for the existence of any definite border at all in the first 
case. 


The second consideration rests on the previously mentioned fact that the two types of disturbances differ 
not only in their early and final phases, but also in the very character of their progress. Moreover, it follows 
from Figures 1a, and 1b that after the entry of the earth into a stream, its situation appears to be identical in 
both cases, because the following (after the commencement of the disturbance) progress of events is determined 
by the velocity of flight of the particles from the sun (velocity of atom movements in Figure 1, b along the 
tangent to the earth's orbit is insignificantly small in comparison with the velocity of the particles moving away 


from the sun). 


The author of this article believes that the differences between the two investigated types of disturbances 
are conditioned by different physical characteristics of the corresponding streams. As such characteristics we must 


list primarily: 
a) velocity _v of the condensations and of the separate atoms in the streams, 
b) concentration n of atoms in a stream and in the condensations, 
c) magnetic field of the condensations in a stream, 


These three basic characteristics will be considered below. 


First let us pass to the question of the velocities of condensations in the streams, A number of data indi- 
cate that velocity v in the streams producing disturbances with sudden commencements is, on the bictete 
greater (in many cases considerably greater) than velocities v in the streams producing disturbances with the 


nonsudden commencements, 


As the disturbances with sudden commencement, we will investigate those disturbances which e related 
to the chromospheric flares. By basing our judgment on the retardation of geomagnetic disturbances Hie respect 
to the appearance of the flare) we can estimate velocity Vv fairly accurately. On the average, this period of 
retardation is equal to 20 hours [11], which gives us the velocity of about 2000 km/sec. (although greater veloc- 
ities can also be observed), On the other hand, velocity of the ‘corpuscles in the streams producing disturbances 
with a nonsudden commencement is undoubtedly lower. In the same way, a study of disturbances produced by 
the passage of the flocculi over the center of the solar disc gives us the period of retardation equal to a few days, 
and sometimes up to ten days. * 


Similar results, pertaining to the differences in the velocities of condensations in the streams of the two 
types, are obtained from a study of the shift of line Hq in the spectrum of the Aurora Borealis (see [12], p. 412). 


At the same time the difference between the velocities of atoms and condensations in the streams is of a 
great significance in determining the character of geomagnetic disturbances. If we should inspect a simple model 
depicting the first phase of a geomagnetic disturbance, a model in which kinetic energy in the corpuscles of a 
stream is expended on changing the earth's magnetic field, we will be able to represent the change in the inten- 
sity of this field with respect to the time, as [13]: 


di mF nv 


a 1 
dt 2H a8 (1) 


Here m is the mass of atoms constituting the stream; F is the effective area of the magnetic "diameter" of the 
earth; a is the earth's radius; n is the concentration of atoms in a stream at the distance of the earth from the 
- sun; Hp is the horizontal component of the earth's magnetic field. 


We can see from Equation (1) that the rate of change of the earth's magnetic field during the starting phase 
of a magnetic disturbance depends greatly on the velocity of the corpuscles in a stream, It is especially true, on 
the basis of the velocity data introduced above, that we can expect particularly fast variations in the earth's 
magnetic field to be observed at the intrusion of the corpuscles producing disturbances with a sudden commence- 
ment. From the same data on velocities and from Formula (1) it becomes obvious that among the weak distur-~ 
bances, those of sudden commencement appear to be most intense, as was noted at the start of this article. 


The significance of the atom concentration in the stream is the same as that of the velocity. Here however, 
into the various factors which are used in determining the magnitudes of v and n enters the product of these two 
magnitudes (or to be more exact, the product of v and n with various exponents). For this reason, up to this 
time we have not obtained any definite data pertaining to the differences in the magnitudes of n in the distur- 


bances with a sudden and nonsudden commencement, Only the order of magnitudes n is known for the given 
velocities v (see [19]). 


Let us pass to the third factor — the magnetic field of those condensations which compose the streams. 


A great many factors indicate that such magnetic fields do exist in the streams moving away from the sun; 
the knowledge of this fact is specifically used in explaining some variations in the cosmic rays. Our basic prob- 
lem here lies in determining which of the condensations have the stronger field — the ones which produce the 
disturbances with a sudden or those which produce the disturbances with a nonsudden commencement, It can be 
said, with a fair degree of certainty, that a substantially more intense field may be expected to exist in a case 
of those condensations which produce the disturbances with a sudden commencement. This can be deduced from 
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the fact that these disturbances are related to the spots characterized by strong magnetic fields. Moreover, we 
have previously stated that there exists a basis for the belief that, in considering the emissions of the condensa- 
tions producing such disturbances, we are dealing with a nonstationary process of electromagnetic nature. For 
this reason, we can expect that at the emission of the condensations from a given solar region, these condensa- 
tions would carry with them noticeable magnetic fields "frozen" into the matter of the condensations. The part 
played by these "frozen" magnetic fields at the approach of the condensations to the earth should be considered 
as an auxiliary factor altering the earth's field at the time of a disturbance. Apparently, these magnetic fields 
may either intensify or weaken the variations of hydrogen, which are determined by intrusion of a condensation 
into the outer layers of the atmosphere. This problem calls for a special investigation. 


As to the condensations which produce disturbances with a nonstationary commencement, here the "frozen" 
magnetic fields must be less intense. This is proven by the fact that these disturbances are not related to the 
spots and, furthermore, they are observed at the time when sunspots are absent. Also, the disturbances with non- 
sudden commencements are not related to the intensifications of the solar radio emission, in which an important 
part is undoubtedly played by electromagnetic effects, etc. It follows from what has been said that the differ- 
ences in a specific "frozen" field of the condensations may constitute a very important factor in the problem 
of explaining the differences in the character of the two types of the disturbances. 


In conclusion we must note the following conditions. Even though division of the disturbances into those 
of sudden and those of nonsudden commencement appears, on the basis of the whole immense scope of material 
collected by now, to be necessary and well grounded, nevertheless such a division is difficult to carry out in 
many instances, Thus, for instance, in some cases the consecutive disturbances are characterized by sudden 
commencement, etc. These cases, however, represent exceptions to the general rule. Furthermore, these anom- 
alous factors can be always explained if we take into consideration the fact that, in accordance with all that has 
been said in this article, the character of the disturbances depends on a great many conditions, For this reason 
we will refrain from conducting a special investigation of this problem. 


_ Crimean Astrophysical Observatory, 
Academy of Sciences, USSR Received November 1, 1956 
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A POSSIBLE NEW TYPE OF FLUORESCENCE OF THE EARTH'S 
ATMOSPHERE 
I. 8. Shklovskii 


As a result of cascade transitions during the resonance absorption of the solar 
emission line Lg, by oxygen atoms of the earth's atmosphere (transition °P2 — 3D) 
infrared lines AA 11299 - 11287 and A 8446 will be emitted. This should lead to 
a twilight flash of these lines. The expected intensity of the infrared lines of oxygen 
during the twilight flash is calculated. For an undisturbed sun it is of the order of 
108 photon/cm? sec. The line » 8446 is in the region of the band (6, 2) of the vibra- 
tion-rotation spectrum of hydroxyl in the earth's atmosphere. This makes the obser- 
vations of the twilight flash of the oxygen line difficult but not impossible. The 
author emphasizes that systematic observations of this flash give a new possibility 
of investigating the emission lines of the Lyman series in the solar spectrum as a 
function of solar activity. 


The distinct flash of the yellow sodium doublet in twilight is well known. We are also familar with the 
twilight flash of the red oxygen doublet AA 6300-6364A and of the A 3914A Ny” band; it is also possible that 
the oxygen band at A 8600A is also strengthened in twilight. 


The lines and bands which have been mentioned are enhanced in twilight by resonance fluorescence. In 
the present article we shall discuss the possibility of twilight flash of certain infrared lines under the influence 
of ultraviolet radiation from the sun. This concerns the possibility of ordinary, rather than resonant, fluores- 
cence of the earth's atmosphere. ! 


It has been known for a long time that the ultraviolet spectrum of atomic oxygen contains a strong line 
whose wavelength is remarkably close to the wavelength of the familar hydrogen line Lg. For example, Hop- 
field [1] in a laboratory study of the OI spectrum pointed out the very intense line A 1025.84A which is close 
toLg at 1025.73A. The intensity of this oxygen line was estimated at 9 points on a ten-point scale. This line 
originates in the transition from the “D321 level to the ground level of oxygen, 3p. Later data [2] give the 
exact wave number of the 3P, ~ *Dp transition as Vo = 97488.14 cm”, whereas for Lg we have Vy = 97492.20 


cm7*, The wave number difference between these two lines is only ~ 4 cm", while the wavelength difference 
isAA © 0.04A 


Because of the extremely close wave numbers of the hydrogen and oxygen lines and the fact that the °P, 
level of oxygen is lowest we can expect resonance effects in the earth's atmosphere due to solar emission in the 
Lg region. Let us consider this in greater detail. After resonance absorption of an Lg photon an atom of atmos- 
pheric oxygen rises to the p§21 state. It can then drop into the ground state 3Pz or into 3P2,,9. In the first 
case, because of the large absorption coefficient at the center of the resonance line the emitted photon is immed-~ 
iately absorbed by an atom of atmospheric oxygen. This occurs in planetary nebulae. Therefore after Lg absorp- 
tion the ground state *P, of an oxygen atom will be reached preferentially by the cascade 3p? —> Pohgeet, So 
—> 5p. Inthe °p}.o44 > PR transition the infrared oxygen lines XA 11294 and 11287A will be emitted. The 
transition °P9,149 —> oS gives the well known infrared line A 8446 A which is characteristic of the auroral spectrum, 
The transition from *S; to the ground level =P, gives the ultraviolet line in which we are interested. 


As a result of the fluorescence process which has been described every Lg photon that is absorbed by 
atmospheric oxygen yields the two infrared photons: AA 11294-87A and X 8446A. Thus we can expect from 
theory that the upper layers of the earth's atmosphere will 


fluoresce, emitting infrared oxygen lines, under solar 
irradiation. 


We shall now discuss the intensity of this fluorescence and discuss the Possibility of observing it. Rockets 
have been used to investigate solar ultraviolet radiation, especially Ly, at high altitudes. According to [3], the 
La flux from the undisturbed sun is 0.10 erg /cm*-sec and the line width is less than 1A. The Lg line in the 
solar spectrum was recently observed [4] by rockets at 130-150 km, although no direct intensity data were given 
in this report. According to Goldberg's recent theory of the production of the La solar emission line [5] the in- 
tensity at the center of the line is determined by the Planck function for the kinetic temperature of the upper 
chromosphere layers where the optical density at the center is unity. This temperature as estimated from the 
Lq intensity [5] is 6000°. This is approximately the temperature which should exist in the region where Lg is 
emitted. The chromospheric layers where the Lyman lines originate is at about 4000 km. From this it can be 
calculated that the Lg flux should be about one order of magnitude smaller than for Lg. It is possible, however, 
that Goldberg's theory neglects some factors, so that the Lg flux may be larger. But for our further calculations 
we shall assume the Lg flux in solar radiation at the edge of the earth's atmosphere is FLg = 0.01 erg /cm?-sec = 
= 5-10® photon/cm? sec, regarding this as a minimum. 


We write the selective absorption coefficient of oxygen in the usual form 


ee / , 


me 4m (v -~ Vo)? -+ (y/41)?’ 


_ Bre? 


ae alia 


= 2.22-10° sec 1, from which, for yv—y>y7/4n, 


f OG. F 
Sa 1258-108) 


(ve) 


(2) 


Considering the character of the 3p, —*D,,1 transition, the oscillator strength of the oxygen line » = 1025.83A 
must be quite large. According to the most recent data, in a column of unit cross section at atmospheric alti- 

ao 

( 7 -2 = 
tude ho > 120 km the number of oxygen atoms is given by \ (Olfidkse Nig 0740? oma~, Atte a0 

No 
150 and 175 km, respectively, Np, = 3- 107 cm7”, 10 cm7? and 2.5-10* cm7*, When v — v9 = 1.42. 10" sec™4 

f : ; = wien es a 
(which corresponds to AA = 0.05A) according to Equation (2), Sy /f{ = 0.16- 10522 cms; for AX = 0,10/Ay 
s,, /f = 11.89.1072" cm?, and for A A = 0.15A, S,,/f = 8.4-107!8 cm*. We can therefore conclude that at 175 km 
Vv 2 . ’ . , ; 

the absorption of solar Lg by atmospheric oxygen cuts out a band of width ~ 0.1A (iff~ 1) inthe Lg line; at 
150 km the width is~ 0.2 A; at 130 km it is~ 0.3A and at 120 kmitis ~0.5A. The width of Lg must be 
less than 0.8 A, since the width of Ly is less than 1A. Therefore at hg = 175 km we have absorption (or more 
precisely, scattering) of more than 1/8 of the flux of solarLg, at hg = 150 km more than 1/4, and at 120 km 
more than 1/2. When f ~ 0.01 the width of this band is ten times less and the relative number of scattered photons 


is diminished in the same ratio. Hereinafter we shall assume f~ 1. 


Thus in a column of unit cross section in twilight every second there will Be produced from 7-107 = 
3-10° infrared photons at AA 8446 and 11294-87 A for 175 > No > 120 km. The intensity pe oxygen ne 
d 6363 A in twilight flash is ~ 3-10 photon /cm? sec. We can thus conclude that ie fheoperieal oie fo) 
infrared oxygen lines in twilight flash is not large but still sufficient to raise the eiot of wat ae ae 
Such observations, which can be very important (see below), encounter sev erabeleot Sunn GS We shall, o 
A 8446 A, because modern radiation receivers (such as electron- 


i ider the shorter wavelength | 
EUR eee as ; at A 8500.A. In the region of the oxygen 


j i itivi t A 11300A than 
optical image converters) have much less sensitivity a th . ne 
a r ane we find the intensity (6, 2) rotation-vibration band of atmospheric hydroxyl. The emission of the 


night sky in this spectrum region can be described by Table 1 [6]. 
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TABLE 1 The characteristics of the suggested oxygen line 
d 8446A in the twilight flash spectrum are shown in paren- 
theses. Its intensity is based on the assumption that the 
Identification intensity of A 6300A in the night sky spectrum is 10° photon 


AA | I (6308 = 8) 
per cm?.sec. Intensities 5-20 are assumed to cover the anti- 
cipated range for different solar angles. We can thus expect 
pas sn tea is that in twilight flash the intensity of the oxygen line will 
8446.0 (5—20) 01 be about the same as in weak satellites of P branch lines in 
iy ae oo P the (6, 2) hydroxyl band. But we must keep in mind that in 
8497.3 7 (6,2) Py, twilight the exposure time is limited and the brightness of 


the sky is quite high. All of these circumstances interfere 

seriously with observation of the twilight flash of \ 8446A 

but do not make it impossible. The sought line is 14A 

from the strong P, component and 11A from the weak P2 
component of the (6, 2) band. Therefore successful observation would be assured by a spectrograph of sufficiently 
high dispersion and resolving power. The sensitivity of the receiver must also be very high. Such spectrographs 
are known to exist at the present time. 


The foregoing calculations for the intensity of infrared oxygen lines in twilight flash refer to an undisturbed 
sun. During large flares and perhaps other manifestations of solar activity the Lg flux may be increased by a 
factor of a few tens. On the other hand, the fluctuations of the Lg and Lq solar lines must be synchronous. 
Therefore, the intensity of the oxygen lines of the earth's atmosphere in twilight flash would be a very sensitive, 
practically ideal indicator of the solar ultraviolet flux in Lyman lines. 


There are strong reasons for believing that Lyman radiation is of great geophysical importance. For 
example, it is possible that one cause of the formation of the D layer is ionization of NO molecules by Ly. 
Therefore during the course of the International Geophysical Year it is proposed to make a special systematic 
search by sending up rockets with spectrographs during flares. If it becomes possible to detect twilight flash 
of oxygen lines we shall have a very efficient as well as a relatively inexpensive and simple method of study- 
ing the Lyman lines in solar radiation. 
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OBSERVATIONS OF HYDROGEN EMISSION IN AURORAE 


Iu. I, Gal'perin 


_ Data are given on 17 spectrograms containing Hg, taken during 
observations at the magnetic zenith and the magnetic horizon. Profiles 
of Hg have been determined from the three best spectrograms taken at 
the zenith. These give a maximum value of the velocity of approach 
of the order of 2000 km/sec. The velocity corresponding to maximum 
density, equals 350 + 100 km. 


Investigations over many years by Vegard [1] and Gartlein [2], who used low-dispersion spectrographs, 
have proved that aurorae contain broadened Balmer lines of hydrogen. Spectrograms taken by Meinel [3] and 
Gartlein [4] at the magnetic zenith and horizon, with high dispersion and high resolving power, have provided 
convincing evidence that protons move approximately along the magnetic lines of force of the earth's field with 
velocities of the order of thousands of kilometers per second. This has also been confirmed by Vegard [1]. 


However, to our knowledge, besides these spectrograms which took the magnetic field into account, the 
literature contains no other similar hydrogen line profiles with sufficient resolution. 


During the period from December 1955 to April 1956 at the northern station of the Institute’ of Atmospheric 
Physics of the Academy of Sciences USSR, which is located in the auroral zone (y = 68° 38's A = gh13M_\3 )y 
observations were made for the purpose of preparing to obtain profiles of hydrogen lines in the aurora during the 
International Geophysical Year. 

This work was done with an SP-48 GOI spectrograph,* of 1:0.8 aperture ratio and 83.5A/mm dispersion 
for Hy and resolution limit about 2A, using pre-exposed Panchrome X film [6]. During this period 17 spectra 
were obtained containing Ha, of which two were taken at the magnetic horizon and the others at the magnetic 
zenith (more precisely, 8° north of the magnetic zenith), The magnetic horizon photographs were taken without 
reinforcing optics (the collimator aperture ratio was 1:4.7). In the magnetic zenith direction an arc of 4° was 


projected on the spectrograph slit. ( 
The Hq line was obtained from a quiet HA arc, from diffuse DS and PS emission and from coronal forms 
C and D (this notation for the forms is taken from [7] ). 
All of the hydrogen lines belonged to weak and average aurorae. Table 1 gives the K index of magnetic 
disturbances during exposure averaged for all observatories of the Soviet Union; the K index is seen not to have 
exceeded 5. 


The Hy intensity varies strongly compared with the intensity of the first positive Nj system. In some 
spectrograms (Nos. 1, 2, 14, 15) Hg is pronounced but the first positive nitrogen system is almost absent. There 
are some instances of the reverse. An absolute calibration was not performed, but it can be stated that the Hq 


intensity did not once exceed the intensity of [OI] A 6363.8 A. 


* This spectrograph is of the type described in [5]. 
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TABLE 1 


Geomagnetic 
coordinates on 


Exposure ‘ 
Auroral form , celestial sphere 
No. uate (Greenwich time) | K index P 
|. | ee ee 
Dec, 4955 
4 17 46630™—20530™ | HA 2 Magnetic 
2 47—18 21 30 —01 30 | DS. HA 2 Harigon 
Feb. 4956 
3 3S 20 00 —22 00 PA. G. DS, PS 4 
4 5 20 00 —23 10 HB. HA. PS 3 
5 41 46 45 —19 10 HB. HA. RA bs) 
6 41 19 10 —19 30 HA. RA. G 5 
7 41—12 22 20 —02 30 RA. C. PS 9) 
8 42—13 47 15 —00 15 HA. PS 3 
9 13 00 15 —02 00 HA. C. PS Ye 
40 16 46 45 —18 15 HA 4 
11 16 18 15 —20 30 HB. RB. PS 3 Magnetic 
12 27 20 00 —23 30 HA 5 wanith 
Mar. 1956 
13 6 47° 30: —19) 30 HBS Cabs 4 
14 41 49 15 —21 40 HA os 
15 SA 21 40 —22 40 CD 3 
46 41 —12 22 40 —02 30 Cc. PS oS 
17 16—17 24 00 —00 40 HB. RB 2 
AeA BL Eee 


Velocity in km/sec 


A comparison with solar data [8] shows that in all cases about one day before exposure a calcium flocculus 
was located within the interval 8-15° from the point of projection of the earth on the solar disk. 


Some spectrograms, particularly Nos. 1 and 2, were distorted by nonuniform pre-exposure. The most 
reliable Hg profiles at the magnetic zenith were obtained from three well developed spectra (Nos. 14, 15, 16) 
with the highest ratio of Ha intensity to the intensity of the interfering bands of the first positive Nz system. 


Figure 1 shows traces of spectra Nos, 14, 15, 16 obtained with an MF-4 microphotometer ( 21 x magnifi- 
cation, microphotometer slit width 0.65 mm, resolution limit r 3A). For comparison Figure 2 shows Meinel's 
traces (resolution limit 7A) and Gartlein's traces (resolution limit about 15 A). 


In order to obtain profiles the trace of the faint background of the first positive Np system was subtracted 
from the Hg trace. Three different background traces were subtracted and the results were averaged. Figure 3 
shows the profiles Nos. 14, 15 and 16 obtained in this manner. Densities were not converted to intensities. These 
profiles permit us to determine the velocities of protons (VA, VB, Vc, Vp) corresponding to the following points of 
the profile: A is the maximum observed velocity of recession from the observer ("red shift"); B is the maximum 
density; C is the density 0.1 of maximum; D the maximum observed velocity of approach to the observer 
("violet shift”). These velocities are given in Table 2. 
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Fig. 1. Traces Nos. 14, 15 and 16 and background without Hq. 
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Fig. 2. Traces obtained by Gartlein (No. 1) and Meinel (No. De 
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For spectrograms numbers 1 and 2 obtained at the magnetic horizon 
we excluded as far as possible the effect of nonuniform exposure and also, 
as described above, the influence of the faint background from the first 
Fig. 3. Hg profiles inrelative positive Ny system. ‘The Hq profile at the magnetic horizon is broadened 
densities. Magnetic zenith. but its center is not shifted. The half width of the profile in spectrograms 


Nos. 1 and 2 is 14 and 17A, respectively, which corresponds to velocities 
+ 300 and + 400 km/sec along the line of sight. When the effect of the 
a” apparatus on the profile is taken into account we obtain 11 and 14A, 
NZ 


The error in determining velocity vp is 4 100 km/sec; for va, V¢ 
and vp the error is + 200 km/sec. The result obtained for the maximum 
NI observed proton velocity of approach v,;, depends on the film sensitivity 
and background fluctuations; therefore the values obtained for vp must 
not be regarded as exact. In this respect vc was more reliable. When we 
take into account the influence of the apparatus on the profile the value 
for va is somewhat reduced; therefore the reality of the "red shift" re- 
ee mains doubtful. * 
We note only small deviations with regard to vg. The values obtained 
for Vc and vp were smaller than those obtained by Meinel (K index 9) but 
larger than Gartlein’s (K index 5). 
NIG 


+ 500 0-500-000-1800 2006 km/sec 


corresponding to + 250 and + 300 km/sec. As a result of nonuniform 
exposure the accuracy of the profile obtained from these spectrograms is 
less than for Nos. 14, 15 and 16. 


The author is deeply grateful to V. I. Krasovskii for his continued 
interest and direction of this work andto M. A. Ermolaev and V. F. 
Musatova of the staff of the northern station for their assistance with 
Fig. 4. Hq profiles inrelative the experimental work. 
densities. Magnetic horizon. 
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*As was first noted by I. S. Shlovskii [9], in high Hq emission the mean free path is much smaller than the 
radius of curvature of the spiral, and the time during which the incoming particle is charged is much shorter 
than the time during which it is neutral, so that the particle trajectory will be a broken line rather than a spiral. 


It appears to us that the "red shift" may be simply the result of particle deflection during multiple elastic 
scattering. 
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BRIEF COMMUNICATIONS 


ON THE UTILIZATION OF A THIN LUMMER-GEHRKE PLATE 


Lee rer din 


A thin Lummer — Gehrke plate is described. The formula for the resolving power 
is derived and calculated parameters of the plate are given. These show that the con- 
clusions reached by A. V. Merkulov in his paper "A Stellar Spectrograph with a Lummer — 
Gehrke Plate” are erroneous. 


Characteristics of a thin Lummer — Gehrke plate and considerations showing the 
inadvisability of its utilization are given. 


A.V. Merkulov [1] has discussed a thin Lummer-Gehrke plate, of the order of 0.2 mm, placed in optical 
contact with a thick glass base. Table 1 shows data obtained by Merkulov and by the present author for the Hy 
line. 


TABLE 1 


Data of present author 


Material of plate Flint L-23 Flint TF4 
fig = 1.78243 
Material of base plate Crown L-20 Crown K3 
Initial Ne = 1.50760 
data Thickness t of plate 0.2 mm 0.2 mm 
Length L of plate 25.4 mm 30 mm 
Angular dispersion oe 0.004 rad/A 0.016 rad/A 


dA 


Spectrum order m 


"Plate constant” g 11.48A 
Calcu- Number of interfering beams N 43 
lations Resolving power R 2.43 -104 


Small differences in the results can be attributed to the different kinds of glass employed. However, the 


difference in the resolving power by a factor of more than a thousand results from Merkulov's error in deriving 
the formula for the resolving power (see [1], p. 151, Formula (16) ). 


A correct derivation can be obtained by appropriately modifying the derivations given by Frish [2] or 
Tolanskii [3]. 


Let ny be the refractive index of the ‘plate and ng the refractive index of the base, with My > Ng. Denoting 
the thickness of the plate by t (Fig. 1), we obtain for the difference between two successive beams 
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Fig. cli 


2 tny 
cosr 


A == n,(AB + BC)—n,AD = 


— 2in,tan rsini. 


Using ngsin i = nysinr, we obtain 


N=2in, (= —tanr sin r) == 2 tht; COST’: 


2 
. ° 1S FRR n . . a 
By substituting cosr = YV1—sin?r = Ve - sin?z,' we obtain 


1% 


A = 21 Vr? — ri sin. 


The condition for the formation of a principal maximum is 


md = 2t V rn? — n2 sini, 


uy ; 
or, replacing the angle i by the Bragg angle € = Stee 


md = 2tV rn? — ni? cos*e. (1) 


This expression is satisfied by the angle € at which an m-th order principal maximum is formed; the 
angle € + Ae corresponding to the (m + 1)-th order for the same wavelength satisfies the condition 


(m+ 1)X\ = 2t Vr — nz cos? (< + Ae). (1a) 


The order is determined from Equation (1): 
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re 
m=! Vni— ny. (2) 


Squaring (1) and (1a), subtracting term by term, and assuming € to be small, we obtain the angular separ- 


ation of two successive orders: 


pa Ve —n} (3) 


Differentiating (1), we obtain the angular dispersion: 


dn dn 

2 3.5 fd Vode, sage 
dene bins (my Te aes =) (4) 
dx nam Ange : 

Hence 
AnZe 

OI et eee es 

n?—n2—r(n en any (9) 

irae : ay edn 


Replacing de by the angular separation Ae of two successive orders according to (3), we obtain the separa- 
tion of successive maxima expressed in wavelengths (the so-called "plate constant”): 


x2 Vr = n} 
dn dn 
2 2 ea ae Apel) 
at|n?— nk —r(m Dr Ng 2) 


f 
} 


(6) 


The wavelength difference 6 A between two lines that can barely be separated with N interfering beams is 


cena 
Cs ae a SVG. O15 (7) 
; rae dn dno 
2 1 
2tN [— ng— (ns rn Ng a )| 


For plate length L the number of interfering beams is 


L 15 
N ne 
AC  2t-tanr 


Since the angle r is very close to the angle of total internal reflection, we can assume approximately 


ny 


7 in =Landtan r = ——- 


hence 
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LY et 
Wis we) (8) 
2t ri 


n 
EG 4 UE SOON Sie er etc We =) 
Ve 1 [nt "a a(n ie 


Vn? — 2 


(9) 


F ‘ A c 
The solution of Hartmann's dispersion formula ” =n, + J23;3.~= «Which is required for our calculation, gives 


. do 
for K3 glass the following constants: c = 76.6615, no = 1.49218, Xo = 1595.74A, and differentiation of the formula 
gives the dispersion for the red hydrogen line A, = 6562.8 A: 


ae 284073 40 


dx 
; ; és u Ye 1 aan =6 
For TF4 glass we obtain the following values: c = 162.317, no = 1.69382, Apo = 2359.86A, aN = — 9.1918 - 107°, 


Thus such a thin Lummer-Gehrke plate loses its principal advantage, that of high resolving power, because prism 
spectrographs possess resolution of the order of 2- 104. 


The dispersion of this plate is almost indistinguishable from the dispersion of a prism spectrograph, with 
which Merkulov compared the plate. The plate retains one advantage over a prism spectrograph, its higher trans- 
mission. 


The plate has another shortcoming by comparison with a prism spectrograph: the interval of the spectrum 
which can be studied does not exceed 10-15A. 


Moreover, it is very difficult to prepare such a thin plate with excellent surfaces (to 1/80 of a wavelength, 
as shown by Tolanskii [4]) in optical contact with a base plate. It must also be kept in mind that a plate is very 
sensitive to temperature fluctuations and requires thermostatic control to within a few hundredths of a degree 


(Korolev [5] ). 
It is thus clear that Merkulov was in error in asserting the advantage of a thin plate, and the advisability 
of using it is at least doubtful. 


Received October 15, 1955 
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A NUMERICAL METHOD FOR DETERMINING THE STELLAR DENSITY 


IN 


A SPHEROIDAL CLUSTER 


A. S. Sharov 


P, N. Kholopov recently worked out and successfully applied a numerical method for determining the 
spatial density of stars in a spheroidal cluster {1]. In this method it is assumed that the cluster consists of n 
spheroidal shells of equal density, formed by concentric, similar, and similarly arranged oblate spheroids, with 
the ratio of the major to the minor semi-axis represented by 7. The principle of this method, first developed 


by Wallenquist [2] for spheres, is as follows. We calculate the volume 
of that part of the outermost spheroid which in projection on the celes- 
tial sphere gives the first outer elliptical ring. By counting the number 
of stars in the ring, we find the mean density in the layer. Then we 
calculate the volumes of those parts of the second and first layers which 
give in projection the second more inward ring. By counting the stars 
in the second ring and knowing their density in the first shell and also 
the volume of the region in which they are found, we obtain the density 
in the second shell. Proceeding to the center, we can determine the 
stellar densities in all layers. The essence of this method is the calcu- 
lation of the volumes of the corresponding parts of the spheroidal shells. 
For the calculation of these volumes it is possible to suggest a method 
based on the theory of affine transformations. This method is much 
simpler than the one applied by Kholopov. 


Let us consider first a spherical cluster in which spheres with radii 
Am-1» 4m and aj-y, aj form any two shells. We describe around the 
spheres with radii ay,-, and ay, two right cylinders with generating 
lines parallel to the line of sight. In Figure 1 a cross section of such a 
cluster is given, That part of the volume of the shell (ay)-1, a) lying 


between cylinder surfaces can be calculated as the difference between the sphere volume (radius a,,), and the 
cylinder volume (radius ayj-1) and two polar caps. Well known formulas of elementary geometry give 
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a 


Sy aas een ae 
a3, — 2na? | Va, Oe 


oo] BS 


£5) (en VERA) = Sete, — a 


The volume formed by the same cylinders in the shell (a;-4, aj ) is obviously: 


{( 


ey | 5) 
a? Gre) 2 


(at — 08)" — [(at_, — a8,_,— (at_, — a2). 


We now consider a spheroidal cluster and apply the theory of affine transformations. As we know, a 
spheroid is an affine transformation of a sphere and can be obtained from it by contraction along the coordinate 
axes and a rotation, Thus all volumes change proportionally to the product of the contraction coefficients. 

Let the spherical cluster contract along the line of sight with a coefficient 7. Then the volumes V and W change 
by a factor n. Thus, the formulas for the volumes of interest in the case of a cluster contracted along the line 
of sight differ from the preceding formulas only by the factor n. Here a; is the major semi-axis of the con- 
tracted cluster. The formulas do not change further if the minor semi-axis of the cluster makes some angle i. 
with the line of sight. In fact, all volumes change in an affine contraction by a coefficient 7, independently 

of the angle formed by the cylinder axis with the line of sight. We select this angle in advance in such a way 
that by contraction of the cluster along the ultimate minor axis (which makes an angle i with the line of sight) 
the cylinder axis will fall along the line of sight. In this way, the formulas for the considered volumes are 
independent of the angle i and the geometric picture created a ey this transformation will correspond to the real 
conditions of observation of a spheroidal cluster. 


P. K, Shternberg State Astronomical Institute Received October 15, 1956 
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CHRONICLE 


SEVENTH CONFERENCE ON METEORITES 


From the 14th to the 17th of November, the Seventh Conference on Meteorites was held in the House of 
Scientists of the Academy of Sciences of the USSR, convened by the Committee on Meteorites of the Academy. 
The presidency of the Academy of Sciences had invited to the conference foreign scientists, who took part in 


its work, 


The Conference was opened with a report by V.G. Fesenkov on “Progress of the study on meteorites.” The 
author cited the most outstanding works in the field of meteorites during the interval after the sixth Conference, 
and described the great role of this study in the solution of a number of important cosmogenic problems. 


The report by A. A. Iavnel' on “Regularities in the composition of meteoritic material and the conditions 
of formation of meteorites” aroused great interest among the members and gave rise to a lively discussion. 
Thorough analysis of data principally gathered from chemical analysis of meteoric material, including mineralogi- 
cal and petrographic characteristics of meteorites, allowed the author to find features common to all varieties 
of meteoritic material. The division of meteorites according to these features into 5 groups and 6 sub-classes, 
sufficiently different from each other, made it possible for the author to draw a number of far-reaching conclu- 
sions, particularly concerning the origin of the whole complex of meteoritic material coming from five asteroids 
having a layered structure. The regular variations in the chemical composition of meteorites, and consequently 
of the asteroidic bodies, leads to the concept of a common origin of these asteroids (and possibly comets) in an 
original proto-planetic cloud, ; 


V.G, Fesenkov touched, in his report on "Some general problems of meteorics and astronomy,” upon a 
number of problems connected with the origin of lunar craters, asteroids, and meteorites in the light of contempo- 
rary concepts of the origin of the solar system. Many data of investigation of lunar forms are more in favor of 
the volcanic than of the meteoric hypothesis of their origin. In any case, there remains the given problem, 
which is still open to discussion. Investigations of meteorites from the pgint of view of their luminosity, color, 
period of revolution, and forms make it possible to conclude that the asteroids are fragments of greater bodies 
and cannot have originated from the destruction of only one parent planet, which is confirmed by the different 
periods of revolution of the asteroids and by the division of meteorites into certain groups according to their 
chemical composition (Urey and Craig, Iavnel'). The origin of asteroids and consequently of meteorites is con- 
nected with the conceptions of the initial material from which our solar system originated. In this connection, 
the problem of the origin of the chondrulic structure of stony meteorites is a basic one, in the consideration of 
which the author came to the conclusion that the most probable hypothesis for the origin of the chondrulic struc- 
ture is "that a similar structure was already inherent in the planet itself and was inherited by it from the original 
gas-dust nebula formed simultaneously with the material of the planet-asteroid itself.” 


Later on, Prof. N. Bonev (Bulgarian National Republic) in his lecture on "A meteoric hypothesis of the origin 
of lunar craters," starting from definite premises on the formation of craters as a result of the fall of meteorites, 
proved mathematically by a statistical calculation of the distribution of lunar forms between its eastern and 
western hemispheres that the craters are not related to a meteoric origin, at any rate not substantially. 


The whole morning session of the second day of the Conference was devoted to papers on the Sikhote- 
Alin* iron meteorite.A paper was read by N. B. Divar’ on "Final elements of the atmospheric trajectories of the 
Sikhote-Alin’ meteoritic rain,” in which more precise data were given on the trajectories of the meteorite in 
the terrestrial atmosphere. In his report, "The pattern of the fall of the Sikhote-Alin’ meteoritic rains” Ess 
Krinov gave an exhaustive picture of the fall of the meteorite on earth, obtained by an analysis of the distribu- 
tion of the meteoritic masses and the character of the craters formed. From a detailed review of the characteris- 


tics of the ellipse of scattering of the meteoritic rain, a picture was established of the motion and disintegration 
of the meteoritic body in the atmosphere. 


* Transliteration of Russian — Publisher's note, 
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The paper on “Chemical composition of the Sikhote<Alin’ meteorite” read by M. I. D'iakonova was de- 
voted to the results of chemical analysis of the meteorite. Because of an apparent homogeneity in the compo- 
sition, it was possible to determine the total composition of all the meteoritic material and also the average 
chemical composition of nickel-iron of some samples of the meteorite. The total nickel content in the meteorite 
was shown to be 5.94% and in nickel-iron 6.00% (by weight). In the second part of the paper, data were given 
of chemical analyses of some of the most abundant minerals in the meteorite: troilite FeS, schreibersite (Fe, Ni)sP, 


chromite (Fe, Mg)Cr,04, kamacite (w-iron), together with a short description of the methods used for their de- 
termination. 


L. G. Kvashi in his report on "Mineral composition and structure of the Sikhote-Alin' iron meteorite" 
gave detailed information on the petrographic characteristics and the mineral composition of the meteorite. 


Starting from the chemical composition of separate minerals, based on their quantitative mineralogical 
analysis, the total composition of the meteorite was calculated and was shown to be approximately the same 
as the average composition determined by chemical analysis. The mineral composition of the meteorite con- 
sisted principally of kamacite (a-iron), taenite (y ~iron) and the additional minerals troilite, schreibersite, 
chromite, etc. From a detailed review of the structural characteristics, it was established that the Sikhote- 
Alin’ meteorite belongs to a group of very coarse~structured octahedrites. From the data of petrographic and 
mineralogical investigations of the melting crust and of meteoritic dust of the Sikhote-Alin meteorite, the 
authors V. D. Kolomenskii and 1. A. Iudin came to the conclusion that the crust has a double-layered structure 
and its outer layer mainly consists of magnetite. The meteoritic dust also consists of magnetite. 


Here the series of problems concerning this fall ended. 


The participants of the Conference watched with interest the scientific documentary film "The Sikhote- 
Alin’ meteorite.” 


The report of E. K. Gerling and L. K. Levskii on "The origin of inert gases in stony meteorites” aroused 
general interest. The isolation from some meteorites of the inert gases Ar, Ne and He and the determination of 
their isotopic composition led to the following conclusions: the large Ar” content in meteorites is connected 
with the radioactive decay of K*’; the ratio Ar’’/ Ar which was established and the isotopic composition of 
Ne in meteorites differ from those in air; the recrystallized chondrules contain the same isotopes of the inert 
gases as the original ones, which is found to be a result of the accumulation of these isotopes after the recrys- 
tallization of the meteoritic material caused by nuclear reactions with cosmic rays. The large quantity of 
inert gases, as compared with other meteorites, and the different isotopic composition of the meteorite "Staroe 
Pes*ianoe" indicate that the latter meteorite once belonged to a very large cosmic body. 


In the report of M. M. Shats on the work of I, E. Starik and M. M. Shats.,, who, by two parallel but dif- 
ferent methods, deterined the uranium content of meteorites and isotopic composition, he drew the conclusion 
that the uranium content in meteorites is generally lower than on earth, while the isotopic composition is the 
same. The data obtained by them showed again the constant occurrence of uranium in meteorites. 


In the lecture of K. P. Staniukovich and N. I. Shchetinina "On the problem of the inpact of a solid at 
high velocity on a compact medium" the character of shockwave propagation in a compact medium was con- 
sidered, as well as the determination of the velocity, density, and pressure dependent on the impact-center 


distance. 

Some considerations based on experimental material concerning the formations of chondrules from the 
gas phase, in a way similar to colloidal particles, were put forth in the paper of B. Iu. Levin and G. L. Slonimskii 
"On the problem of the origin of chondrules." 

In a report by A. Aaloe, the results of geological research and excavation of craters on the island of Saarema 
(Estonia) were given and these data fully confirmed the meteoritic nature of the craters. About 0.5 kg of 
meteoritic matter was collected. The study of the structural characteristics and the composition of meteoritic 
dust particles found in the region of the crater were of particular interest for the investigators. 


The paper by E. S. Burkser and K. N. Alekseeva was devoted to the physical properties of meteorites, to 
which little study has yet been given. The results of the study of some physical properties of stony meteorites 
showed that, as compared to the peridotites and pyroxenites which have about the same composition as stony 


145 


meteorites, the latter have a lower heat capacity than terrestrial crystallized rocks, a higher specific gravity and 
a higher electroconductivity, which can be explained by its nickel-iron content. On this content depend also the 
magnetic properties of meteorites. Stony meteorites have a relatively low melting point and high porosity, 
caused by their particular chondrulic structure. On the basis of the results of investigation of the above-men- 
tioned physical properties, an attempt was made to draw a picture of the origin of meteorites from asteroids. 


I. A. Iudin reported the results of petrographic investigations of some nontransparent minerals of primary 
and secondary origin which he had observed in stony meteorites, e.g., nickel-iron, non-nickel iron, native 
copper, troilite, iosite, ilmenite, chromite, Their structural features were studied in detail. 


In his paper concerning the petrographic investigations of the "Nikol'sk" stony meteorite which fell in 
1954, L. G. Kvashi characterized this meteorite as a typical chondrite, distinguished by the abundance of chon- 
drules of diameters varying from microscopical to 5 mm, and also by the special porosity of the basic mass, 
which contained about 10% metallic iron. 


The x-ray analysis, considered by V. I. Mikheev, is based on a definite dependence of the variation of 
crystal-lattice parameters on the mineral composition, which gives it a wide application for minute analysis 
of the composition of meteorites. 


F, Heide (German Democratic Republic) considered the origin of the craters of Carolina Bays. From 
aerial photographs and observations, F. Heide believed these craters to have been originated by the action of 
the sea on dry land. There are similarformsin the USSR near the town of Osipenko. 


G. S. Pokshivnitskii (Poland) reported on meteorites in Poland and on new findings of the iron meteorite 
"Morasko." 


The Seventh Conference on Meteorites concluded its work with a discussion of the papers and information 
given from the floor. 


The very wide range of problems connected with the appearance and falling of meteorites and with the 
study of their composition, clarified at the Conference, illustrates the great interest among the scientists both 
in meteorites and in cosmic material, the investigation of which will contribute to the solution of a series of 
important geochemical and cosmogonic problems. 


A display of literature on meteorites was arranged for the participants in the Conference. 
M. D'iakonova 
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